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This  report  emphasizes  results  from  completed  research  in  seismic  wave 
attenuation.  We  also  describe  and  illustrate  current  signal  analysis  methods 
using  a  variety  of  seismic  time  series  data. 


The  results  reported  include  relative  amplitudes,  at  a  large  number  of 
standardized  seismic  stations  in  the  U.S.  and  Southern  Canada,  for  a  large 
number  of  explosions  and  earthquakes  observed  at  three  fixed  azimuths  for 
many  stations.  The  average  at  1  Hz,  for  the  upper  300  km  of  the  crust- 
mantle  beneath  the  study  area  varies  from  about  65  at  Golden,  Colorado  to 
about  2000  at  Red  Lake,  Ontario.  The  mean  9a  is  140  for  the  continental  U.  S— 
These  absolute  9a  values  are  obtained  from  the  relative  amplitude  data  by 
requiring  that  the  largest  positive  amplitude  anomalies  for  the  U.S./Canada 
correspond  to  an  upper  mantle  Qa  with  a  value  greater  than  1000.  The  9„ 
values  obtained  correspond  to  t£  values  ranging  from  .02  up  to  .5,  with  t£  =  .3 
the  mean  value  for  the  continental  U.  S.  The  results  are  sufficiently  dense  spa¬ 
tially  to  allow  a  contour  map  for  Qa  to  be  constructed  which  incorporates  the 
azimuthal  variations  observed  at  many  of  the  stations.  This  representation 
shows  general  correlation  of  low  with  high  heat  flow  and  low  Pn  velocity  for 
example,  although  on  a  fine  spatial  scale  there  appears  to  be  numerous  excep¬ 
tions  which  may  be  related  to  unaccounted  for  focusing  and  scattering  effects 

rather than  anelastic  effects. 

'"Other  work  completed  has  provided  new,  high  resolution,  Q  determinations 
from  low  order  free  oscillation  excitation,  and  these  results  provide  strong  con¬ 
straints  on  lower  mantle  attenuation.^ 

^Tfe  describe  new  signal  analysis  methods  devised  for  use  in  automatic 
discrimination  and  yield  estimation  in  a  multivariant  approach.  In  this  context 
we  have  the  current  ability  to  rapidly  and  automatically  obtain  spectral  magni¬ 
tudes  for  a  large  number  of  seismic  phases,  which  can  be  used  in  both  telese- 
ismic  and  regional  discrimination  and  for  yield  estimation. 
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Research  on  Deterministic  Methods  of 
Seismic  Source  Identification 

Report  No.  4  -  Summary 

The  objectives  of  this  research  program  are  to:  (1)  Develop  and  perfect 
methods  of  synthesis  of  theoretical  seismic  radiation  fields  at  regional  and 
teleseismic  distance  ranges  in  order  to  provide  a  basis  for  developing,  testing 
and  understanding  seismic  source  discrimination  methods;  (2)  Develop  formal 
inversion  procedures  to  be  used  to  obtain  quantitative  estimates  of  physical 
parameters  of  seismic  sources  and  the  medium  and  to  use  these  estimates  to 
determine  underground  explosion  yield,  and  as  a  basis  for  source  identification 
and  discrimination;  (3)  Develop  the  means  and  methods  to  analyze  seismic  data 
for  source  identification  and  propagational  medium  characteristics  for  event 
discrimination  and  yield  estimation  purposes;  and  finally  to  (4)  Develop,  test 
and  evaluate  microprocessor  systems,  incorporating  software  and  hardware 
components  designed  to  provide  advanced  interactive  seismic  time  series 
analysis  for  RST  (Remote  Seismic  Terminal)  applications. 

In  the  first  year  annual  report  of  this  two  year  contract,  we  described 
research  results  in  the  areas  of  discrimination  using  inferred  physical  parame¬ 
ters  of  seismic  sources;  the  synthesis  of  complete  seismic  wave  fields  in  the 
near  and  regional  distance  ranges  from  the  seismic  source  and  the  utilization 
of  this  capability  in  source  and  path  of  propagation  studies;  basic  theoretical 
developments  in  seismic  source  representations  and  wave  propagation  studies 
end  the  implementation  and  application  of  advanced  signal  processing  methods. 

In  the  present  report  we  emphasize  results  relating  to  completed  research 
in  seismic  wave  attenuation  studies.  We  also  describe  and  illustrate  current  sig¬ 
nal  analysis  methods  using  a  variety  of  seismic  time  series. 
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Our  attention  to  studies  of  seismic  attenuation  is  motivated  by  the  strong 
spatial  variability  of  observed  body  wave  amplitudes  and  the  corresponding 
variability  and  uncertainties  when  seismic  body  wave  data  is  used.  Thus  in 
investigating  and  quantifying  the  magnitude  and  spatial  dependence  of  the 
adsorption  of  seismic  energy  in  the  earth,  we  expect  to  be  able  to  infer  the 
mechanisms  responsible  for  the  energy  loses  and  their  relationships  to  other 
geophysical  variables,  such  as  heat  flow,  seismic  velocities,  seismicity  and 
tectonic/geologic  setting,  and  thereby  be  able  to  infer  attenuation  properties 
of  particular  regions  on  the  earth  based  on  other,  more  easily  measured  and/or 
currently  known  seismic/geophysical  variables.  More  directly,  however,  we  wish 
to  quantify  the  magnitude  and  variability  of  attenuation  losses  in  the  context  of 
explosive  yield  estimation  and  also  for  discrimination,  particularly  for  applica¬ 
tions  of  spectral  discriminants. 

The  results  reported  here  describe  the  analysis  of  body  wave  amplitude 
data,  from  a  large  number  of  explosions  and  earthquakes,  observed  at  WWSSN, 
LRSM  and  other  standardized  seismic  stations  throughout  the  continental  U.  S. 
and  Southern  Canada.  The  average  Qa  (dissipation  function  for  P  waves)  in  the 
upper  300  km  of  the  mantle  beneath  the  U.S.  is  140  at  1  Hz.  This  mean  value  is 
obtained  from  the  relative  amplitude  values  at  the  stations  in  the  U.S.  and  So. 
Canada,  under  the  condition  that  the  largest  relative  amplitudes  observed  (at 
the  station  at  Red  lake  Ontario)  correspond  to  an  upper  mantle  Qa  of  about 
2000.  In  this  case  the  variation  in  the  upper  mantle  ranges  from  above  1000 
to  about  65,  the  latter  at  Golden,  Colorado  and  the  former  at  Red  Lake  Ontario. 
These  values  of  Qa  correspond  to  an  average  t£  =  .3,  for  one  way  transit  of  a  P 
wave  through  the  upper  mantle,  and  variations  from  approximately  =  .5  to 
tS  =  .02.  The  station  sampled  variations  of  the  mean  upper  mantle  and  crust 
Qm  throughout  the  U.  S.  are  used  to  construct  a  contour  map  which  shows  the 
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variation  of  the  mean  Q *  in  the  U.S.,  and  its  correlation  with  geophysical  pro¬ 
vinces  and  independent  geophysical  observations,  such  as  heat  flow,  gravity  and 
Pn  velocities.  There  is  broad  correlation  of  low  Qa  with  high  seismicity  and  high 
heat  flow  and  with  low  Pn  velocity  and  gravity,  although  on  a  fine  spatial  scale 
there  appear  to  be  numerous  exceptions.  (i.e.  local  anomalous  Qa  values,  of 
small  dimension,  not  correlating  with  known  average  values  of  seismicity,  heat 
flow,  etc.  These  may  be  due  to  focusing  and/or  scattering  effects  in  the  upper 
mantle  or  crust,  rather  than  anelastic  effects.) 

In  other  reported  work  on  anelasticity  and  associated  seismic  wave 
attenuation  we  have  used  the  free  oscillations  of  the  earth  to  investigate  the  Q 
for  the  lower  mantle.  In  this  work  we  show  that  the  Q  for  0S2  and  0S3  are  399 
and  358,  respectively.  The  accurate  Q  estimates  for  low  order  free  oscillations 
provide  us  with  the  means  of  determining  the  proper  attenuation  corrections 
for  P  wave  propagation  in  the  lower  mantle,  to  be  used  to  obtain  accurate  yield 
estimates  from  teleseismic  data.  While  the  lower  mantle  absorption  is  relatively 
small,  the  knowledge  of  the  lower  mantle  Q  allows  us  to  more  accurately  assess 
the  large  and  variable  upper  mantle  effects. 

We  describe  new  signal  analysis  methods  devised  for  use  in  automatic 
discrimination  and  yield  estimation  in  a  multivariant  approach.  In  this  context, 
we  have  the  current  ability  to  rapidly  and  'automatically'  obtain  spectral  mag¬ 
nitudes  for  a  large  number  of  seismic  phases,  which  can  be  used  in  both  telese¬ 
ismic  and  regional  discrimination,  such  as  frequency  dependent  body  and  sur¬ 
face  wave  magnitudes.  We  expect  to  implement  some  of  these  methods  at  the 
Center  for  Seismic  Studies  and  on  a  remote  seismic  terminal  (RST  system)  for 
extensive  and  systematic  testing  in  the  near  future. 


1.  Introduction 


Hie  objectives  of  the  research  being  conducted  under  the  current  2  year 
contract  are  to:  (l)  Develop  and  apply  methods  of  seismogram  synthesis  using 
mode  superposition  and  similar  methods  as  aids  in  discrimination/yield  estima¬ 
tion  studies;  (2)  Develop  and  apply  source  inversion  methods,  to  infer  physical 
properties  of  seismic  sources  for  discrimination  purposes;  (3)  Develop  regional 
discrimination  techniques  based  on  the  physical  properties  of  seismic  sources; 
(4)  Develop/apply  methods  of  determining  anelastic,  frequency  dependent, 
absorption  characteristics  of  the  medium;  (5)  Pursue  studies  designed  to  pro¬ 
vide  a  physical  basis  for  seismic  discriminants;  (6)  Apply  and  evaluate 
automated  methods  of  signal  analysis  for  event  detection  and  discrimination  in 
the  regional  and  teleseismic  distance  ranges;  (7)  Pursue  studies  of  wave  propa¬ 
gation  along  oceanic  paths  with  the  objectives  of  understanding  the  propaga¬ 
tion  out  to  large  distances  and  developing  discrimination  and  yield  estimation 
methods  utilizing  the  broad  band  nature  of  the  signals  observed;  and  (B)  Imple¬ 
ment  advanced  time  series  analysis  software  on,  both,  large  main-frame  and  the 
RST  (Regional  Seismic  Terminal)  microprocessor.  Test  and  evaluate  the  main 
frame  operation  and  the  RST  system,  and  recommend  final  hardware/software 
systems  for  general  and  specialized  use. 

In  the  previous  annual  report,  for  the  period  October  1,  19B1  to  September 
30,  1982,  we  summarized  results  relating  to:  (1)  the  use  of  inferred  physical 
parameters  of  seismic  sources  (stress  changes  and  rupture  zone  dimensions)  as 
a  means  of  event  discrimination;  (2)  modal  synthesis  of  the  complete  seismic 
wave  Field  in  the  near  and  regional  distance  ranges  for  source  and  path  studies; 
(3)  theoretical  developments  in  seismic  source  representations  and  wave  propa¬ 
gation  studies;  (4)  applications  of  advanced  signal  processing  methods  to  high 
frequency  near  field  data  for  inferences  of  source  properties;  and  (5)  impie- 
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mentation  of  signal  processing/analysis  software  on  an  RST  microprocessor  and 
on  a  CRAY- 1  main  frame  computer. 

In  the  present  report  we  emphasize  research  results  relating  to  completed 
research  in  seismic  wave  attenuation  studies  (Sections  II  and  III).  We  also 
include  a  summary  of  the  signal  processing  capability  that  has  been  developed, 
along  with  examples  of  its  application  in  Section  IV. 

The  importance  of  the  attenuation  studies  relates  to  the  strong  depen¬ 
dence  of  yield  estimates  on  the  attenuation  characteristics  of  the  upper  mantle 
and  crust,  particularly  at  the  source  and  receiver  sites.  In  addition,  spectral 
discriminants,  such  as  the  VFM  (Variable  Frequency  Magnitude  or  m6(/)) 
discriminant,  are  strongly  influenced  by  local  variations  in  attenuation.  For 
applications  of  this  discriminant  then,  it  is  important  to  be  able  to  isolate  the 
zones  of  high  attenuation  in  terms  of  the  magnitude  of  the  attenuation,  the 
lateral  extent  of  such  zones  and,  as  importantly,  the  depth  distribution  of  the 
high  attenuation  zones  (upper  mantle  versus  crust  in  particular).  Specifically, 
if  high  attenuation  zones  are  confined  to  the  upper  mantle,  world-wide,  then 
regional  discrimination  using  VFM  methods  should  be  extremely  effective  v-hen 
(only)  crustal  phases  are  used.  Finally,  the  frequency  dependence  of  the 
attenuating  mechanism,  as  well  as  its  depth  distribution,  is  of  considerable 
importance  for  both  yield  estimation  and  discrimination  since  high  frequency 
spectral  magnitudes  may  be  much  less  affected  by  attenuation  processes  in  the 
earth  and  therefore  more  reliable  for  both  yield  estimation  and  discrimination. 

The  emphasis  placed  on  the  advanced  signal  analysis  methods  is  based  on 
the  requirement  that  fast  and  systematic  methods  of  considerable  sophistica¬ 
tion  are  needed  to  implement  yield  estimation  as  well  as  discrimination  pro¬ 
cedures  in  practice.  Our  results  in  this  area  are  particularly  encouraging  and 
we  hope  to  implement  the  methods  developed  on  the  new  RST  system,  cs  well  as 
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at  the  Center  for  Seismic  Studies,  in  the  near  future. 

II.  Attenuation  of  Seismic  P  waves  and  Amplitude  Variations  in  the 

Continental  U.  S. 

This  study  (Butler,  1983),  is  largely  concerned  with  amplitude  effects  that 

may  be  ascribed  to  the  locale  of  the  seismic  station.  Naturally  to  obtain  this 

end.  a  great  deal  of  care  must  be  taken  to  minimize  amplitude  effects  which  are 

attributable  to  the  seismic  source  or  the  propagation  path.  The  source  term 

for  each  event  used  should  be  nearly  the  same  at  all  stations  of  the  network  - 

this  constraint  is  relatively  easy  for  explosions,  but  can  be  quite  difficult  for 

earthquakes.  A  consistent  amplitude  measure  from  the  P  waveform  must  be 

made  at  all  stations,  and  it  is  favorable  to  have  one  seismologist  make  all  the 

measurements.  The  distance  range  between  sources  and  receivers  should  be 

such  that  the  seismic  waves  bottom  in  the  relative  smooth  lower  mantle  of  the 

earth,  avoiding  upper  mantle  triplications  and  core-mantle  boundary  diffrac- 

p  o 

tion.  This  corresponds  to  a  distance  range  of  30  to  95  ~  3300  to  10500  km 

For  S  waves,  interference  with  SKS  and  ScS  beyond  B5*  must  be  avoided.  The 
amplitude  data  are  corrected  for  the  seismometer  gain  and  corrected  for 
geometric  spreading  to  a  reference  distance  -  the  mean  distance  of  the  data  set 
is  reasonable.  The  amplitude  data  should  be  processed  separately  for  different 
azimuths,  as  azimuthally  dependent  station  amplitudes  have  been  observed  for 
WWSSN  stations  in  the  United  States  (Butler  and  Ruff,  1980)  and  at  the  large 
seismic  arrays  LASA  (Chang  and  von  Seggern,  1980)  and  N0RSAR  (Berteussen, 
1975). 

A  map  of  the  United  Stales  in  Figure  1  shows  the  location  of  WWSSN  sta¬ 
tions  in  relation  to  physiographic  provinces.  Table  1  lists  the  names  and  loca¬ 
tions  of  the  WWSSN  stations.  The  locations  of  three  digital  SDCS  (Special  Data 
Collection  System)  stations  are  also  noted  -  RK0N,  Red  Lake,  Ontario,  HNME, 
Houlton,  Maine;  and  0B2NV,  Oak  Springs  Butte,  Nevada.  These  digital  SDCS 
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Table  1 

Seismic  Stations 


Station  Code 

Name  and  Location 

AAM 

Ann  Arbor.  Michigan 

ALQ 

Albuquerque,  New  Mexico 

ATL 

Atlanta,  Georgia 

BEC 

Bermuda 

BKS 

Berkeley.  California 

BLA 

Blacksburg,  Virginia 

BOZ 

Bozemann,  Montana 

COR 

Corvallis.  Oregon 

DAL 

Dallas,  Texas 

DUG 

Dugway,  Utah 

FLO 

Florissant,  Missouri 

GEO 

Georgetown.  Washington.  D.  C. 

GOL 

Golden,  Colorado 

GSC 

Goldstone,  California 

HNME 

Houlton,  Maine 

JCT 

Junction  City.  Texas 

LON 

Longmire,  Washington 

LUB 

Lubbock,  Texas 

MDS 

Madison,  Wisconsin 

MSO 

Missoula,  Montana 

0B2NV 

Oak  Springs  Butte,  Nevada 

OGD 

Ogdensburg,  New  Jersey 

OXF 

Oxford,  Mississippi 

RCD 

Rapid  City,  South  Dakota 

RKON 

Red  Lake,  Ontario 

SCP 

State  College,  Pennsylvania 

SKA 

Spring  Kill.  Alabama 

TUC 

Tucson,  Arizona 

WES 

Weston,  Massacahuselts 

B 


as 
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stations  are  filtered  to  match  the  WWSSN  short  period  instrument  response  and 
may  be  used  to  tie  LRSM  studies  to  the  WWSSN.  The  WWSSN  stations  consist  of 
■ix  seismometers  -  a  vertical  and  two  horizontal  instruments  with  matched 
responses,  recording  short  period  (1  sec)  and  long  period  (20  sec)  ground 
motion.  The  SDCS  sites  are  also  equipped  with  three  component  short  period 
and  long  period  seismometers.  In  this  study  all  amplitudes  are  measured  from 
the  short  period  vertical  instrument.  Most  of  these  stations  are  located  on  hard 
rock  sites.  Stations  located  on  sediments  on  the  west  coast  and  in  the  central 
United  States  -  BKS,  RCD,  LUB,  JCT,  DAL,  OXF,  SHA,  AAM,  FLO  -  are  reviewed  in  a 
following  section  for  signal  amplification  effects  due  to  the  sediments. 

The  amplitude  measure  used  in  this  study  is  the  peak  to  peak  amplitude  of 
the  first  cycle  of  P  wave  motion.  This  measure  was  found  to  be  more  stable  in 
the  presence  of  seismic  noise  than  a  peak  amplitude  measurement  of  the  P 
wave  first  motion.  Crustal  reverberations,  near  source  surface  reflections,  and 
P  to  S  conversions  near  the  receiver  may  contaminate  the  amplitude  through 
constructive  or  destructive  interference,  but  the  first  cycle  of  P  motion  is  most 
free  of  these  difficulties. 

Earthquake  and  underground  nuclear  explosion  sources  are  used  in  this 
amplitude  study.  Underground  nuclear  explosions  have  nearly  the  same 
waveform  at  all  stations  with  an  impulsive,  distinct  shape  allowing  an  unambigu¬ 
ous  measurement  of  amplitude.  The  earthquake  data  is  not  so  simple  and  con¬ 
siderable  care  was  used  in  event  selection  for  this  study. 


-10- 


Specifically,  amplitude  data  from  earthquake  sources  was  selected  by  the 
set  of  criteria  of  Butler  (I979a,b),  Butler  and  Ruff  (1980),  and  Lay  and  Helm- 
berger  (19B1).  The  following  procedure  is  used:  Earthquakes  in  South  America 
and  the  Asian-Pacific-Aleutian  arc  were  viewed  on  seismograms  from  three 
widely  spaced  stations  in  the  United  States.  Events  whose  waveforms  passed 
the  set  of  criteria  at  the  three  stations  were  examined  at  all  stations.  Out  of 
the  set  of  events  which  fulfilled  the  selection  criteria  at  all  stations,  a  quality 
subset  of  events  were  measured  for  amplitude.  Even  under  these  stringent 
selection  procedures,  some  P  waveform  vagaries  remained  in  the  data  set.  In 
some  instances,  one  station  would  have  a  slightly  different  frequency  content 
than  other  stations.  In  other  instances  the  waveform  at  a  station  might  have  an 
additional  inflection,  or  the  relationship  between  peak  amplitudes  might  be 
slightly  different,  or  the  signal  to  noise  ratio  at  a  station  would  be  borderline. 
These  vagaries  were  handled  in  the  following  manner.  An  amplitude  w as  meas¬ 
ured.  but  a  note  of  the  questionable  behavior  was  made.  After  all  events  bad 
been  measured,  reduced,  processed,  and  preliminary  station  amplitude  effects 
determined,  then  the  amplitude  data  were  reviewed.  Amplitudes  which  were 
greater  than  two  standard  errors  away  from  the  mean  were  checked  to  see  if  a 
note  of  questionable  behavior  had  been  made.  If  a  note  had  been  made,  the 
amplitude  measure  was  deleted  from  the  jata  set.  If  a  no  questionable  behavior 
had  been  noted,  the  amplitude  measure  was  retained  in  the  data  set.  Following 
this  method,  four  P  wave  amplitude  measurements  from  waveforms  with  ques¬ 
tionable  behavior  were  deleted. 
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Three  azimuths  of  approach  were  examined  in  determining  the  station 
amplitude  effects  for  the  WWSSN  stations  in  the  United  States.  For  a  northern 
azimuth,  nuclear  explosion  sources  in  the  U.S.S.R.  were  analyzed.  Most  of  these 
data  were  discussed  in  Butler  and  Ruff  (1980).  The  explosion  data  set  were  dis¬ 
tributed  among  five  test  sites:  northern  and  southern  Novaya  Zemlya,  eastern 
and  western  Semipalatinsk,  and  western  Kazakh.  Three  additional  explosions  in 
the  eastern  Semipalatinsk  site  in  the  period  1976-1977  were  measured  at  the 
WWSSN  stations  and  at  the  three  SDCS  Stations,  RKON,  HNME,  and  0B2NV  to  tie 
the  SDCS  stations  to  the  W'WSSN  network.  A  pattern  of  relative  amplitudes  in 
the  United  States  was  determined  for  each  test  site,  and  the  results  from  the 
five  test  sites  were  averaged  to  yield  a  relative  amplitude  pattern  for  a  northern 
azimuth  from  the  United  States.  The  apparent  period  of  the  first  cycle  of  P 
wave  motion  for  the  explosion  sources  measured  on  the  WWSSN  short  period 
instrument  was  usually  about  1  second,  with  little  variation  from  one  event  to 
another  (a  range  of  1.2  to  0.8  sec).  A  special  exception  was  the  SDCS  station  • 
filtered  to  have  a  WWSSN  response  •  RKON,  which  had  an  apparent  period  of  0.5 
seconds  for  the  first  cycle  of  P  wave  motion  from  three  Semipalatinsk  explo¬ 
sions.  This  anomalous  observation  will  be  discussed  later.  Aside  from  the 
strange  behavior  of  RKON,  the  amplitude  data  from  explosions  is  quite  stable, 
and  for  a  single  test  site  the  relative  amplitude  pattern  shows  little  variation 
from  event  to  event  (Butler  and  Ruff,  1980).  Amplitude  pattern  variation 
between  the  test  sites  indicating  near  source  heterogeneity  is  discussed  in 
Butler  and  Ruff  (1980).  Core-mantle  boundary  diffraction  effects  at  distances 
greater  than  95.5*A  are  discussed  in  Ruff  and  Helmberger  (1981).  Amplitude 
data  which  suffered  diffraction  at  the  core-mantle  boundary  were  excluded  in 
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this  present  study  of  receiver  amplitude  effects.  ' 

The  location  and  origin  times  of  the  36  nuclear  explosions  analyzed  are 
listed  in  Table  2  and  grouped  by  their  respective  test  sites.  The  amplitude  data 
measured  from  the  first  cycle  of  P  wave  motion  are  listed  in  the  Appendix  in 
Tables  A1-A6.  The  data  have  been  corrected  for  seismometer  gain  and  are 
expressed  in  nanometers  of  ground  displacement.  No  geometric  spreading 
corrections  were  made  as  the  effect  is  ±10%  for  this  geometry  of  sources  and 
receivers.  However,  in  transforming  the  amplitude  variations  into  apparent 
variations  of  the  dissipation  parameter  Q  or  attenuation,  geometric  spreading 
corrections  are  included. 

While  the  relative  amplitude  for  a  northern  azimuth  from  the  United  States 
comes  exclusively  from  explosions,  the  relative  amplitude  patterns  of  the 
WWSSS  stations  in  the  United  States  for  both  a  northwest  azimuth  and  a  south- 
southeast  azimuth  come  exclusively  from  earthquakes.  Events  at  a  northwest 
azimuth  include  earthquakes  in  the  Bonin  Islands,  Japan,  Kurile  Islands,  and  the 
Aleutians.  Events  at  a  south-southeast  azimuth  are  comprised  by  earthquakes 
in  South  America  where  most  of  the  events  lie  along  the  north  and  central 
Chile-Argentine  border.  Amplitude  data  from  some  of  these  events  were 
presented  in  a  preliminary  form  by  Butler  and  Ruff  (I960),  who  were  largely 
concerned  with  establishing  to  what  extent  were  azimuthal  variations  of  station 
amplitude  effects  important.  Finding  that  azimuthal  variation  is  relatively 
important,  this  present  study  follows  to  consider  these  effects  in  more  detail. 
The  earthquakes  were  not  sufficiently  localized  to  be  divided  into  subgroups  in 
the  fashion  that  the  explosions  were  divided  among  test  sites.  However,  though 
the  earthquakes  were  subdivided  among  broad  source  regions  to  look  for  sys¬ 
tematic  source  region  effects  such  as  focusing/defocusing  by  subduction  zone 
structure,  in  the  final  analysis  the  amplitude  data  from  the  earthquakes  are 
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Table  2 

Explosion  Data  Set* 


Northern  Novaya  Zemlya 


27 

Oct 

66 

5:57:58 

73.44N 

54.75E 

21 

Oct 

67 

4:59:58 

73.37N 

54.B1E 

7 

Nov 

68 

10:02:05 

73.40N 

54.B6E 

14 

Oct 

69 

7:00:06 

73.40N 

54.B1E 

14 

Oct 

70 

6:02:57 

73.31N 

55.15E 

27 

Sept 

71 

5:59:55 

73.39N 

55.10E 

28 

Aug 

72 

5:59:57 

73.34N 

55.08E 

12 

Sept 

73 

6:59:54 

73.30N 

55.16E 

29 

Aug 

74 

9:59:56 

73.37N 

55.09E 

23 

Aug 

75 

8:59:58 

73.37N 

54.64E 

21 

Oct 

75 

11:59:57 

73.35N 

55.09E 

Southern  Novaya  Zemlya 


27  Oct 

73 

6:59:58 

70.76N 

53.97E 

2  Nov 

74 

4:59:57 

70.B2N 

54.06E 

18  Oct 

75 

8:59:56 

70.B4N 

53.69E 

Semipalatinsk  East 

15  Jan 

65 

5:59:59 

49.89N 

7B.97E 

30  Nov 

69 

3:32:57 

49.92N 

79.00E 

2  Nov 

72 

1:26:58 

49.9  IN 

7B.94E 

23  Jul 

73 

1:22:58 

49.99N 

7B.B5E 

14  Dec 

73 

7:46:57 

50.04N 

79.01E 

31  May 

74 

3:26:57 

49.95N 

78.94E 

4  Jul 

76 

2:56:56 

49.91N 

7B.95E 

Semipalatinsk  East-Additional  Data 


23  Nov  76 
7  Dec  76 
29  May  77 


5:03:00 

4:57:00 

2:59:00 


50.00N 

49.90N 

49.90N 


79.00E 

78.90E 

79.90E 
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Table  2 

Explosion  Data  Set* 


Semipalatinsk  Vest 


19 

Oct 

66 

3:57:58 

49.75N 

78.03E 

20 

Apr 

87 

4:07:58 

49.74N 

78.12E 

17 

Oct 

87 

5:03:58 

49.82N 

7B.10E 

29 

Sept 

68 

3:42:58 

49.77N 

7B.19E 

28 

Jun 

70 

1:57:58 

49.83N 

78.25E 

22 

Mar 

71 

4:32:58 

49.74N 

7B.19E 

25 

Apr 

71 

3:32:5B 

49.B2N 

7B.09E 

30 

Dec 

71 

6:20:58 

49.75N 

7B.13E 

20 

Feb 

75 

5:32:58 

49.82N 

72.08E 

Kazakh 


6 

Dec 

69 

7:02:57 

43.B3N 

54.7BE 

12 

Dec 

70 

7:00:57 

43.B5N 

54.77E 

23 

Dec 

70 

7:00:57 

43.B3N 

54.B5E 

'Locations  and  origin  times  (1966-1978)  from  Dahlman  and  Israelson  (1977), 
after  July.  1976  from  Seismic  Data  Analysis  Center 


APRIL  25,  1967  SOUTH  AMERICAN  EARTHQUAKE 


Figure  2.  P  waves  recorded  on  short  period  vertical  seismometers  of  WW’SSN 
stations  In  the  United  States  are  shown  for  an  earthquake  in 
South  America  on  April  25,  1967.  The  amplitude  of  the  first 
cycle  of  P  wave  motion  is  given  in  nanometers,  corrected  for 
geometric  spreading  to  a  distance  A  ■  60*. 
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Table  3 

Earthquakes  in  a  Northwest  Azimuth 


Date 

Origin  Time 

Location* 

Depth  (km) 

Kurile  Islands  Earthquakes 

1 1/22/66 

6:29:52.4 

48.0  N  146.BE 

441 

3/20/87 

12:31:34.0 

45. 6N  151.4E 

51 

B/ 10/67 

11:21:22.3 

45.4N  150.3E 

37 

2/10/68 

10:00:05.8 

46.  ON  152.3E 

67 

4/28/68 

4:18:15.7 

44. BN  174.5E 

39 

7/25/68 

10:50:31.5 

45.7N  146.7E 

16 

10/26/76 

5:58:56 

46.  IN  159. 7W 

130 

3/19/77 

10:56:06 

43.  ON  149.0E 

0 

Japanese  Earthquakes 

12/31/78 

9:16:37 

40.0N  145.0E 

0 

1/1/77 

11:33:42 

30.  BN  137.2E 

483 

2/18/77 

20:51:26 

34.  ON  143.0E 

0 

6/12/77 

8:48:05 

4  3.  ON  142. 3E 

241 

• 

Bonin  Islands  Earthquakes 

9/22/76 

8:20:28 

23.3N  142. IE 

110 

12/5/76 

22:01:22 

23.0N  140. OE 

393 

12/22/76 

1:01:42 

24.0N  145. OE 

0 

1/5/77 

22:44:57 

23.3N  143.BE 

0 

Date 

Origin  Time 

Location* 

Depth  (km) 

Region 

Other  Events 

10/22/76 

18:35:24 

56.  IN  153.31V 

0 

Kodiak  Islands 

4/22/77 

5:58:56 

52.5N  13B.BE 

408 

Kamchatka 

4/23/77 

0:52:05 

75.  ON  134. 9E 

0 

New  Siberian  Is. 

7/20/77 

14:49:06 

50. 6N  161.9V 

0 

Alaska  Pen. 

6/7/77 

23:26:55 

52.2N  176.21V 

125 

Andreanof  Is. 

'Locations  before  1975  from  International  Seismological  Center,  after  1975 
from  Seismic  Data  Analysis  Center 


-17- 


r? 


Table  4 

South  American  Earthquakes* 


Date 

Origin  Time 

Location 

Depth  (km) 

4/25/67 

10:26:14.3 

32.8S  69. 0W 

39 

11/15/67 

21:35.51.5 

2B.7S  71.21V 

15 

2/6/67 

11:19:23.1 

28. 5S  71. 0W 

23 

4/21/6B 

9:24:35.5 

23.4S  70.51V 

41 

4/30/88 

23:51:17.9 

3B.4S71.1W 

40 

9/30/76 

8:04:11 

24.2S  69.21V 

0 

12/3/76 

5:27:34 

21. OS  69.01V 

79 

12/4/76 

12:32:35 

20.0S  69.0W 

103 

3/8/77 

22:46:44 

B.OS  63.01V 

0 

3/13/77 

4:55:55 

2.0S  58. OW 

0 

4/15/77 

23:35:38 

22.9S  68.B1V 

109 

6/  2/77 

16:50:36 

29.9S  88. 61V 

94 

6/  5/77 

2:46:07 

24. OS  70.51V 

30 

6/  B/77 

13:25:16 

22. IS  67.31V 

135 

6/18/77 

16:49:42 

21. OS  69.71V 

125 

'Locations  before  1975  from  International 
Seismological  Centre,  after  1975  from  Seismic 
Data  Analysis  Center 
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divided  only  between  northwest  and  south-southeast  azimuths.  The  earthquake 
sources  had  slightly  longer  apparent  periods  of  the  first  cycle  of  P  wave  motion 
than  were  observed  for  the  explosion  sources.  The  apparent  period  of  the  first 
P  wave  cycle  for  individual  events  ranged  between  1.6  and  1.0  seconds,  with  a 
mode  of  1.2  seconds.  The  individual  stations  exhibited  a  range  of  variation  of 
apparent  period  of  ±0.2  seconds  about  the  average  apparent  period  for  the 
whole  network.  Stations  exhibiting  a  greater  frequency  variation  were  either 
noted  as  such,  or  else  the  event  itself  was  discarded.  While  the  individual  sta¬ 
tions  occasionally  exhibited  slightly  longer  or  shorter  apparent  periods  than 
the  network  average,  no  station  other  than  RKON  showed  a  consistent,  sys¬ 
tematic  deviation  of  frequency  content.  RKON  at  Red  Lake.  Ontario,  in  the 
Canadian  shield  often  but  not  always  was  shorter  period  than  other  stations. 

The  origin  times  and  locations  of  the  earthquakes  analyzed  are  presented 
in  Tables  3  and  4.  Twenty-one  earthquakes  to  a  northwest  ezimuth  from  the 
United  States  are  listed  in  Table  3.  and  are  subdivided  among  four  groups: 
Kurile  Islands,  Japan,  Bonin  Islands,  and  other  events.  Fifteen  earthquakes  in 
South  America  are  listed  in  Table  4  and  are  not  subdivided.  The  amplitude  data 
obtained  from  these  events  are  presented  in  the  Appendix  in  Tables  A  7  through 
All.  The  amplitude  data  have  been  corrected  for  seismometer  gain  and 
corrected  for  geometric  spreading  to  a  common  distance  of  C0d£,  and  ere 
expressed  as  nanometers  of  ground  displacement.  The  P  waveforms  of  4/25/67 
South  American  earthquake  are  plotted  in  Figure  2.  While  the  wave  shape  of 
the  P  wave  maintains  a  relative  simplicity  and  constancy  over  the  whole  WWSSN 
array  in  the  United  States,  the  absolute  amplitudes  vary  by  nearly  a  factor  of 
twenty  between  LON,  Longmire,  Washington,  and  OXF,  Oxford,  Mississippi.  Low 
amplitudes  are  seen  at  DUG.  Dugway,  Utah:  TUC,  Tucson,  Arizona;  and  GOL,  Gol¬ 
den,  Colorado.  High  amplitudes  are  observed  at  I.UB,  Lubbock,  Texas;  DAL,  Dal- 


las,  Texas;  JCT,  Junction  City.  Texas;  and  ATL,  Atlanta.  Georgia.  Even  though  the 
data  shown  for  the  4/15/67  event  in  Figure  2  represent  a  paradigm  of  ampli¬ 
tude  variation  for  a  simple  earthquake,  these  amplitude  data  must  be  averaged 
with  many  other  events  to  properly  yield  amplitude  effects  which  may  be  attri¬ 
butable  to  the  region  of  the  receiving  stations.  Waveforms  for  many  of  the 
events  in  Tables  3  and  4  are  plotted  in  a  report  by  Butler  and  Hart  (1979). 

Amplitude  Data  Reduction 

The  absolute  amplitudes  of  the  events  in  Appendix  Tables  A1  through  All 
measured  at  the  WWSSK  and  filtered  SDCS  stations  in  the  United  States  range  in 
size  from  tens  of  nanometers  to  thousands  of  nanometers.  While  the  explosions 
and  earthquakes  in  Tables  2,  3,  and  4  may  have  minimal  amplitude  variations 
across  our  network  of  seismic  stations  which  are  due  to  source  effects,  the 
events  have  different  scale  sizes;  i.e..  some  events  are  bigger  or  smaller  than 
others.  Let  us  presume  that  the  scale  size  of  an  event  may  be  characterized  by 
a  constant  a.  The  body  wave  magnitude,  mb ,  is  such  a  constant  and  is  deter¬ 
mined  from  averaging  mb  measurements  made  at  stations  all  over  the  world.  In 
our  case  we  would  like  to  determine  the  scale  size  constants  dj  for  i  events 
directly  from  our  own  amplitude  data  set.  If  we  divide  our  absolute  amplitude 
measurements  for  event  i  by  the  scale  size  constant  oq.  we  have  then  effec¬ 
tively  removed  the  size  of  the  source  and  are  left  with  a  residual  variation  from 
station  to  station  which  is  dimensionless  and  expresses  the  relative  amplitudes 
of  the  stations.  The  relative  amplitudes  at  a  single  station  over  many  events  is 
a  distribution  which  may  be  characterized  by  the  mean  and  the  variance.  The 
mean  relative  station  amplitude  characterizes  the  amplitude  variation  which 
may  be  attributable  to  the  locale  of  the  station.  If  the  event  scale  size  con¬ 
stants  at  are  larger  or  smaller  than  necessary,  the  total  variance  of  the  relative 
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amplitudes  summed  over  all  stations  will  increase.  A  simple  method  of  obtain¬ 
ing  event  scale  size  constants  at  is  to  let  at  be  the  average  amplitude  measured 
for  event  i.  However,  each  station  does  not  yield  an  amplitude  measurement  on 
every  occasion  -  reason:  poor  signal  to  noise,  no  record,  station  down,  ampli¬ 
tude  off  scale,  etc  -  and  thus  the  stations  used  in  the  network  average  would 


vary  from  one  event  to  another.  A  more  efficient  method  is  to  determine  the  a, 
while  simultaneously  assuring  minimum  total  variance  of  the  remaining  relative 
amplitudes.  An  iterative  procedure  was  adopted  to  obtain  this  end. 

Let  Ay  be  amplitude  measurements  from  the  i  events  at  the  j  stations  for 
a  given  azimuth  for  the  earthquake  sources  or  for  a  given  test  site  for  the 
explosion  sources.  The  Ay  have  already  been  corrected  for  seismometer  gain 
and  geometric  spreading.  From  the  t  events  choose  a  master  event  k  and 
divide  the  A*j  by  the  master  event  average  value.  pk 

**  (2) 

where  nk  is  the  number  of  stations  reporting  the  At,  amplitude  observations, 
n*  <  j,  for  event  k.  Now  for  events  i  Jt  k  determine  scale  factors  a,  such  that 
for  each  i 


min  £ 
i 


Oj  Ay  — 


(3) 


is  minimum.  Using  these  scale  factors  a,1  (where  superscript  k  simply  denotes 
reference  to  master  event  k)  the  squared  variation  of  i  events  are  minimized 
relative  to  the  master  event  k.  The  total  amount  of  squared  variation  is 


t4=  E  E 

*  i 


(4) 


As  not  all  stations  record  amplitude  data  from  each  event,  the  total  squared 
variation  14  for  master  event  k  is  scaled  by  the  number  of  elements  m*  in  the 
sum  expression  in  (4) 

V*  =  “  Uk  (5) 


If  the  unsealed  variation  L'k  were  not  rescaled  to  V*.  then  the  master  event  k 
which  had  the  fewest  observations  would  tend  to  produce  the  minimum  total 


-22- 


variation;  since  if  event  k  bas  no  observation  for  station  j,  observations  at  sta¬ 
tion  J  for  events  i  cannot  be  used  in  the  sums  in  (3)  and  (4).  Next  iterate  the 
master  event  k  over  the  events  t  and  choose  the  master  event  U  which  minim¬ 
izes  Vk 

Vu  ~  min(  Vk  )xmk  #vtr  ^  j  (6) 

We  now  form  scaled  observations  C(Jf 

Ofj  *  V-  Akj  ,k  -  U  (7) 

0%  =  Axj  ,iP  U  (S) 

At  is  defined  in  (2),  a»  is  determined  from  condition  (3).  the  superscript  M  on 
Q*j  denotes  an  observation  rescaled  with  the  master  event  k  =  M.  Each  station 
j  now  has  P}  observations  from  the  t  events,  Pj  si.  The  mean  relative  ampli¬ 
tude  for  each  station  may  be  calculated 

<9> 

r)  i 

Note  that  the  0%  were  calculated  relative  to  the  master  event  M  which  pro¬ 
duced  a  minimum  average  variance  V*.  The  estimate  of  mean  relative  station 
amplitude  pt;  can  be  improved  by  substituting  our  preliminary  estimate  of  - 

i.e..  nf  •  into  (3)  for  Akj 

min  £  [a,**  ^  -  ^j'f]  0°) 

i  4 

and  solving  for  a/*  which  minimizes  (10).  New  scaled  amplitude  observations  are 

formed 

0$  -  o.\  Ai}  .for  all  i,j  (11) 

The  improved  estimate  of  the  mean  relative  station  amplitude  is 


-23- 


=  IT  £  (12) 

'  i  » 

where  Pj  is  the  number  of  events  i  for  which  the  amplitude  is  measured  for  sta¬ 
tion  j,  Pj£i. 

For  this  amplitude  data  set.  it  was  found  that  the  estimator  ^  is  stable; 
such  that  further  iterations  are  unnecessary.  The  sample  standard  deviation  of 
the  relative  station  amplitudes  is  then 


and  is  defined  when  there  is  more  than  one  observation.  P}. 

For  each  event  we  now  have  a  scale  size,  Qt.  and  for  each  station  there  is  a 
mean  relative  amplitude.  fi},  and  its  standard  deviation.  S;.  A  measure  of  the 
uncertainty  of  the  mean.  /*.  is  the  standard  error  of  the  mean 

SEM  =  S/VZ  (14) 

where  n  is  the  number  or  measurements  at  the  station,  and  5  the  standard 

deviation  of  the  distribution  of  relative  amplitudes  at  the  station.  The  957,  con¬ 
fidence  limits  on  the  mean  is  roughly  twice  the  standard  error  of  the  mean 

957,  confidence  limits  =  1.96  SEM  (15) 

The  variance  of  the  distribution  of  relative  amplitudes  at  a  station  is  the  square 
of  the  standard  deviation 

Var  =  S*  (16) 

In  compiling  the  results  from  the  amplitude  reduction,  the  station  mean, 
standard  error  of  the  mean,  and  number  are  listed.  From  these  quantities,  the 
others  may  be  easily  calculated  For  the  explosion  data  the  mean  relative  sta¬ 
tion  amplitudes  are  determined  for  each  test  site.  To  obtain  the  mean  relative 


station  amplitudes  for  a  northern  azimuth,  the  results  of  the  five  test  sites  are 
averaged.  This  process  is  accomplished  by  the  same  reduction  as  discussed 
above;  where  now,  instead  of  averaging  amplitudes  from  events,  we  are  averag¬ 
ing  relative  amplitudes  characteristic  of  each  test  site.  We  could  have  obtained 
relative  station  amplitudes  for  a  northern  azimuth  by  averaging  the  amplitude 
measurements  without  regard  to  test  site.  Averaging  by  test  site  first  and  then 
averaging  the  test  sites  was  chosen  for  two  reasons,  (l)  Some  test  sites  had 
more  events  than  others.  If  the  amplitudes  were  averaged  without  regard  to 
test  site,  some  test  sites  would  be  weighted  more  heavily  in  the  average.  (2) 
For  the  east  and  west  Semipalatinsk  and  western  Kazakh  sites,  stations  in  the 
southern  United  States  lie  within  the  region  of  diffraction  at  the  core  mantle 
boundary.  One  of  the  spinoffs  of  this  amplitude  study  is  a  more  precise  deter¬ 
mination  of  the  core-mantle  shadow  zone  boundary  at  95.5°A  (Butler  and  Ruff. 
19B0;  Ruff  and  Helmberger,  1981).  For  these  reasons  the  Russian  test  sites 
were  averaged  separately. 


The  scale  sizes  a,  of  the  events  are  not  actually  used  in  this  study.  Though 
the  a*  are  potentially  good  estimators  of  the  relative  sizes  of  events,  they  meas¬ 
ure  only  a  small  solid  angle  about  the  sources  and  are  therefore  of  restricted 
utility. 

The  mean  relative  station  amplitudes  for  a  given  lest  site  or  azimuth  are 
determined  relative  to  the  group  and  not  to  a  fixed  baseline.  We  may  then 
choose  a  baseline.  The  baseline  chosen  is  that  the  geometric  mean  value  of  the 
mean  relative  amplitudes  over  all  stations  shall  be  equal  to  one.  The  geometric 
mean  7  of  the  mean  relative  station  amplitudes  is 


jn  M,]" 

7=b 


O7) 


where  is  the  product  of  the  means  and  j  goes  from  1  to  n.  the  number  of 


stations.  The  mean  relative  amplitude  is  then  defined 


m,  =  fiij/y  (IB) 

With  the  mean  relative  station  amplitudes  in  this  format,  we  may  then  speak  of 
station  amplitude  as  being  a  factor  of  two  greater  than  the  mean  (twice  the 
amplitude)  or  a  factor  of  two  less  than  the  mean  (one-half  the  amplitude), 
where  the  mean  refers  to  the  geometric  mean  amplitude  over  all  stations. 

Amplitude  Results  for  Three  Azimuths 

For  each  seismic  station  there  are  many  relative  amplitude  determinations, 
each  determination  from  a  particular  source  event.  The  individual  vagaries  of 
the  sources  will  contribute  scatter  to  relative  amplitude  determinations  even 
though  we  have  minimized  visible  waveform  variations.  Unknown  scattering 
bodies  within  the  earth  and  lateral  heterogeneity  will  cause  the  relative  ampli¬ 
tudes  determinations  to  scatter.  Whatever  the  cause  of  the  scatter,  we  have  a 
distribution  of  relative  amplitude  determinations  at  each  seismic  station  for  the 
many  events  with  respect  to  three  azimuths.  This  distribution  of  relative  ampli¬ 
tude  determinations  at  each  station  may  be  characterized  by  the  mean  and  the 

standard  deviation  of  distribution.  The  mean  of  the  distribution  is  an  unbiased 
1 

estimator  of  the  relative  amplitude  of  the  station  with  respect  to  the  network. 
The  standard  deviation  measures  the  scatter  of  the  relative  amplitude  determi¬ 
nations  about  the  mean.  The  variance  is  the  square  of  the  standard  deviation: 
the  sum  of  the  variances  over  all  stations  was  minimized  in  the  amplitude 
reduction  procedure.  Thus  the  scatter  in  the  relative  amplitude  determinations 
is  real,  and  minimum  in  a  least  squares  sense  with  respect  to  all  stations. 

The  relative  amplitude  patterns  for  three  azimuths  for  twenty-five  WWSSN 
and  SDCS  stations  in  the  United  Slates  are  graphed  in  Figures  6.  ?  and  8. 
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BKS  LON  082NV  DUG  B07  SOL  JCT  MON  SNA  ATL  SCP  OGD  HNWF 


IKS  LON  0B2NV  DUG  BOZ  SOL  JCT  RKON  SNA  ATL  SCP  OGD  HKW* 

Figure  6.  Relative  amplitude  of  short  period  P  waves  from  Russian  nuclear 

test  sites  are  graphed  above  for  WWSSN  stations  in  the  United  States. 
The  seismic  stations  are  arranged  by  location  in  a  west  to  east 
direction.  The  corresponding  mean  and  standard  error  of  the  mean  are 
graphed  beneath.  The  vertical  amplitude  scale  is  logarithmic  with 
unity  equal  to  the  geometric  Man  over  all  stations. 
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Figure  7.  Relative  amplitude*  of  short  period  P  waves  from  earthquake*  to  a 

northwest  azimuth  from  the  United  States  are  graphed  above  for  VL’SSN 
stations.  The  seismic  stations  are  arranged  by  location  in  a  vest  to 
east  direction.  The  corresponding  mean  and  standard  error  of  the 
swan  are  graphed  beneath.  The  vertical  amplitude  scale  is  logarithr.ic 
with  unity  equal  to  the  geometric  mean  over  all  stations. 
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EARTHQUAKES  •  SOUTH  AMERICA 


COR  CSC  MSO  TUC  ALO  LUB  DAL  OXF  AAM  BLA  GEO  WES 
aneKS  LON  0B2NV  DOG  BOZ  601  JCT  RKON  SHA  ATI  SCP  OGD  HNWE 

•On — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — n 


COR  CSC  MSO  TUC  ALO  LUB  DAL  OXF  AAM  BLA  GEO  WES 
BKS  LON  082NV  DUG  BOZ  60L  JCT  RKON  SHA  ATL  SCP  OGD  NNWE 

Figure  B.  Relative  amplitudes  of  short  period  P  waves  from  earthquakes  in  South 
America  are  graphed  above  for  WWSSN  stations  in  the  United  States. 

The  seismic  stations  are  arranged  by  location  in  a  west  to  east 
direction.  The  corresponding  mean  and  standard  error  of  the  mean 
are  graphed  beneath.  The  vertical  amplitude  scale  is  logarithmic 
with  unity  equal  to  the  geometric  mean  over  all  stations. 
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Amplitude  data  from  a  northern  azimuth  come  from  mean  relative  station 
amplitudes  determined  for  five  Russian  nuclear  test  sites  and  are  graphed  in 
the  upper  part  of  Figure  6.  The  mean  relative  amplitudes  were  obtained  using 
the  amplitude  measurements  from  Butler  and  Ruff  (19S0).  To  include  the  sta¬ 
tions  RKOX,  HNME,  and  0B2NV,  three  additional  explosions  from  the  Semipala- 
tinsk  east  test  site  were  integrated  into  the  Butler  and  Ruff  (1930)  amplitude 
data  set.  Even  though  this  presents  RKON,  HNME.  and  0B2NV  with  one  to  three 
measurements,  the  overall  stability  of  relative  amplitudes  from  individual 
nuclear  test  sites  suggests  that  the  individual  amplitude  variations  will  be  close 
to  the  mean  relative  amplitude  patterns  determined  from  many  events.  The 
mean  relative  amplitudes  for  the  twenty-five  stations  are  listed  left  to  right  in  a 
west  to  east  fashion  as  they  are  located  in  the  United  States  for  a  northern 
azimuth.  The  vertical  scale  is  logarithmic  and  the  axis  of  unity  is  the  geometric 
mean  of  the  distribution  of  data.  The  dashed  parallel  lines  indicate  the  factor 
of  two  variation  about  the  mean,  or  relative  amplitudes  of  2  and  0.5.  To  make 
clear  the  number  of  test  sites  when  the  data  overlap,  the  X’s  and  A  s  indicate 
two  and  three,  respectively,  test  sites  share  a  relative  amplitude  (to  within  the 
size  or  the  symbol).  Some  stations  have  five  determinations,  some  less  due  to 
distance  range  or  variations  of  reporting  (Butler  and  Ruff.  1980).  Figure  6  also 
graphs  the  mean  and  the  standard  error  of  the  mean  for  the  data  for  the 
northern  azimuth.  If  the  standard  error  of  the  mean  was  smaller  than  the  sym¬ 
bol  size,  it  was  not  graphed.  Individual  measurements  do  not  have  a  standard 
error  defined  and  are  shown  as  open  circles. 

This  graphical  representation  of  the  data  was  chosen  to  best  illustrate  the 
interrelationship  of  the  data,  its  mean  and  its  standard  error  among  the  many 
stations.  The  west  to  east  arrangement  of  the  stations  was  chosen  to  provide 
geographical  context  to  the  data  and  to  make  visible  possible  west  to  east 
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amplitude  variations.  The  three  letter  seismograph  code  names  refer  to  the 
locations  on  the  map  of  Figure  1;  Table  1  lists  the  station  names. 

Two  stations  are  a  factor  of  two  less  than  the  mean  for  all  test  sites:  ALQ. 
Albuquerque.  New  Mexico,  and  GOL.  Golden.  Colorado.  Three  stations  have  mean 
amplitudes  consistently  a  factor  of  two  greater  than  the  mean:  DAL.  Dallas. 
Texas;  RKOX.  Red  Lake.  Ontario;  and  KNME.  Houlton,  Maine.  The  stations  0B2XV, 
Oak  Springs  Butte,  Nevada,  BOZ.  Bozeman.  Montana,  LUB,  Lubbock.  Texas.  SKA, 
Spring  Hill.  Alabama  and  AAM.  Ann  Arbor,  Michigan,  measure  above  the  mean. 
Stations  consistently  measuring  below  the  mean  are  DUG.  Dugway.  Utah.  TUC, 
Tucson,  Arizona.  JCT.  Junction  City,  Texas,  ATL.  Atlanta.  Georgia,  OGD. 
Ogdcnsburg.  New  Jersey,  and  HITES,  Weston.  Massachusetts.  The  rest  of  the  sta¬ 
tions  have  relative  amplitudes  falling  on  both  sides  of  the  mean:  BKS.  Berkeley, 
California.  COR,  Corvallis.  Oregon.  LON.  Longmire,  Washington.  CSC.  Goldstone, 
California.  MSO,  Missoula.  Montana.  OXF.  Oxford,  Mississippi.  BLA.  Blacksburg. 
Virginia.  SCP,  State  College.  Pennsylvania,  and  GEO,  Georgetown.  District  of 
Columbia.  The  stations  which  show  the  greatest  scatter  and  consequently  have 
the  largest  standard  errors  of  the  mean  are  LON,  BLA,  and  SCP.  The  maximum 
amplitude  ratio  in  the  amplitude  pattern  for  a  northern  azimuth  is  between 
RKOX  und  GOL.  where  the  mean  amplitude  of  RKON  is  15  times  greater  than  GOL 
The  amplitudes  at  DAL  and  HNME  are  6  to  7  times  greater  than  the  mean  ampli¬ 
tude  at  GOL. 

Amplitude  data  for  a  northwest  azimuth  from  the  United  States  are 
graphed  in  Figure  7.  The  formats  or  all  of  the  amplitude  figures  in  this  paper 
arc  all  the  same  to  allow  relative  ease  in  comparing  one  set  of  data  to  another 
The  data  in  Figure  7  are  from  21  earthquakes  in  Kurile-Japan-Bonin  arc.  Kam¬ 
chatka.  Eastern  Siberia.  Aleutian  arc.  and  Alaska.  The  data  in  Figure  7  have 
similarities  and  differences  with  the  amplitude  pattern  for  a  northern  azimuth 
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in  Figure  6.  LON,  TUC.  and  ALQ  consistently  report  amplitudes  less  than  the 
mean,  while  COL  observes  amplitudes  consistently  a  factor  of  two  less  than  the 
mean.  The  amplitude  observations  at  RKON  and  OXF  are  consistently  greater 
than  the  mean.  High  amplitudes  are  noted  at  AAM.  but  also  with  relatively  large 
scatter.  Other  stations  with  large  scatter  are  LON  and  TUC.  In  contrast  to  LON 
and  AAM,  most  of  the  amplitudes  at  TUC  are  centered  near  the  mean  relative 
amplitude  and  consequently  the  standard  error  of  the  mean  at  TUC  is  small. 

Several  stations  show  significant  azimuth  amplitude  variation  between 
north  and  northwest.  DUG,  Dugway,  Utah  and  WES,  Weston.  Massachusetts  have 
mean  relative  amplitudes  that  are  a  factor  of  two  less  than  the  network  mean 
for  Russian  explosions  to  a  northern  azimuth,  and  shift  to  amplitudes  greater 
than  the  network  mean  for  the  earthquakes  to  a  northwest  azimuth.  For  LON, 
Longmire,  Washington,  the  situation  is  the  reverse;  where  greater  than  average 
amplitudes  are  observed  for  nuclear  explosions  to  the  northern  azimuth,  an 
amplitude  a  factor  of  two  less  than  the  network  mean  is  observed  for  earth¬ 
quakes  to  a  northwest  azimuth.  The  stations  ATL  and  BLA  show  azimuthal  varia¬ 
tions.  but  to  a  lesser  extent  than  do  LON,  DUG,  and  WES.  The  range  of  ampli¬ 
tudes  is  slightly  less  for  the  northwest  azimuth  than  for  the  northern  azimuth. 
The  amplitude  ratio  between  RKON  and  GOL  is  5.2  and  is  exceeded  by  the  ampli¬ 
tude  ratio  of  8.1  between  OXF  and  COL.  While  ALQ  is  nearly  as  low  as  COL  for  the 
explosion  data,  the  relative  amplitude  of  ALQ  increases  for  earthquakes  to  a 
northwest  azimuth. 

The  earthquakes  to  a  northwest  azimuth  cover  a  broad  region  These 
regions  were  subdivided  in  Table  3  and  in  Appendix  Tables  A7  through  A10  to 
tee  if  systematic  source  region  effects  could  be  resolved.  Four  figures  are 
shown  in  the  Appendix  -  Figures  A1  through  A4  -  which  jointly  graph  the  mean 
relative  amplitude  pattern  from  all  data  to  the  northwest  versus  the  amplitude 


patlcrns  from  the  four  regional  groupings  in  Table  3.  No  systematic  source 
region  effects  were  apparent. 

The  amplitude  pattern  for  a  south  to  southeast  azimuth  from  the  United 
States  is  graphed  in  Figure  B.  Data  from  fifteen  South  American  earthquakes 
are  plotted  in  the  upper  part  of  Figure  B.  The  seismic  stations  COR.  0B2NV.  LUB. 
and  RKON  consistently  have  amplitudes  greater  than  or  equal  to  the  network 
mean.  The  two  measurements  at  HNME  are  greater  than  the  network  mean. 
Lower  than  average  amplitudes  are  consistently  observed  at  BOZ,  ALQ,  COL.  SCP. 
and  WES.  The  single  measurement  at  BOZ  is  less  than  half  the  network  mean, 
and  was  the  only  measurement  available  from  the  earthquake  data  set  in  Table 
4.  Lay  and  Kelmberger  (19B1)  obtained  three  amplitude  measurements  at  BOZ 
from  four  deep  focus  earthquakes  in  Peru;  low  relative  amplitudes  were  deter¬ 
mined  at  BOZ.  However.  TUC,  ALQ.  and  COL  recorded  smaller  amplitudes  than 
BOZ  in  the  Lay  and  Kelmberger  (1981)  data  set.  The  amplitude  determinations 
for  LON  in  Figure  7  are  again  generally  low,  with  large  scatter.  Other  stations 
exhibiting  relative  large  scatter  of  amplitude  observation  from  the  data  set  of 
South  American  earthquakes  are  DAL,  OXF.  ATL,  BLA,  and  OCD.  It  appears  that 
more  stations  of  the  WWSSN  in  the  United  Stales  exhibit  large  scatter  for  a 
south-southeast  azimuth  than  for  a  northwest  azimuth. 

A  number  of  azimuthal  amplitude  variations  are  seen  at  single  stations. 
COR  observes  its  highest  amplitudes  to  a  southeast  azimuth.  LON  shows  low 
amplitudes  for  the  South  American  earthquakes,  in  agreement  with  a  northwest 
azimuth  and  in  contradistinction  to  high  amplitudes  observed  from  Russian 
nuclear  events.  DUC  observed  amplitudes  near  the  netw-ork  mean  for  a 
southeast  azimuth.  BOZ  shows  low  amplitudes  to  the  south,  high  amplitudes  to 
north  and  northwest.  COL  shows  low  amplitudes  to  the  southeast,  but  not  to 
the  extent  a  factor  of  two  less  than  the  network  mean  -  as  was  the  case  for  GCL 
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to  the  north  and  northwest.  LUB  shows  high  amplitudes  for  a  southern 
azimuth.  WES  has  low  amplitudes  to  the  south  and  north,  while  higher  ampli¬ 
tudes  are  observed  to  the  northwest. 

The  lower  amplitudes  seen  at  BOZ.  Bozeman.  Montana  from  South  American 
earthquakes  are  more  than  likely  associated  with  the  Yellowstone  hot  spot. 
High  temperatures  beneath  Yellowstone  will  lead  to  greater  attenuation  of 
seismic  waves  passing  through  the  region.  The  large  regional  geoid  high  in  the 
western  United  States  is  centered  on  Yellowstone.  However,  though  a  straight¬ 
forward  cause  may  be  noted  for  the  low  amplitudes  at  BOZ  for  a  southern 
azimuth,  equally  low  amplitudes  are  observed  at  COL,  ALQ.  and  TUC.  In  inter¬ 
preting  the  amplitude  patterns  for  the  network,  consideration  must  be  given  to 
individual  station  effects,  but  effects  noted  in  common  at  many  or  several  sta¬ 
tions  must  be  considered  together. 

The  relative  amplitude  pattern  for  a  south-southeast  azimuth  in  the  lower 
part  of  Figure  3  shows  a  smaller  range  of  amplitude  variation  between  the  high 
and  low  stations  than  was  seen  for  the  northern  azimuth.  The  relative  ampli¬ 
tude  ratio  between  DAL  and  BOZ  for  the  South  American  data  is  6.7,  while  the 
ratio  of  RKON  to  COL  is  3.2.  This  range  of  amplitudes  is  comparable  to  the 
range  for  northwest  azimuth. 

The  amplitude  data  for  a  northwest  azimuth  were  subdivided  into  four 
regions  to  see  if  near  source  regional  effects  might  be  delected  in  the  data  No 
systematic  relationships  were  apparent.  While  there  were  no  obvious  regional 
subgroups  into  which  to  divide  the  South  American  amplitude  data,  t^e  indivi¬ 
dual  events  were  plotted  in  relation  to  the  mean  relative  amplitude  pattern 
from  all  the  data  to  see  if  any  of  the  events  showed  systematic  trends.  None  of 
the  individual  South  American  events  evinced  an  apparent  biasing  trend  with 
respect  to  the  mean  relative  amplitude  pattern  for  the  south-southeastern 
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Station  Kean  Relative  Amplitude  Summary 

A  summary  table  of  the  station  mean  relative  amplitudes  is  presented  in 
Table  5  for  the  three  azimuths  in  this  study.  The  standard  error  of  the  mean 
and  the  number  of  observations  is  given.  For  the  northwest  and  south- 
southeast  azimuths  the  number  of  observations  refers  to  the  number  of  earth¬ 
quakes.  For  the  northern  azimuth  the  number  of  observations  refers  to  the 
number  of  Russian  nuclear  test  sites.  The  earthquake  data  have  been 
corrected  for  geometric  spreading  to  a  common  distance  of  BOcf  A.  The  explo¬ 
sion  data  have  not  been  corrected  for  geometric  spreading,  and  thus  Lhe  north¬ 
erly  stations  are  about  10%  too  large  and  the  southerly  stations  10%  too  small. 
Individual  summary  lists  for  Russian  test  sites  are  found  in  Tables  A 12  through 
A17  in  the  Appendix  and  were  discussed  in  Butler  and  Ruff  (1900).  Individual 
summary  lists  for  the  four  regional  subgroups  of  earthquake  data  to  a 
northwest  azimuth  are  found  in  the  Appendix  as  Tables  A18  through  A21. 

Mean  Amplitude  Pattern  for  the  United  States  and  Southern  Canada 

To  obtain  a  mean  amplitude  pattern  for  the  United  States  and  soi’U-crn 
Canada,  the  mean  relative  amplitudes  for  the  three  azimuths  -  N,  N’«v,  -ind  SEE  - 
were  weighted  equally  and  averaged  together-.  The  averaging  process  va? 
accomplished  in  the  same  way  that  the  dat  a  were  reduced  -  by  nvnimizai  'or  of 
total  variance.  The  mean  relative  amplitudes  for  the  United  States  are  graphed 
in  Figure  9.  The  error  bars  in  Figure  9  contain  the  standard  error  of  the  mean 
over  all  azimuths.  The  pattern  from  left  to  right  is  one  of  mean  amplitudes 
decreasing  from  the  west  coast  eastwards,  reaching  a  low  at  GOL  on  the  Rocxy 
Mountain  Front  Range,  jumping  to  greater  than  the  mean  in  Texas,  reaching  its 
height  at  RKON  on  the  Canadian  shield,  and  decreasing  toward  the  mean  toward 
the  east  coast  stations.  The  striking  feature  of  Figure  9  is  the  amplitude 
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MEAN  P  WAVE  AMPLITUDE  VARIATION  FOR  U.  S.  A. 


Figure  9.  The  relative  amplitude  variation  for  the  N,  NW,  and  SSE  azimuths 
are  averaged  to  obtain  a  mean  P  wave  amplitude  variation  for  the 
United  States  and  southern  Canada.  The  error  bars  represent  the 
total  range  of  the  standard  errors  of  the  means  for  the  three 
azimuths . 


Table  5 

Mean  P  Wave  Amplitudes  and  Standard  Error  of  the  Mean 


Russian  Explosions  All  Earthquakes  South  American 
to  North  to  Northwest  Earthquakes 


Station 

Mean 

S.E.M 

N 

Mean 

S.E.M 

N 

Mean 

S.E.M. 

N 

(west) 

BKS 

1.10 

0.14 

4 

1.26 

0  16 

13 

1.09 

0  06 

10 

COR 

0.B9 

0.10 

3 

1.30 

0.1? 

6 

1.63 

0.16 

7 

LON 

3.2? 

0.26 

5 

0.3B 

0  10 

7 

0.52 

0.13 

6 

CSC 

0.75 

Oil 

4 

0.76 

0.06 

13 

0.99 

0.09 

12 

032 MV 

1.22 

10.10]* 

1 

0.91 

0.06 

7 

1.39 

0.16 

4 

MSO 

0.67 

0  16 

3 

0.75 

0.09 

12 

0.6? 

0.10 

6 

DUG 

0.54 

0.09 

2 

1 .16 

0.10 

17 

1.0? 

0.12 

12 

TUC 

0.66 

0  02 

2 

0.51 

0  04 

19 

0.69 

0.07 

14 

B02 

1.44 

[0  10]* 

1 

1.24 

0  13 

4 

0.38 

[014]* 

1 

ALQ 

0.37 

0.04 

4 

0.5? 

0.05 

12 

0.68 

0  06 

7 

COL 

0  32 

0  06 

3 

0  34 

0  03 

12 

0.58 

0  03 

6 

Lua 

1.15 

0.12 

2 

1.02 

0  09 

6 

2  15 

0  31 

4 

JCT 

0.65 

0.08 

2 

1.22 

0.24 

4 

1.26 

0  76 

2 

DAL 

2  23 

[0.10]* 

i 

3 .42 

0  28 

4 

2.53 

1.20 

3 

RXOS 

4.69 

[0  10]* 

i 

1.7? 

0.20 

6 

1.65 

0.30 

4 

OX? 

1.19 

0  08 

4 

2.09 

0  16 

4 

2.10 

0  59 

4 

Slit 

1.14 

0.15 

2 

1.14 

10  12]* 

1 

- 

- 

- 

AAM 

1.48 

0  10 

5 

1.64 

0  49 

5 

0.61 

0  22 

2 

ATL 

0.60 

0.07 

5 

1.14 

0.11 

6 

0  94 

0  14 

11 

BLA 

0.6) 

0.20 

5 

1.44 

0.16 

6 

1.11 

0.23 

6 

SCP 

1.26 

0.25 

5 

0.95 

0.11 

7 

0  75 

0.04 

5 

CEO 

1.40 

0,20 

3 

0  92 

0  15 

3 

0.94 

0  17 

6 

OCD 

0  65 

0.09 

5 

0.65 

0.22 

4 

0  58 

0  12 

6 

WES 

0.52 

0.1" 

4 

1.2? 

0.13 

5 

0  09 

0  02 

3 

1INME 

(east) 

2.14 

1010]* 

1 

1.00 

[0  12]* 

1 

110 

0.00 

2 
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variation  between  the  Rocky  Mountain  Front  Range  station  GOL  and  the  Supe¬ 
rior  Province  Canadian  shield  station  RKON.  The  change  between  low  ampli¬ 
tudes  at  GOL  and  at  ALQ  in  the  Rio  Grande  rift  and  the  moderate  to  high  ampli¬ 
tudes  in  Texas  is  sudden.  This  is  the  boundary  between  the  late  Cretaceous  and 
Cenozic  orogeny  of  the  Rocky  Mountains  and  the  Great  Plains  -  a  craton  stable 
since  Phanerozoic  times.  The  amplitude  jump  also  marks  a  change  of  seismom¬ 
eter  siting  from  hard  rock  sites  in  the  west  to  sedimentary  siting  in  the  central 
United  States.  The  highest  amplitude  in  Figure  9  lies  on  a  section  of  the  Cana¬ 
dian  shield  with  a  2.5  billion  year  age  (Stockwell,  1964).  This  amplitude  pattern 
is  similar  to  previous  LRSM  amplitude  studies,  where  low  amplitudes  are  found 
in  the  west  and  high  amplitudes  in  the  east  (Cleary,  1967;  Evernden  and  Clark, 
1970;  Booth  et  al.,  1974).  As  in  the  previous  studies  there  is  variability.  In  Fig¬ 
ure  9  the  only  consistently  low  amplitude  stations  are  MSO,  TUC,  ALQ.  and  GOL. 
The  consistently  high  amplitude  stations  are  DAL.  RKON,  QXF,  and  SHA.  Three  of 
the  four  high  amplitude  stations  -  DAL,  OXF,  and  SHA  -  lie  along  the  Gulf  Coast 
and  are  situated  on  low  velocity  sediments  which  effectively  amplify  the  ground 
signal.  Four  high  stations  not  situated  on  sediments  are  COR.  0B2NV,  RKON,  and 
HNME.  The  high  amplitudes  at  0B2NV  lie  within  the  Great  Basin  province  of 
Nevada  in  an  area  of  low  observed  heat  flow  (Sass  et  al.,  1971,1976). 

The  mean  relative  amplitudes  from  Figures  6,  7,  and  B  are  presented  in  Fig¬ 
ure  10  in  relation  to  simplified  physiographic  provinces  of  the  United  States.  In 
general  three  symbols  are  plotted  at  each  station  representing  the  mean  rela¬ 
tive  amplitudes  for  the  three  azimuths.  The  symbols  are  located  N,  NW,  and  SSE 
of  the  station  location  to  indicate  direction  of  azimuth.  Upward  pointing  trian¬ 
gles  represent  mean  relative  amplitudes  which  are  factors  greater  than  the 
network  geometric  mean  amplitude.  Downward  pointing  triangles  represent 
mean  relative  amplitudes  which  are  factors  less  than  the  network  geometric 
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SHORT  PERIOD  P  WAVE  RELATIVE  AMPLITUDES 


Figure  10.  Short  period  P-wave  amplitude  anomalies  are  plotted  relative  to  station 

location  and  approximate  physiographic  province  within  the  United  States. 
Upward  pointing  triangles  are  values  greater  than  the  mean,  downward 
pointing  are  values  less  than  the  mean.  A  circle  indicates  a  mean  ampli¬ 
tude.  The  amplitude  factor  relative  to  the  mean  is  noted  at  the  lower 
right.  The  position  of  the  triangle  relative  to  each  station  indicates 
azimuth:  north,  northwest ,  southeast.  The  data  to  the  north  are  from 
explosions  in  Russia.  The  data  to  the  northwest  and  southeast  azimuths 
are  from  earthquakes  in  Asia  and  South  America  with  simple  sources. 


I 


mean.  The  circles  represent  mean  relative  amplitudes  which  are  at  the  network 
geometric  mean.  The  size  of  the  symbol  is  scaled  to  the  deviation  from  the  net¬ 
work  geometric  mean.  For  example,  to  a  NW  azimuth.  HNME  in  Maine  records  an 
amplitude  at  the  network  geometric  mean.  TUC  in  Arizona  records  an  amplitude 
of  0.5  or  a  factor  of  two  less  than  the  geometric  mean,  while  RKON  in  Ontario 
records  an  amplitude  a  factor  of  two  greater  than  the  network  mean.  Supple¬ 
mental  data  for  RCD,  Rapid  City,  South  Dakota;  MDS,  Madison,  Wisconsin;  FLO, 
Florissant,  Missouri;  and  BEC,  Bermuda  were  taken  from  Butler  and  Ruff  (i960). 
Many  of  the  features  of  this  figure  have  been  discussed  in  previous  sections: 
lowest  amplitudes  in  the  Rockies  at  GOL  and  south  of  BOZ,  and  in  the 
southeastern  Basin  and  Range  at  TUC  and  ALQ;  lower  amplitudes  to  the 
northwest  of  GSC  and  0B2NV,  north  DUG,  and  between  COR  and  LON;  scattered 
high  amplitudes  in  the  west;  high  amplitudes  in  the  center  of  the  continent;  and 
moderate  high  and  low  amplitude  variation  in  the  Appalachias  and  at  the  east 
coast.  Note,  however,  that  the  amplitudes  reported  have  not  been  corrected  for 
any  site  amplification  effects  for  the  seismometers  situated  on  low  velocity  sed¬ 
iments. 

Amplification  Effects  of  Shallow  Sediments 

The  amplification  of  seismic  waves  recorded  at  seismometers  situated  cn 
soft  sediments  relative  to  seismometers  situated  on  hard  bedrock  has  been 
known  for  some  time.  Literature  references  include  Gutenberg  (1956,  1957), 
Carpenter  et  al.  (1967),  Borcherdt  (1970),  Booth  et  al.  (1974),  and  Butler  and 
Ruff  (1980).  The  effect  may  be  understood  fairly  easily  by  comparing  two  elas¬ 
tic  springs  with  different  spring  constants  -  one  relatively  stiff  and  one  rela¬ 
tively  flexible.  Let  us  suppose  that  the  springs  are  free  at  one  end,  as  an  analog 
to  the  free  surface  of  the  earth  where  the  seismometer  is  placed.  Now  corsider 


elastic  waves  with  the  same  energy  propagating  through  both  springs.  As  the 
elastic  waves  reach  the  free  ends  of  the  springs,  the  relatively  flexible  spring 
has  more  "give"  and  stretches  further  than  the  stiff  spring.  The  greater  stretch 
in  the  flexible  spring  relative  to  the  stiff  spring  is  equivalent  to  the  larger 
amplitudes  seismic  waves  produced  in  sediments  versus  hard  rock.  For  the  P 
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waves  measured  in  this  amplitude  study  the  apparent  period  is  about  1  second 
and  the  phase  velocity  (at  a  representative  distance  of  60"A)  is  about  17 
km/sec,  such  that  the  wavelength  of  the  P  waves  is  on  the  order  of  17  km.  The 
question  then  is:  given  the  relatively  long  wavelengths,  how  much  effect  will 
shallow  sediments  have  upon  the  P  wave  amplitudes?  Butler  and  Ruff  (1950) 
performed  numerical  experiments  using  Thomson-Kaskell  propagator  matrices 
and  a  synthetic  explosion  waveform  to  derive  the  receiver  effect  of  a  variety  of 
different  sediment  structures  overlying  bedrock  relative  to  a  simple  bedrock 
site.  For  shallow  sediments  having  an  average  velocity  of  2  km/sec,  a  200  meter 
thick  surface  layer  amplified  the  P  waves  by  a  factor  of  1.1  while  thicknesses 
300  meters  and  more  amplified  the  P  waves  by  a  factor  of  1.4  to  l.B.  Shallow 
sediments  with  higher  velocities  produced  less  amplification.  The  relative 
amplification  of  sediment  sites  versus  bedrock  seismometer  sites  is  a  function 
of  the  sediment  thickness  and  sediment  velocity.  Table  11  lists  the  results  of 
the  numerical  experiments  of  Butler  and  Ruff  (19B0)  for  sediment  layers  of  con¬ 
stant  velocity  and  velocity  gradients  for  various  velocities  and  thicknesses.  The 
average  velocity  in  the  top  1/2  kilometer  essentially  determines  the  relative 
amplification  of  the  P  waves  in  this  study.  As  most  of  the  central  and 
southeastern  portions  of  the  United  States  are  covered  by  sediments,  the  effect 
of  sediment  amplification  is  important  in  understanding  the  variation  of  P  wave 
amplitudes. 

To  provide  a  standard  bias  of  comparison  of  P  wave  relative  amplitudes, 
those  seismic  stations  situated  on  sediments  must  be  corrected  for  the  sedi¬ 
ment  amplification  effect  relative  to  the  stations  situated  on  hard  bedrock 
Seismic  velocity  well  logs  provide  accurate  velocity  control  of  surface  sediment 
layers,  but  this  type  of  information  was  obtainable  only  at  SHA.  Spring  Hill.  Ala¬ 
bama.  Short  (<  1 0  km)  seismic  survey  lines  are  sensitive  to  uppermost  crustal 
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Table  1 1 

Sediment  Amplification  Effects  for  P  Waves* 


Sediment 

amplification** 

Sediment  description*** 

P  wave  velocitv.  thickness 

l.B 

unconsolidated  sediments. 

2.0  km/sec  at  surface  and 

thicknesses  of  2+  km  with 

increasing  velocity  w/depth 

1.6 

2.0  km/sec  at  surface  and 

thickness  of  400-500  meters 

2.5  km/sec  at  surface  and 

thicknesses  of  1-2  km  with 

increasing  velocity  w/depth 

1.5 

3.0  km/sec.  2  km  thick 

1.4 

2  km/sec,  300  meters  thick 

3  km/sec.  3  km  thick  with 

increasing  velocity  w/depth 

1.3 

3.5  km/sec,  2+  km  thick  with 
increasing  velocity  w/depth 

1.2 

3.0  km/sec.  500  meters 

1.1 

2  km/sec,  s2DD  meters  thick 

•  at  a  period  T  ~  1  sec. 

••  relative  to  a  bedrock  site. 

•••  sediment  layers  overlay  Vp  =  G  O  km/sec  bedrocl 


layers;  however,  only  at  BKS.  Berkeley,  California  was  such  a  survey  carried  out 
and  made  available  (T.  McEvilly,  personal  communication,  19B2).  Long  seismic 
refraction  survey  lines  provide  only  general  control  as  sediments  are  grouped 
in  a  single  layer  characterized  by  the  highest  velocity  arrival.  To  estimate  shal¬ 
low  sediment  velocities  and  hence  from  Table  11  the  relative  sediment  amplifi¬ 
cation  where  seismic  refraction  lines  or  well  log  control  is  not  available  -  this  is 
most  often  the  case  -  recourse  must  be  taken  from  geologic  site  information 
and  systematics  between  velocity,  and  lithologic  type  and  age  of  the  sediments. 
Limestone  velocities  are  sensitive  to  the  extent  of  crystallization  and  range 
between  1.7  and  6.1  km/sec  (Press.  1966).  Faust  (1951)  developed  empirical 
systematics  of  shale  and  sandstone  seismic  velocities  as  a  function  of  depth  and 
geologic  time  from  velocity  data  from  over  500  wells  in  the  United  States  and 
Canada.  Marrying  the  shale  and  sandstone  velocity-depth-age  systematics  of 
Faust  (1951)  to  the  qualitative  sediment  velocity- amplification  relationship  in 
Table  11.  an  approximate  geologic  age  versus  relative  amplification  relationship 
was  derived  for  sandstone  and  shale  sediments  with  thicknesses  of  300  meters 
or  greater.  Table  12  lists  this  relationship  of  shale  and  sandstone  geologic  age 
and  amplification  of  P  waves  relative  to  hard  bedrock  sites.  Naturally,  the  rela¬ 
tionships  in  Table  12  must  be  regarded  as  only  approximate,  and  were 
developed  only  from  the  necessity  to  apply  corrections  for  sediment  (shale  and 
sandstone)  amplification  in  a  systematic  fashion.  Velocity  information  on  other 
sedimentary  materials  -  chalk,,  gypsum,  mud.  etc.  -  were  obtained  from  Press 
(19CC)  and  appropriate  amplification  corrections  were  deduced  from  Table  11. 
In  interpreting  the  P  wave  amplitude  variations  in  the  upcoming  section, 
corrections  for  the  amplification  effects  of  surface  sediment  cover  at  the 
seismometer  sites  will  be  adopted  from  Tables  11  and  12. 
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Table  12 

Seismic  Amplification  for  Seismometer 
Site  on  Sandstone  and  Shale  Relative  to 
Hard  Bedrock* 


Geologic  Age 
of 

Sandstone  and  Shale 

Relative  P  Wave 
Amplification** 
(ADDroximatel 

Precambrian 

1.0 

Lower  Paleozoic 
Cambrian 
Ordovician 
Silurian 
Devonian 

1.3 

Upper  Paleozoic 
Mississippian 
Pennsylvanian 
Permian 

1.4 

Mesozoic 

1.6 

Cenozoic 

l.B 

Vp  =  6.0  km/sec 

Stratigraphic  thicknesses  Z  0  meters  or  greater. 


■  Focussing /DcfocussiDg  and  Scatter  Effects  at  Seismometer  Sites 

Fluctuations  of  P  wave  amplitudes  due  to  local  elastic 
focussing/defocussing  and  scattering  effects  are  obviously  likely  contributors 
to  observed  P  wave  amplitude  variations  in  the  United  States.  In  the  amplitude 
stability  section  of  this  paper  the  variability  of  P  wave  relative  amplitudes 
among  seismic  stations  in  a  region  was  examined.  This  section  focusses  on 
more  site  specific  effects. 

Dipping  planar  structures  beneath  a  seismometer  site  tend  to  diminish  the 
amplitude  of  a  P  wave  recorded  on  the  vertical  component  of  a  seismometer  by 
transferring  elastic  energy  to  the  horizontal  components  of  P  wave  motion.  The 
existence  of  transverse  or  tangential  components  of  P  wave  move  is  diagnostic 
of  dipping  structure  beneath  a  particular  site  (Burdick  and  Langston,  197?; 
Langston.  1977).  Model  studies  of  dipping  planar  structure  by  Langston  (197?) 
show  that  a  10°  dip  on  the  Moho  beneath  a  seismic  station  can  reduce  the  P 
wave  amplitude  on  the  vertical  component  of  motion  of  the  long  period 
seismometer  by  10%.  Anomalous  tangential  components  of  P  wave  motion  have 
been  noted  at  COR,  Corvallis.  Oregon  (Burdick  and  Langston.  1977;  Langston. 
1977;  Langston.  1931)  and  at  LON,  Longmire.  Washington  (Langston.  1979).  This 
suggests  that  the  mean  amplitudes  of  P  waves  at  COR  and  LON  -  which  were 
measured  from  the  vertical  components  of  motion  of  the  short  period  V.V.SSN 
seismograhs  -  are  slightly  underestimated  relative  to  other  WWSSN  stations  in 
the  United  States. 

Seismic  stations  situated  on  sedimentary  basins  can  experience  elastic 
focussing/dcfocussing  effects.  Calibrated  nuclear  explosions  at  Yucca  rials  ir. 
Nevada  *  a  long,  trough-like  sedimentary  basin  offset  by  several  large  bloc* 
faults  in  its  basement  -  show  a  systematic  trend  in  their  global  average  magni¬ 
tudes,  aij,  from  one  side  of  the  basin  to  the  other  (Hart  et  a!.,  1979).  The  distri- 


-46- 


bution  of  magnitudes  bas  a  standard  deviation  of  (0.13)  which  is  equivalent  to 
an  amplitude  variation  of  about  a  factor  of  1.3  about  the  mean  for  the  basin. 
By  reciprocity  arguments  (see,  for  example,  the  appendix  of  Chang  and  von 
Seggern,  1980)  mean  relative  amplitudes  of  seismometers  located  at  different 
places  in  Yucca  Flats  will  also  vary.  This  variability  has  been  observed  and 
quantified  using  minimum  entropy  deconvolution  techniques  and  modeled  by- 
Hart  et  al.  (1979)  using  the  glorified  optics  method  of  Kong  and  Kelmberger 
(1979).  Hart  et  al.  (1979)  demonstrated  that  changes  of  seismic  station  loca¬ 
tions.  at  Yucca  Flats  by  several  kilometers  can  change  apparent  magnitudes  by 
(0.1),  equivalent  to  a  factor  of  1.25  in  amplitude.  Note  that  while  the  apparent 
source  strengths  of  the  underground  explosions  exhibited  a  factor  of  1.3  aver¬ 
age  variation  about  the  mean  (i.e.,  greater  than  and  less  than  the  mean),  the 
receiver  studies  of  seismic  stations  at  the  surface  of  the  basin  showed  a  total 
-variation  of  about  a  factor  of  1.25  between  the  largest  and  smallest  amplitude. 
The  receiver  function  studies  et  the  Yucca  Flats  sedimentary  basin  indicate 
that  seismic  stations  situated  on  complex  sedimentary  structures  may  exhibit 
amplitude  variations  of  a  factor  of  1.1  to  1.2  about  the  mean  amplitude  value 
for  the  basin. 

Haddon  and  Kusebye  (1976)  have  modeled  a  smoothed  version  of  the  NOR¬ 
SAR  amplitude  anomalies  with  a  two  dimensional  heterogeneous  layer  at  depths 
around  150  to  200  km  or  the  bottom  of  the  lithosphere.  The  velocity  perturba¬ 
tions  required  to  account  for  the  NORSAR  amplitude  observations  were  of  the 
order  of  a  few  per  cent.  Haddon  and  Husebyc  (197B)  indicate  that  the  thin  lens 
model  does  a  good  job  matching  P  wave  travel  time  anomalies  at  NORSAR.  This 
result  was  said  to  demonstrate  that  travel-time  and  amplitude  anomalies 
observed  across  the  large-apcrturc  array  NORSAR  arc  intimately  correlated  and 
consequently  have  the  same  origin  (Haddon  and  Husebyc.  1978).  In  their 
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azimuthal  study  of  LASA  amplitude  variations  Chang  and  von  Seggern  (1980) 
state  that  the  LASA  amplitude  anomalies  are  linearly  related  to  travel  time 
anomalies,  with  large  amplitudes  being  related  to  late  travel  times.  These 
amplitude  and  travel  time  variations  at  LASA  were  discussed  in  relation  to 
crust-mantle  models  of  Greenfield  and  Sheppard  (1969)  and  Iyer  (1971)  which 
involve  a  varying,  dipping  Moho  beneath  LASA.  As  the  Haddon  and  Kusebye 
(1979)  study  also  showed  a  positive  correlation  between  amplitudes  and  travel 
time  for  P  waves  at  both  LASA  and  NORSAR,  then  large  amplitudes  relate  to  slow 
travel  times  and  small  amplitudes  relate  to  fast  times. 

The  success  of  Haddon  and  Husebye  (1979)  in  mapping  NORSAR  P  wave 
amplitude  and  travel  time  anomalies  into  varying  elastic  structure  beneath 
NORSAR  suggests  that  some  part  of  the  amplitude  variations  in  the  United 
States  found  in  this  present  study  and  in  others  are  due  to  elastic 
focussing/defocussing  effects  of  laterally  varying  uppermost  mantle  structure. 
Using  the  three-dimensional  block  velocity  model  of  the  upper  mantle  beneath 
the  United  states  of  Romanowicz  (1979),  one  might  expect  to  be  able  to  esti¬ 
mate  areas  of  focussing  or  defocussing.  Interestingly,  one  can  easily  show  that, 
statistically.  P  wave  amplitude  and  travel-time  variations  do  not  correlate  well 
for  either  the  WWS5N  station  or  the  LRSM  stations  in  the  United  States.  In  com¬ 
paring  WWSSN  amplitudes  presented  here  with  P  wave  travel  time  anomalies 
from  Romanowicz  (1979)  and  the  LRSM  amplitudes  of  Cleary  (196~)  with  the 
LRSM  travel  time  anomalies  of  Cleary  and  Hales  (19SG),  both  comparison  yield  a 
correlation  coefficient  of  less  than  D.  1  where  a  perfect  correlation  is  1.0.  The 
correlation  between  Sengupta  (1979)  P  travel  time  station  anomalies  from  deep 
earthquakes  and  the  Wtt'SSX  mean  relative  amplitudes  presented  here  has  less 
than  a  0.5  correlation  coefficient.  This  result  indicates  that  while  P  wave  ampli¬ 
tudes  and  travel  times  may  correlate  over  an  area  of  MO4  km*  for  an  area  the 
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size  of  the  United  States  *  nearly  BOO  times  greater  -  P  wave  amplitudes  do  not 
correlate  statistically  with  travel  times.  This  in  turn  suggests  that  local  varia¬ 
tions  of  amplitude  may  be  influenced  by  elastic  focussing/defocussing  effects 
which  will  also  influence  local  travel  time  variations.  For  areas  of  10®  to  107  km2 
in  size  observed  P  wave  amplitude  variations  cannot  be  ascribed  wholely  and 
simply  to  elastic  focussing/defocussing  effects.  To  explain  the  amplitude  varia¬ 
tions  across  the  United  States  we  shall  have  recourse  in  energy  loss  mechan¬ 
isms  such  as  anelastic  attenuation  or  Q. 

Random  scattering  will  induce  variations  in  the  P  wavefront  as  it  pro¬ 
pagates.  The  seismic  station  LCN,  Longmire,  Washington  situated  on  the  flank 
of  the  volcano  Mt.  Ranier  exhibits  a  large  standard  deviation  in  measured  P 
wave  amplitudes,  extreme  azimuthal  amplitude  variability,  and  noticeable 
waveform  complexity  of  P  waves  relative  to  other  stations  in  the  United  States. 
These  characteristics  may  be  attributed  to  some  scattering  effect,  but  little  is 
gained  in  doing  so.  The  question  is  how  such  an  appellation  will  increase  our 
knowledge  as  to  what  is  actually  happening.  The  appellation  scattering  begs  a 
whole  series  of  further  questions:  where  is  it  occurring?  What  are  the  velocity 
perturbations?  What  are  the  size  of  the  scatterers?  Is  part  of  the  effect  ane- 
lastic?  How  does  it  affect  mean  travel  times  or  amplitudes?  How  is  it  related  to 
the  standard  deviation?  What  does  it  imply  in  regards  to  earth  processes?  What 
is  its  frequency  dependence?  How  important  are  scattering  losses  in  relation  to 
anelastic  attenuation  or  in  relation  to  elastic  focussing/defocussing  effects 
from  deterministic  structure?  None  of  these  are  easy  questions.  Some 
researchers  (e.g..  Aki,  1973,  Capon.  1974;  Berteussen  el  al.,  1975)  have  used  the 
Chernov  (19(30)  approach  to  random  media.  Butler  and  Ruff  (1930)  compared 
amplitudes  of  P  wave  secondary  arrivals  to  the  amplitude  of  the  first  cycle  of  P 
wave  motion  of  Russian  nuclear  explosions.  If  scattering  was  the  predominant 
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mcchanism  underlying  the  P  wave  amplitude  variations,  then  the  testable 
hypothesis  was  whether  low  amplitude  P  waves  showed  large  secondary  (scat¬ 
tered)  arrivals  while  high  amplitude  P  waves  exhibited  minor  secondary  arrivals. 
The  results  Butler  and  Ruff  (19B0)  found  were  that  the  secondary  P  wave 
arrivals  followed  in  relation  to  the  primary  P  wave  amplitude  and  that  the  sim¬ 
ple  test  of  the  scattering  hypothesis  did  not  fit  the  data. 

Scatter  in  the  amplitude  data  does  exist  in  Figures  6.  7,  and  B  -  yet  it  is 
unclear  where  along  the  path  between  source  and  seismometer  it  is  most 
important.  Obviously,  at  LON  the  receiver  region  is  important.  The  average 
standard  deviations  of  amplitudes  from  the  individual  Russian  test  sites  are 
lower  than  are  the  average  standard  deviations  of  amplitudes  measured  from 
collections  of  earthquake  sources.  However,  the  standard  deviations  of  ampli¬ 
tudes  Trom  the  northern  azimuth,  which  are  averaged  over  the  five  Russian  test 
sites,  have  values  comparable  to  the  average  standard  deviation  of  amplitude 
for  the  northwest  and  south-southeast  azimuths  comprised  of  earthquakes. 
This  shows  that  variation  among  the  source  regions  contribute  to  the  scatter  of 
P  waves,  while  P  wave  amplitudes  from  a  small  source  region  have  relatively  low 
scatter.  Scatter  is  an  important  phenomena  which  may  be  qualitatively  dis¬ 
cussed,  but  not  quantitatively  addressed  easily.  In  interpreting  the  mean  rela¬ 
tive  amplitudes  at  individual  seismic  stations  the  influence  of  scatter  in  relation 
to  other  seismic  amplitude  effects  must  be  considered  and  weighed,  but  a  more 
quantitative  evaluation  must  await  a  more  refined  seismic  experiment. 


llopping  P  wave  Amplitude  Variations  into  Qa 

The  density  and  elastic  properties  of  the  seismometer  site  enter  into  the 
energy  amplitude  relationship.  It  was  Tor  this  reason  that  site  corrections  for 
sediment  amplification  were  derived  earlier.  By  correcting  the  amplitudes  of 
seismometers  situated  on  sediments  to  values  appropriate  for  hard  bedrock, 
the  seismometer  site  effects  are  normalized  to  a  common  basis.  Note,  however, 
that  estimated  sediment  amplitude  corrections  are  squared  in  transformation 
to  energy.  Therefore,  the  uncertainties  in  the  true  corrections  also  squared. 
With  the  corrections  for  sediment  amplification  applied,  the  square  of  the  rela¬ 
tive  amplitude  variation  in  the  United  States  is  equivalent  to  the  variation  of 
relative  energy  content  in  the  P  wavefronts  recorded  by  the  seismic  stations. 

By  transforming  the  amplitude  variations  into  seismic  energy  variations, 
effective  Q  variations  naturally  follow.  Let  us  assume  for  the  moment  that  the 
energy  variations  are  due  to  energy  losses  from  attenuation.  From  the  defini¬ 
tion  of  the  quality  factor  Q  (e.g.  Anderson,  1967), 

LE/  E  =  2tt/  Q  (20) 

where  A E  is  the  energy  loss  per  cycle.  In  an  attenuating  medium  characterized 
by  the  dissipation  parameter  Q,  in  each  cycle  of  motion  the  energy  decreases 
by  a  factor  of  [l  -  2rr/  ^],  and  over,  n  cycles  the  decrease  is  a  factor  of 
[l  -2 -/  Q~\n .  Since  the  amplitude  variations  in  the  United  States  arc  relative, 
then  the  energy  variations  and  hence  the  energy  losses  are  relative.  Since  the 
amplitude  variations  were  developed  in  relation  to  the  geometric  mean  ampli¬ 
tude  fur  the  United  States,  then  the  energy  variations  are  in  relation  to  the 
mean  energy.  If  the  P  waves  at  each  seismic  station  i  have  incurred  energy 
losses  from  attenuation  due  to  a  local  Qt,  then  the  geometric  mean  amplitude  in 
the  United  States  has  incurred  an  energy  loss  from  attenuation  due  to 
average  Q 
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Ulustrative  example:  let  station  i  record  a  P  wave  amplitude  A,  and  station 
m  record  amplitude  An  with  both  stations  situated  on  the  same  material,  and 
let  the  P  waves  recorded  at  both  stations  have  the  same  apparent  period.  Then 
the  energy  £  at  station  i  is  related  to  the  amplitude  Al  by  £  -vt*,  and  at  sta¬ 
tion  m  the  same  relation  in  En  ~  Alii.  Since  the  two  seismic  stations  are 
situated  on  the  same  material  and  the  P  waves  recorded  have  the  same 
apparent  periods,  the  squared  amplitude  ratio  equals  the  energy  ratio: 

(A,/ An)2  =  £;/£•„  (21) 

Now  assume  that  the  energy  in  the  P  wavefronts  at  both  stations  i  and  m  were 
equal  to  £0  at  some  depth  and  that  energy  losses  were  incurred  in  propagation 
through  an  attenuating  region.  If  the  energy  loss  occurred  over  n  cycles,  then 

£  =£c[l  (22) 


£n  =r°[i  -  £-r  (23) 

where  Qx  and  Qm  are  the  average  quality  factors  controlling  the  attenuation  at 
station  i  arid  m,  respectively.  Substituting  (22)  and  (23)  into  (21)  and  cancel¬ 
ing  the  common  energy  E0.  a  relation  between  rr-lative  amplitude  and  relative  Q 


is  derived 
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Note  that  in  (24)  only  the  ratio  of  the  measured  P  wave  amplitude  enters. 
Since  only  the  ratio  is  important,  let  Am  -  1.  that  is,  station  m  records  ampli¬ 
tudes  at  a  level  equal  to  the  geometric  mean  amplitude  for  the  United  States. 
Then,  Qm  a  the  average  level  of  attenuation  beneath  the  United  Slates  over 
n  cycles  of  motion. 


In  mapping  the  variations  of  P  wave  amplitudes  into  variations  in  Qa 
beneath  the  seismic  stations  following  the  method  of  the  previous  section  all 
other  effects  influencing  the  P  wave  amplitudes  must  be  corrected.  Amplifica¬ 
tion  effects  for  seismometers  situated  on  sediments  relative  to  hard  bedrock 
sites  are  most  important  and  most  easily  recognized.  Elastic 
focussing/ defocussing  effects  are  second  in  importance,  but  their  specific 
influence  is  unknown  at  almost  all  seismic  stations  at  this  time.  Following 
amplitude  stability  results  in  a  previous  section,  amplitude  variability  among 
stations  of  a  region  averages  at  about  a  factor  of  1.25  around  the  regional 
mean,  suggesting  elastic  focussing/defocussing  effects  of  this  order  are  credi¬ 
ble.  Amplitude  variations  between  azimuths  average  at  a  factor  of  1.5  for  the 
WWSSN  stations.  Random  scattering  probably  increases  the  standard  dev.ation 
of  the  relative  amplitude  measurements,  but  the  effect  of  random  scattering 
upon  the  mean  relative  amplitude  is  unknown  at  this  time. 

In  transforming  the  variations  of  amplitude  into  variations  of  Qa,  only  sedi¬ 
ment  amplification  effects  v  ill  be  corrected.  Part  of  the  amplitude  variations 
are  probably  due  to  scatter.  In  mapping  the  variation  into  attenuation,  the 
resultant  values  must  be  considered  to  be  effective  or  apparent  Qa.  The  fun¬ 
damental  direction  is  to  obtain  Qa  variations  which  are  consistent  with  the 
observed  amplitude  variations.  The  standard  errors  of  the  mean  relative  ampli¬ 
tudes  will  be  used  to  estimate  the  uncertainty  of  the  effective  Qa  at  each  sta¬ 
tion. 

Geologic  site  information  W'as  listed  in  earlier  and  estimated  sediment 
amplitude  corrections  where  appropriate  are  listed  in  Appendix  Table  A22  for 
the  WWSSN  stations  in  the  United  States.  Seven  stations  out  of  29  required  sed¬ 
iment  corrections:  BKS,  LUB,  DAL,  FLO,  OXF,  AAM,  SHA,  and  BEC.  These  ampli¬ 
tude  corrections  neither  effect  the  highest  amplitude  station,  RKON,  nor  the 


lowest  amplitude  station,  GOL 


The  greatest  amplitude  ratio  in  the  data  set  is  between  RKON  and  GOL  - 
ratios  of  3.2  for  the  south-southeast,  5.2  for  northwest  azimuth.  The  amplitude 
value  of  RKON  to  north  will  be  discussed  in  a  moment.  Following  the  previous 
section  an  amplitude  ratio  of  5.2  is  equivalent  to  an  energy  ratio  of  27.  The 
energy  loss  of  GOL  relative  to  RKON  may  be  related  to  Qa  if  the  number  of  cycles 
over  which  the  energy  loss  occurred  is  estimated.  As  an  assumption  the  size  of 
the  attenuation  zone  must  be  estimated.  In  this  analysis  it  shall  first  be 
assumed  that  the  attenuation  variation  takes  place  in  the  upper  300  km  of  the 
earth.  This  data  set  has  little  resolution  of  depth,  and  the  300  km  range  was 
chosen  to  be  consistent  with  free  oscillation  Q  models  which  have  relatively 
higher  Q  below  this  depth  (Anderson  and  Hart,  1978;  Sailor  and  Dziewonski, 
1973).  The  apparent  lack  of  a  distance  trend  in  the  body  wave  attenuation  data 
between  30°A  and  85°A  further  suggests  attenuation  occurs  primarily  in  the 
upper  mantle  (e  g.,  Hart  and  Butler,  1977;  Sengupta  and  Toksoz,  1977).  The 
actual  attenuation  zone  could  occur  over  a  narrower  or  broader  depth  range 
with  a  corresponding  respective  decrease  or  increase  in  the  apparent  effective 
Qa.  The  basic  data  are  the  apparent  energy  losses,  which  are  uneffected  by  our 
choice  of  where  the  attenuation  occurs. 
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To  define  the  starting  estimate  of  Q*v  *  the  average  Qa  for  the  United 
States  -  the  geometric  mean  amplitude  y  was  calculated  for  each  azimuth, 
excluding  the  sediment  amplitude  corrected  stations  above  and  also  LON.  BOZ, 
and  RKON'  for  reason  of  high  azimuthal  amplitude  variability.  The  geometric 
mean  y  was  then  divided  into  the  observations  in  each  azimuth.  The  resulting 
baseline  is  then  nearly  free  of  sediment  amplitude  correction  uncertainties  and 
instability  due  to  the  azimuthal  highly  varying  stations.  After  this  renormaliza¬ 
tion  all  stations  were  included  in  the  further  analysis. 

Assuming  the  average  Qa  of  =  140  for  the  renormalized  baseline,  the 
variation  of  effective  Qa  for  the  seismic  stations  may  be  calculated  from  equa¬ 
tion  (24).  using  n  =  35  cycles  for  northwest  and  south-southwest  azimuth  and 
n  =  34  cycles  for  the  northern  azimuth.  The  effective  values  of  Qa  model  KMRl 
for  the  VfWSSN  stations  in  the  United  States,  derived  under  the  restrictive 
assumptions  mentioned,  are  listed  in  Table  13  for  the  three  azimuths  studied. 
The  effective  Qa  are  computed  from  the  mean  relative  amplitudes,  and  the 
ranges  of  effective  Qa  are  computed  from  the  standard  error  of  the  mean  rela¬ 
tive  amplitudes.  Using  the  =  140  yields  a  range  of  relative  Qa  values  which 
are  finite  and  positive  to  within  a  standard  error  in  the  mean  relative  amplitude 
measure  Tor  all  data  except  the  unique  observations  to  a  northern  azimuth  for 
RKON  and  HNME.  HNME  has  an  effective  Qa  >  1000  estimated  from  the  high 
amplitude  and  from  uncertainty  of  standard  errors  of  the  means  at  other  sta¬ 
tions.  If  Qy  equals  150  or  greater  certain  stations  -  RKON,  HNME.  DAL,  COH  - 
will  show  negative  effective  <?#,  implying  relative  amplification  in  relation  to  the 
energy  loss  differential  available  from  attenuation  over  34  to  35  cycles  C»i'  for 


*- 


Table  13 


Qa  Model  KMRl 


Northern 

Northwestern 

South- South eastern 

Azimuth 

Azimuth 

Azimuth 

Station 

5* 

Range* 

5* 

Range* 

Q* 

Range* 

(west) 

BKS 

166 

132-215 

157 

126-193 

125 

113-139 

COR 

166 

141-199 

262 

197-367 

916 

494-4100 

LON 

296 

167-661 

70 

60-81 

83 

69-93 

CSC 

147 

123-178 

118 

106-132 

153 

135-173 

0B2XV 

430 

312-563 

145 

134-158 

279 

209-401 

MSO 

147 

113-194 

116 

103-130 

130 

115-149 

DUG 

97 

65-113 

215 

184-255 

171 

145-203 

TUC 

127 

123-131 

84 

79-89 

104 

95-114 

BOZ 

451 

324-596 

239 

193-302 

69 

54-93 

ALQ 

77 

72-83 

91 

85-98 

102 

94-110 

COL, 

68 

82-77 

65 

62-68 

90 

87-94 

LUB 

122 

101-149 

91 

80-103 

203 

133-349 

JCT 

167 

161-215 

228 

150-360 

225 

81-4600 

DAL 

950 

533-1600 

140 

96-211 

415 

69-** ••• •••• 

RKON 

>1000 

•  •• 

956 

470-13500 

999 

385-  '* 

OXF 

122 

10B-140 

209 

160-285 

195 

89-677 

SHA 

161 

125-210 

99 

84-116 

AAM 

212 

174-256 

270 

12B-1170 

B9 

65-118 

ATL 

111 

100-125 

201 

170-270 

143 

120-173 

BLA 

145 

105-207 

346 

254-513 

180 

132-237 

SCP 

327 

188-702 

153 

131-179 

112 

106-117 

CEO 

463 

262-1000 

147 

120-182 

143 

115-151 

OCD 

111 

97-129 

133 

99-182 

90 

77-105 

WES 

88 

75-105 

250 

202-319 

104 

101-107 

HNME 

>1000 

Mil 

165 

139-197 

198 

162-246 

(cast) 

•  From  standard  error  of  the  mean  relative  amplitudes 

••  Qa  slightly  out  of  range  of  finite  positive  value,  implying  very  high 

•••  RKON  has  higher  frequency  content  to  the  northern  azimuth  and 
extremely  high  amplitude,  implying  very  high  Qa. 

••••  Estimated  from  uncertainty  of  standard  errors  of  the  means  at  other  sta¬ 
tions 


P  waves  may  be  less  than  140.  which  represents  an  upper  bound.  As  Qtv  -  140 
for  1  sec  P  waves  in  the  United  States  is  already  low  in  relation  to  free  oscilla¬ 
tion  and  surface  wave  Q  models  of  the  upper  mantle,  there  seems  to  be  little 
reason  to  chose  a  lower  value. 

The  mean  and  range  of  Qa  for  Qa  model  KMRl  naturally  follow  the  pattern 
of  the  amplitude  data,  with  corrections  included  in  Qa  model  KMRl  for  sediment 
amplification.  The  range  of  65  <  Qa  <  105  is  evident  for  all  azimuths  along  the 
Rock  Mountain  front  and  Rio  Grande  Rift  Zone  at  GOL  and  ALQ.  Effective 
9«>950  is  consistent  with  the  amplitude  observations  on  the  Canadian  shield  at 
RKO.N  and  to  the  north  of  Maine  at  HXME.  The  range  450  <  Qa  £  950  is  observed 
for  certain  azimuths  of  DAL  and  COR.  The  range  250  <  Qa  <  450  is  observed  at 
certain  azimuths  of  COR,  LON,  0B2NV,  B0Z,  DAL.  FLO.  AAM.  BLA,  SCP,  GEO.  and 
WES.  At  the  low  end  stations  with  Qa  <  100  for  certain  azimuths  include  LON, 
B0Z.  TUC.  ALQ,  GOL.  LUB.  AAM.  0GD.  and  WES. 

Tying  LRSM  Amplitudes  to  the  WWSSN  Array 

Having  mapped  the  variations  of  mean  relative  amplitudes  into  variations 
of  effective  Qa  for  WWSSN  stations,  we  now  proceed  to  tie  LRSM  amplitude  varia¬ 
tions  to  the  WWSSN  array.  Cleary's  (1967)  amplitude  study  of  LRSM  amplitudes, 
although  not  azimuthal  in  scope,  analyzed  more  than  fifty  LRSM  stations,  pub¬ 
lished  error  bounds,  and  used  a  consistent  amplitude  measure  of  the  P  wave 
first  motion.  Eighteen  out  of  twenty-two  earthquakes  measured  by  Cleary 
(19G7)  lay  to  azimuths  northwest  and  south-southeast  from  the  United  States, 
indicating  that  the  non-ezimuthal  Cleary  (1967)  LRSM  amplitudes  were  influ¬ 
enced  mostly  from  NW  and  SSE  azimuths  of  approach.  While  6  out  of  26  WWSSN 
stations  in  the  United  States  were  situated  on  sediments.  30  out  of  56  LRSM  sta¬ 
tions  were  situated  on  sediments.  Geologic  site  information  and  estimated  sedi¬ 
ment  amplitude  corrections  are  listed  in  the  Appendix  in  Table  A23.  The  WWSSN 
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*  SDCS  stations  HKON  and  HNME  are  located  at  sites  previously  occupied  by 
LRSM  instruments  used  by  Cleary  (1967).  No  other  site  locations  are  shared  by 
the  WWSSN  and  LRSM  systems.  However,  the  siting  of  LCNM  at  Las  Cruces,  New 
Mexico  within  the  Rio  Grande  rift  zone  is  similar  to  the  site  of  ALQ.  Albuquerque. 
New  Mexico.  In  tying  the  LRSM  to  the  WWSSN-SDCS  array,  the  relative  amplitude 
triangle  RKON-HNME-LCNM  is  matched  to  RKON-HNME-ALQ.  Three  amplitude 
ratios  are  formed  combining  the  three  stations  two  at  a  time.  For  Cleary  (1967) 
the  amplitude  ratios  are:  RKON/LCNM  =  2.95,  RKON/HNME  =  1.74,  and 
HNME/LCNM  =  170.  For  the  WWSSN-SDCS  amplitudes  presented  herein,  the 
northwest  and  south-southeast  azimuths  are  used  [NW,  SSE]:  RKON/ALQ  =  [3.1 1, 
2.74],  RKON/HNME  =  [1.76,  1.57],  and  KNME/ALQ  =  [1.77,  1.74].  Thus  the  rela¬ 
tive  amplitudes  of  the  LRSM  stations  RKON-HNME-LCNM  form  a  similar  triangle 
with  relative  amplitudes  of  the  WWSSN-SDCS  stations.  RKON-HNME-ALQ.  To 
obtain  the  best  least-squares  amplitude  match  for  RKON-KNME-LCNM  from 
Cleary  (1967)  to  RKON-HNME-ALQ  amplitudes  presented  here,  a  constant  log 
amplitude  of  (+  0.01)  was  added  to  each  log  amplitude  determination  of  Cleary 
(1967).  This  procedure  yields  mean  relative  amplitudes  of  RKON  2.00,  KNME 
1.15,  and  LCNM  0.6B  for  Cleary's  (1967)  LRSM  stations  which  compare  favorably 
with  northwest  and  south-southeast  [NW.  SSE)  azimuthal  determinations  for  the 
present  study:  RKON  [1.99.  2. 00].  HNME  [1.13,1.27],  and  ALQ  [0.64.0.73].  This 
shows  that  log  amplitudes  of  Cleary  (1967)  with  a  small  baseline  shift  of  (+  0.01) 
are  compatible  with  relative  amplitudes  determined  for  the  WWSSN  array.  From 
the  Cleary  (19G7)  study.  56  LRSM  stations  were  tied  to  the  WWSSN-SDCS  base¬ 
line.  While  the  amplitude  measurements  in  Evernden  and  Clark  (1S70)  and 
Booth  et.  al.  (1974)  were  not  limited  to  the  first  cycle  of  P  wave  motion,  the 
measurements  present  amplitude  variations  at  a  period  T  ■**  1  sec.  To  obtain  an 
amplitude  estimate  for  Qa  at  LRSM  stations  not  covered  by  Cleary  (1967),  the 


station  amplitude  anomalies  from  Evernden  and  Clark  (1970)  and  Booth  et.  al. 
(1974)  were  tied  to  the  Cleary  (1967)  study  through  common  stations.  Sedi¬ 
ment  corrections  are  estimated  from  Tables  11  and  12  where  necessary. 
Appendix  Tables  A24  and  A25  list  amplitude  and  Qa  data  for  LRSM  stations  in 
the  studies  of  Evernden  and  Clark  (1970)  and  Booth  et.  al.  (1974),  respectively, 
for  stations  not  covered  in  the  study  of  Cleary  (1967). 

Effective  Qa  for  the  LRSM  Array 

Applying  the  small  (+  0.01)  baseline  shift  to  the  Cleary  (1967)  log  ampli¬ 
tudes  and  correcting  the  amplitudes  of  stations  situated  on  low  velocity  sedi¬ 
ments,  the  LRSM  log  amplitudes  are  transformed  by  exponentiation  (101)  to 
amplitude  and  then  into  effective  Qa.  Table  14  lists  the  mean  and  range  of 
effective  Qa  for  the  LRSM  stations  of  Cleary  (1967),  where  the  range  is  derived 
from  the  standard  error  of  the  station  mean  amplitude.  As  in  the  case  of  the 
earthquake  sources  for  the  WWSSN  data,  equation  (24)  was  used  assuming 
n  =  35  cycles  and  Q^y  =  140.  The  mean  Q„  for  the  LRSM  stations  range  from  a 
maximum  978  at  RKON  to  a  low  of  64  for  CKBC,  Cache  Creek,  Eritish  Columbia. 
Stations  with  high  amplitudes  having  apparent  effective  Qa  &  200  are:  BDPA, 
Bedford,  Pennsylvania;  FSAZ,  Flagstaff,  Arizona;  HHND,  North  Dakota;  WNSD, 
South  Dakota;  PFM1,  Pickford,  Michigan  upper  peninsula;  RKON,  Red  Lake, 
Ontario  and  RYND,  Ryan,  North  Dakota.  With  exception  to  FSAZ  the  highest  Qa 
stations  lie  in  the  north  for  the  LRSM  data  set.  Stations  with  low  amplitude  and 
apparent  effective  9«^100  are:  CKBC.  Cache  Creek,  British  Columbia;  CPCL, 
Campo,  California;  CVTN,  Centerville,  Tennessee;  GDVA,  Grundy,  Virginia;  G0NB, 
Nebraska;  KGAZ,  Kingman,  Arizona;  MPAR,  Mount,  Arkansas;  MVCL,  Marysville, 
California;  TFCL,  Taft,  California;  W1NV,  Winnemucca,  Nevada,  and  WTTN.  Wart- 
burg,  Tennessee. 


NUMBER  OF  Qa  DETERMINATIONS 


KMR1  Qa  DETERMINATIONS:  130  TOTAL 


MEDIAN  140 
LOWER  25%  <95 


UPPER  300  KM  AVERAGE 


UPPER  25%  >200 


vO°  N^°  cfP  >1000 

APPARENT  EFFECTIVE  Qa 


Figure  15.  A  histogram  of  130  determinations  of  Qa  in  the  United  States 

and  southern  Canada  shows  the  distribution  Qa  for  Qa  model  K>!R1 . 
The  Qa  represent  average  values  over  34  to  35  cycles  of  vibration 
in  the  upper  300  km  of  the  earth.  The  median  is  Qa  -  140,  while 
25*1  of  the  Qfl  are  less  than  95  and  25*  are  greater  than  200. 
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Results  for  Qa 

A  summary  histogram  of  Tables  13  and  14  shows  in  Figure  15  the  distribu¬ 
tion  of  Qa  determinations  for  the  United  States  and  southern  Canada,  where 
numbers  of  observations  are  grouped  versus  Qa  values.  The  median  Qa  is  140 
for  the  130  determinations  of  Qa  from  Qa  model  KMRl  using  the  amplitude  data 
from  Cleary  (1967)  and  here.  This  median  Qa  of  140  equals  the  starting  value 
Qiv  =  140  chosen  for  WWSSN  stations.  The  equality  of  the  median  Qa  with  the 
starting  Q^y  —  140  indicates  that  the  sediment  amplitude  correction  is  a  proper 
procedure  and  that  the  Qa  for  sedimentary  sites  are  not  biased  with  respect  to 
hard  bedrock  sites.  Note  also  in  Figure  15  that  25"  of  the  determinations  are 
Qa  a  200  and  25%  have  Qa  <  95. 

The  apparent  effective  qa  determinations  of  Q  model  KMRl  are  plotted  in 


w'tf'  ifci'jpwjTM-  - xgp--  jigiy- 


Figure  16.  The  WWSSN  and  SDCS  stations  are  marked  by  solid  dots  at  the  sta¬ 
tion  location  -  see  Figure  1  for  comparisons.  Three  Qa  determinations  are  (usu¬ 
ally)  plotted  around  the  site  in  relation  to  the  three  azimuths  of  of  approach  in 
the  study  -  N.  N W,  and  SSE.  The  LRSM  stations  of  Cleary  (1967)  are  marked  by 
open  circles,  with  a  single  Qa  plotted  with  the  station  location.  Estimated  Qa 
from  the  studies  of  Evemden  and  Clark  (1970)  and  Booth  et.al.  (1974)  are  plot¬ 
ted  with  an  "X"  in  Figure  16.  The  effective  are  plotted  with  the  station  name 
code.  The  data  plotted  are  the  mean  values  from  Tables  13  and  14  and  from 
Appendix  Tables  A24  and  A25.  In  the  KMRl  Qa  model  the  amplitude  variations 
are  mapped  into  energy  losses  in  the  upper  300  km  of  the  earth  through 
attenuation  over  34  to  35  cycles  of  vibration.  The  effective  Qa  in  Figure  16  are 
replotted  as  symbols  in  Figure  17  divided  among  the  4  quartiles  of  Figure  15: 
low  Qa  <  95;  low-intermediate  Qa  <  140;  intermediate-high  Qa  fe  140;  and  high 
Qa  2  200.  Qa  values  below  140  are  surrounded  by  stippling  to  draw  the  low  Qa 
areas  together.  The  symbol  area  in  Figure  17  represents  roughly  the  lateral 
distance  the  seismic  ray  travels  in  the  upper  220  km  of  the  earth.  From  a 
depth  of  600  km  a  P  wave  travels  a  lateral  distance  of  about  385  km  when  arriv¬ 
ing  from  an  event  60*A  distant.  The  Qa  values  are  average  values  over  over 
some  lateral  distance. 

The  west  coast  of  the  United  States  shows  the  full  range  of  effective  Qa. 
Some  very  low  Qa  are  seen  in  California  in  the  Imperial  Valley,  Taft,  and  Marys¬ 
ville.  The  Great  Basin  has  both  high  apparent  Qa  at  0B2NV  in  the  Nevada  Test 
Site  and  northwest  of  DUG,  near  Provo,  Utah  and  low  Qa  at  W1NV,  Winnemucca 
and  MNNV,  Mina,  Nevada.  W1NV,  Winnemucca,  Nevada  is  near  the  Battle  Moun¬ 
tain  heat  flow  high.  High  Qa  is  seen  near  COR,  Corvallis,  Oregon  and  north  of 
LON,  Longmire,  Washington.  Low  Qa  extends  northwest  and  south-southeast  of 
LON,  across  the  Columbia  Plateau  and  eastern  Snake  River  plain  to  Yellowstone 
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Table  14 
Qa  Model  KMRl 


Station 

Q. 

Range* 

Station 

Qo 

Rarge* 

ARV.’S 

133 

100-172 

KCAZ 

B4 

69-109 

ATXV 

144 

104-235 

KNUT 

90 

AYSD 

90 

77-109 

LCNM 

95 

BS-104 

BDPA 

324 

170-4255 

MMTX 

101 

91-113 

BFCL 

140 

102-224 

MX  XV 

93 

95-102 

BIWV 

83 

77-90 

MPAR 

72 

6B-75 

BUQB 

117 

99-147 

Xft’CL 

72 

65-77 

CKBC 

64 

55-77 

XfZAR 

110 

65-145 

CPCL 

69 

62-96 

KDCL 

114 

95-154 

CTOK 

112 

92-143 

NGYi'S 

110 

25-130 

CVTX 

66 

59-75 

KPXT 

169 

114-330 

DHNY 

140 

123-163 

PFMI 

347 

193-1752 

DRCO 

82 

76-98 

PMVTY 

90 

93-95 

EYNV 

110 

85-154 

PT0R 

99 

90-110 

F.MUT 

123 

109-140 

RKCN 

975 

299-** 

FRMA 

148 

111-223 

RTNM 

93 

75- 125 

FSAZ 

212 

171-279 

RYXD 

347 

176-40,000 

cdva 

63 

73-96 

SE.'.fX 

144 

125-169 

G1MA 

163 

113-299 

SJTX 

93 

94-  r:  5 

GOXB 

66 

72-105 

SSTX 

105 

94-119 

GYTX 

136 

120-159 

ST  XV 

149 

105-253 

KBOX 

110 

99-123 

TFCL 

70 

65-75 

KHND 

324 

170-4255 

TKVA 

101 

79. 1 39 

HKWY 

169 

122-273 

TUPA 

129 

105-167 

HUD 

103 

9-112 

IV]  XV 

79 

73-94 

KMBC 

169 

122-273 

VTNSD 

199 

157-233 

KXME 

160 

142-209 

vrrrx 

76 

63-55 

HTMX 

195 

135-352 

WWL’T 

101 

£  ]  3*r 

•  »* 

MDS[N] 

235 

194-201 

BEC[X] 

66 

54-79 

RCD[NJ 

162 

1 3~- 1 24 

ussMm 

100 

85-110 

flo[  xv;  3 

321 

255-420 

135mm 

55 

42-G9 

Fi,o[ssri 

299 

225-352 

‘From  standard  error  of  mean  relative  amplitudes 

’*  Qa  slightly  out  of  range  of  positive  finite  value,  implying  very  (ugh  *'a 
‘"The  following  WViSSN  stations  from  Butler  and  Ruff  (19G0)  hav  e  []  a:  h  /  Tj 

determinations. 


Figure  16.  Determinations  of  Qq  for  WUSSN  and  LRSM  seismic  stations  are 

plotted  in  relation  to  physiographic  provinces.  Qa  model  KXRl 
maps  amplitude  variations  into  apparent  Qa  representing  energy 
losses  in  the  upper  3DD  km  of  the  earth.  Azimuthal  Qa  determinations 
at  stations  with  solid  dots  are  plotted  in  relation  to  azimuth 
from  the  station.  The  Qa  values  are  determined  in  relation  to  an 
average  Qa  ■  140. 


A  >  140 
^7-.<  140 
r&  <  95 


Figure  17.  The  variation  of  Qa  is  shown  with  four  symbols  representing  ranges 
of  Qa.  Lower  Qq  areas  are  stippled  and  lower  Qa  values  are  plotted 
as  triangles  or  circles.  Higher  Qq  values  are  plotted  as  solid 
triangles  or  diamonds.  The  motion  of  the  North  American  plate 
for  20  m.y.  is  shown  relative  to  the  underlying  mantle,  the  two 
arrows  indicating  the  uncertainty  of  direction. 
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While  no  data  are  available  for  the  central  Rockies,  low  Qa  are  evident  along  the 
Rocky  Mountain  front,  the  Rio  Grande  Rift,  the  Great  Plains  north  of  LUB,  Lub¬ 
bock,  Texas,  and  in  the  southern  basin  and  range  at  TUC,  Tucson,  Arizona.  With 
the  exception  of  FSAZ,  Flagstaff.  Arizona  the  edges  of  the  Colorado  plateau  are 
characterized  by  low  Qa. 

The  central  United  States  and  Canada  has  among  the  highest  effective  Qa 
at  RKON.  Red  Lake,  Ontario  and  DAL.  Dallas,  Texas.  The  northern  Great  Plains 
has  a  trend  of  intermediate  high  in  eastern  Montana  increasing  eastward  to 
high  Qa  in  North  Dakota  and  the  highest  Qa  at  RKON  on  the  Canadian  shield. 
Intermediate  Qa  occur  across  Nebraska,  South  Dakota,  Minnesota,  and  Wiconsin. 
The  Ozark  plateau  shows  high  Qa  to  the  northeast  and  low  Qa  to  the  south.  The 
low  Qa  at  MPAR,  Mountain  Pine,  Arkansas  lies  very  near  to  the  hot  springs  of 
Arkansas.  Intermediate  to  high  Qa  are  found  in  central  Texas,  with  a  low  Qa  in 
south  Texas  at  SJTX,  San  Jose. 

The  United  States  east  of  the  Mississippi  River  displays  a  range  of  high  to 
low  apparent  effective  Qa  observed  in  the  central  and  western  United  States. 
High  to  intermediate- low  Qa  are  seen  at  SHA,  Spring  Hill.  Alabama  along  the  Gulf 
coast  and  in  Mississippi  embuyment  at  QXF,  Oxford,  Mississippi  High  Qa  is  found 
from  south  of  AAM,  Ann  Arbor,  Michigan,  then  east  to  BLWV,  Beckley,  West  Vir¬ 
ginia  and  GDVA.  Grundy,  Virginia,  extending  through  Tennessee  to  northeastern 
Mississippi  and  bounded  to  the  west  at  FLO,  Florissant,  Missouri  near  Saint 
Louis.  Although  the  connection  between  these  low  Qa  stations  may  be  more 
apparent  than  real,  the  individual  low  amplitudes  are  well  determined.  Inter¬ 
mediate  and  high  Qa  occur  in  the  eastern  Appalachias  and  the  Piedmont.  The 
high  apparent  Qa  raypath  northwest  of  BLA,  Blacksburg,  Virginia  passes 
beneath  the  low  amplitude  station  BLWV,  Beckley,  West  Virginia  at  a  depth  of 
nearly  200  km.  Intermediate  and  low  Qa  are  observed  in  southern  New  England. 


High  and  intermediate-high  effective  Qa  are  seen  at  HNME,  Houlton.  Maine  and 
to  the  northwest  of  WES,  Weston,  Massachusetts  near  Boston.  An  intermediate- 
low  apparent  Qa  is  found  for  BUQB,  Buckingham,  Quebec  on  the  Grenville  pro¬ 
vince  of  the  Canadian  shield.  A  low  Qa  is  observed  at  Bermuda  in  the  Atlantic 
Ocean  east  of  the  United  States. 

The  variation  of  the  mean  Q  in  the  upper  300  km  of  the  mantle  within  the 
continental  U.S.  can  also  be  represented  in  the  form  of  a  contour  map  of  Qa.  Of 
course,  the  station  coverage  of  the  U.  S.  is  such  that  several  areas  are  not  well 
sampled,  however,  if  the  azimuthal  variation  of  Q„  at  the  WWSSN  stations  is 
taken  into  account,  the  inference  of  the  Qa  variation  is  somewhat  more  con¬ 
strained.  In  any  case.  Figure  18  shows  Qa  contours  based  on  all  the  data, 
including  the  azimuth  variations  estimated  at  the  WWSSN  stations  (black  dots). 
While  there  are  uncertainties  in  this  representation,  it  nevertheless  provides  a 
useful  summary  diagram  incorporating  all  the  known  results. 
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III.  Anelasticity  of  the  Lower  Mantle  and  Source  Excitation  Inferred 
from  Low  Frequency  Free  Oscillations. 

INTRODUCTION 

In  recent  years,  the  availability  of  improved  quality  seismic 
records  has  stimulated  research  efforts  directed  at  determining  the 
anelastic  behavior  of  the  earth.  Obtaining  refined  estimates  of  the 
Q  of  the  earth  free  oscillations  has  been  one  of  the  primary  goals  of 
this  research.  Research  efforts  directed  at  measuring  the  Q  of  free 
oscillations  have  included  Alsop  et  a l.  (1961);  Slichter  (1967); 
Norwoozi  (1968,  1974);  Dratler  et  al.  (1971);  Sailor  and  Dziewonski 
(1978);  Buland  and  Gilbert  (1978);  and  Stein  and  Geller  (1970b) 
among  others.  Anderson  and  Hart  (1970)  very  nicely  summarize  the 
results  of  these  research  efforts.  Because  of  these  efforts,  many 
reliable  estimates  of  Q  for  all  but  the  longest  period  (lowest  order) 
free  oscillation  modes  are  now  available.  The  reliable  measurement 
of  Q  for  the  lowest  order  modes  has  suffered  from  the  relative  scar¬ 
city  of  high-quality  very  long  period  (VLP)  seismic  records.  How¬ 
ever,  this  is  not  the  only  obstacle  to  Q  estimation  for  these  modes 
Every  earth  free  oscillation  with  non-zero  angular  order  is  subject 
to  a  spectral  splitting  of  its  resonance  peak  that  results  from  the 
earth’s  rotation,  slightly  aspherical  shape,  and  lateral  hetero- 
genieties  (the  mechanical  analog  of  the  Zeeman  effect).  At  very 
long  periods,  the  spacing  of  the  split  members  (singlets)  of  a  given 
free  oscillation  is  such  as  to  cause  an  interference  that  manifests 
itself,  in  the  time  domain,  as  a  beating,  and,  in  the  frequency 
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domain,  as  a  deformation  of  spectral  peaks.  Stein  and  Geller 
(1078b),  in  an  experiment  to  evaluate  the  effects  of  this  splitting  on 
commonly  used  methods  of  measuring  attenuation,  demonstrated 
that  the  splitting  can  introduce  substantial  inaccuracies  in  the 
results  of  such  measurements  One  such  "conventional"  method, 
that  of  estimating  the  rate  of  decay  of  a  time  series,  is  thrown  off 
by  the  time  domain  beating,  while  another  method,  that  of  measur¬ 
ing  the  width  of  a  resonant  spectral  peak,  is  thrown  off  by  the 
interference  of  the  neighboring  singlet  peaks.  Reliable  measure¬ 
ment  of  Q  at  these  long  periods  requires  techniques  that  take  into 
account  the  spectral  splitting.  Stein  and  Geller  present  such  a  tech¬ 
nique.  Making  use  of  assumed  source  properties,  they  model  the 
excitation  of  the  free  oscillation  for  several  values  of  Q  and  then 
compare  the  envelope  functions!  of  the  resulting  time  series  with 
the  envelope  of  the  actual  time  series  band-pass  filtered  about  the 
frequency  range  of  the  free  oscillation  under  consideration.  The 
value  of  Q  that  corresponds  to  the  envelope  that'  most  closely 
approximates  the  actual  envelope  is  the  value  that  they  present  as 
the  Q  determined  by  their  method,  in  this  study  (Schnapp, 
1983),  I  present  a  method  of  measuring  Q  that,  also, 
involves  modeling  the  split  normal  mode.  However,  I 
deviate  from  the  approach  taken  by  Stein  and  Geller  in 
that  the  modeling  is  done  independently  of  any 

t  The  definition  of  the  envelopeg  function  that  I  use  throughout 
this  paper  is  E(t )=[/(*  )*+/'(*)*  J  where  E(t)  is  the  envelope  func¬ 
tion,  /(f)  is  the  time  series,  and  f(t)  is  the  quadrature  of  the  time 
series,  defined  as  the  inverse  fourier  transform  of  sign(u)i  f  (u) 
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knowledge  of  the  source,  and  the  comparison  or  ’fitting’  is  done  in 
the  frequency  domain.  Further,  estimations  of  the  source  depen¬ 
dent  excitation  amplitudes  are  obtained  in  addition  to  estimates  of 
Q  for  each  mode  inverted,  giving  the  inversion  results  potential 
value  in  an  inversion  for  source  properties. 


-86- 


THEORETICAL  DEVELOPMENT 

The  equation  of  motion  governing  the  oscillations  of  an  elastic, 
rotating  and  self-gravitating  body  is  (Dahlen  1968) 

AS  =  p0^S  (2  1) 

•where  A  is  an  operator  defined  by: 

AS  =  I1S+  2ipc0x  S 

and 

+  p1V(p0+y()  +  V[Sp07(^o+^)l 

-E-V[S(Vt0)]  +  9(SVtc) 
where  the  following  definitions  apply: 

Po  the  equilibrium  density  distribution 

p,  =  -V  (p0S)  is  the  change  in  density  due  to  displacement 

Po  is  the  equilibrium  gravitational  potential 

V*pj  a  4nGp, 

i>{r)  s  -I[n*r*  -  (fi  r)2]  is  the  rotational  potential 
Eis  the  elastic  stress  tensor 
t0  =  To-^tr  T0)I  is  the  static  stress  deviator 
T0  is  the  static  stress  field  tensor 

Thus,  n  defines  the  effects  of  gravity  and  the  relationship  of  stress 
tr.?  strain  in  the  body 

Under  the  assumption  that  the  effects  of  rotation  and  devia¬ 
tions  of  the  earth  from  radial  symmetry  and  lateral  homogene  ity 
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are  small,  it  is  possible  to  deal  with  these  effects  a?  perturbations 
to  the  equation  of  motion  of  a  spherical,  non-rotating,  radially  syrr - 
metric  earth  model.  Dahlen  (1968)  applies  perturbation  theory, 
and,  by  taking  advantage  of  the  hermitian  character  of  A,  he 
derives  an  expression  for  the  first  order  perturbation  to  the  eigen- 
frequency  due  to  the  effects  of  the  earth’s  ellipticity  ar.d  rotation. 
He  shows  that,  to  zeroth  order,  the  perturbed  eigenfunctions  are 
the  some  as  the  eigenfunctions  of  the  unperturbed  equation  of 
motion.  Because  A  is  hermitian,  the  eigenvalues  of  Eqn  2  1.  -J.  are 
real,  so  that  the  normal  mode  solutions  are  unattenuated  in  time 
Following  Liu  and  Archambeau  (1975),  the  effects  of  ar.eiasii- 
city  for  an  isotropic  body  can  be  included  in  the  equation  of  motion 
by  adding  an  anelastic  operator  Z  to  Tl,  where,  for  an  isotopic 
material  with  complex  elastic  contants  X*  and  /A 

(  n  ♦  Z  )S  =  -(X*  +2m*  )V(VS)  ♦  (f  *(7xS) ) 

-<7A*)(vs)-(V)(?s+sr)  {2  2) 

If  we  approximate  the  anelastic  behavior  of  the  body  by  a  Ke' -in- 

Voigt  solid  model,  II  is  the  real  part,  and  Z  is  the  imaginary  pe-t  of 

the  expression  on  the  right  side  of  Eqn  2.2.  If  we  assume  the  s.  me- 

what  more  general  standard  linear  solid  approximation  for  ire 

body,  Z  takes  on  a  real  component  Liu  and  Archambeau  ( *  9~5) 

compute  the  complex  perturbation  to  the  eiger.fr a quency  of 

toroidal  normal  modes  for  the  Kelvin-Voigt  solid  approximate-  ar.d 

in  a  subsequent  paper  (1976),  they  do  the  same  for  a  standard 

linear  solid  approximation  They  determined  that  the  corrections 
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to  the  real  part  of  the  eigenfrequency  for  the  Kelvin -Voigt  care 
were  on  the  order  of  0.1  percen*  which  was  at  the  limit  of  observa¬ 
tional  accuracy  at  that  time,  while  in  the  more  general  standard 
linear  solid  case,  the  corrections  were  considerably  more  signifi¬ 
cant,  between  2  and  5  percent.  In  either  case,  the  effect  of  Z  is  Lo 
make  the  operator  A  non-hermitian  so  that  u  takes  on  en  ir?.i=...arj 
component.  The  result  is  that  the  solutions  to  the  equation  of 
motion,  S(r  .*.?;«).  transformed  into  the  time  domain,  take  the  :crm 

*(*•.*.*;<)  *  f<0*‘ 

where  the  subscript  k  is  an  index  that  distinguishes  the  individual 
solutions  and  0*  and  yk  are  the  real  and  imaginary  parts  of  the  per¬ 
turbed  eigenfrequency.  From  this  expression  for  the  solutions,  it  is 
immediately  evident  that  the  effect  of  the  perturbation  is  to  intro¬ 
duce  the  damping  term  e-’’*'  into  the  solutions  In  terms  of  the 
quality  factor,  Q,  the  solution  takes  the  form 

pit)  =  A (2  5' 
It  is  this  form  of  the  solution  to  the  equation  of  motion  of  er.  isi  tro¬ 


pic,  anelastic,  rotating,  elliptical,  and  self-gravitating  earth  model 
that  forms  the  foundation  for  the  method  of  inversion  presen  .e_  .n 


the  following  discussion. 

The  inversion  technique  for  determining  the  Q  of  the  lc..  cr*.er 
free  oscillations  consists  of  finding  a  suitably  constrained  mode! 
that  has  a  complex  spectrum  that  is  closer  to  the  actuel  spec. rum. 
in  a  least  squares  sense,  than  any  other  similarly  co..r.. 
model  The  model  is  an  ideal  damped  multi-mode  oscillator  *■  uh 
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resonances  at  frequencies  that  correspond  to  the  singlet  frequen¬ 
cies  of  the  free  oscillation.  It’s  time  dependence  is  simply  that  of 
Eqn.  2.3.  The  mathematical  form  of  the  model,  in  the  time  domain, 


is: 


F(t )  *  Real 


(2  4) 


where  the  subscript  m  denotes  the  singlet  of  azimuthal  order  m 
Am  is  the  complex  excitation  amplitude  corresponding  to  the  m"1 
singlet  whose  magnitude  and  phase  are  determined  by  both  the 
source  properties  and  the  lag  between  the  time  when  the  actual 
excitation  occured.  and  the  reference  time  of  the  model  (the  time 
Is o  in  Eqn  2  4).  O,  is  the  angular  frequency  of  oscillation  of  the 
m1*  singlet  and  the  Q  is  the  Q  of  the  free  oscillation.  In  an  effort  to 


minimize  the  number  of  independent  parameters  in  the  inversion,  I 


have  assumed  that  Q  is  sub  tantially  constant  over  the  frequency 
range  of  the  multiplet. 

Eqn.  2.4  defines  the  time  domain  representation  of  a  mod?:  cf  a 
single  normal  mode  multiplet  (radial  order  n  and  angular  order  l 
are  fixed)  In  actual  application,  such  a  multiplet  is  isolated  from 
other  multiplets  by  the  application  of  a  narrow  bar.d-pass  filter  to 
the  data.  The  width  of  the  filter  is  chosen  so  that  all  spectral  sam¬ 
ples  that  display  significant  signal  content  as  compared  to  the 
ambient  noise  level  are  included  in  the  pass-band  Isolatatirr.  of 
the  •  5,  and  „S,  multiplets  from  other  multiplets  is  insured  in  this 
way  because  of  the  fact  that  at  the  long  periods  corresponding  to 


-90- 


these  modes.  neighboring  multiplets  are  well  separated  in  fi  e- 
quency. 

The  Fourier  transform  (abbreviated  FT)  of  F(t)  in  Eqn  2  ~  is 

(2  5) 

MI-I 

where  fm  is  the  complex  function  that  corresponds  to  the  FT  of  the 
exponents  in  Eqn  2  4  The  inversion  technique  involves  find.ng  .he 
model  whose  complex  spectrum,  F,  minimizes  the  difference  '.  no¬ 
tion: 

A  =  |^(«)-Fo(«)  (  du 

where  Fc  is  the  spectrum  of  the  data  to  be  inverted  S.r.ce  the  spec¬ 
tra  that  I'll  be  dealing  with  are  not  continuous  but.  rather, 
discretely  sampled,  a  more  appropriate  form  for  the  difference 
function  minimized  by  the  inversion  is 

A  =  £  !  /‘U  )-/-„(.*)  2  Ay 

*«i  '  1 

where  N  is  the  number  of  sampled  frequencies  .  end  1-  is  the 
spectral  sampling  interval 

An  important  property  of  Eqn  2.5  is  that  f  is  linear  in  the  c*-n- 
plex  excitation  amplitudes  Am  So.  with  given  single4,  frequencies 
Om.  ^  exhibits  a  nonlinear  dependence  on  only  cr.e  unhr.c-  n  3 
The  linearity  of  Eqn  2  5  in  all  unknowns  except  for  Cf  perm  the 
application  of  a  simple,  fast  inversion  method  for  both  the  3  and 
the  An's  simultaneously  This  technique  involves  solving  a  linear 
least-squares  problem  in  the  unknown  s  for  fixed  ve’ces  r'  .?  ?\ 
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comparing  the  resulting  values  of  the  difference  function,  the  Q 
which  minimizes  A  can  be  determined 

At  this  point,  it  is  useful  to  formulate  the  problem  in  vector 

form  Let  £  be  the  vector  J  P(uk),  *=0,1 . .V  J  and  A  be  the  vector 

. I  ]  Since  P  is  linear  in  the  An's,  Eqn  2  5  can  be  : Tit- 

ten  in  vector  form  as 

£  =  M d  (25) 

where  ■  is  the  matrix  whose  elements  are  (from  Eqn.  2  5) 

Jtfvs/t(uj.9.nt)  For  fixed  Q  and  singlet  frequencies  fl,-.  mir.:m  ?ir.c 

the  difference  function  expressed  in  vector  form  as 

A  =  I  £-£ o  !*  =  (£-£,)'  (Z~£ c) 

becomes  a  problem  in  determining  the  least-squares  inverse  of  Eon 
2  6,  that  is.  in  solving  the  vector  equation 

M*H4  =  M'£c 

for  A  UtfM  is  non-singular  the  solution  can  be  vritten  p.s 

A  =  (lTll)-,ll,£c  (? 

Thus,  given  a  sei  of  observed  complex  spectral  amplitudes. £z  er.d  a 

set  of  singlet  frequencies,  f^,  for  any  one  value  of  Q.  the  elements 
of  M  can  be  computed  and  a  least-squares  closes'-/:' P-.p  set  of 
excitation  amplitudes  A  can  b*  determined  by  computing  the  pro¬ 
duct  of  the  matrix  (M*ll)~lll>  with  the  vector  of  the  actual  complex 
spectrum.  £0  Repeating  the  process  for  several  values  cf  >?.  the 
relative  fits  of  the  resulting  cXosesf-/tllmg  models  to  the  actual 
spectrum  (as  measured  by  A)  can  be  compared  to  f’r.d  the  ~  :h:.t 
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minimizes  A  This  Q  and  associated  set  of  excitation  amplitudes  ! 
refer  to  as  the  estimated  or  invented  Q  and  associated  excitation 
amplitudes  The  model  generated  by  this  Q  and  associated  excita¬ 
tion  amplitudes  1  refer  to  as  the  best-fitting  model. 

Throughout  the  preceding  discussion  the  singlet  frecuer.cies 
were  assumed  to  be  known  Uncertainty  in  the  singlet  frequencies 
can  be  taken  into  account  by  allowing  them  to  vary  w.ih:n  their 
respective  uncertainty  limits  over  the  course  of  repeated  applica¬ 
tions  of  the  previously  described  inversion  process  to  a  particular 
data  set  Each  set  of  singlet  frequencies  will  have  associated  v.-ith  it 
a  best-fitting  model  and  corresponding  inserted  Q  Uncertainty  in 
the  Q  that  results  from  an  uncertainty  in  the  singlet  frecyuer.ries 
can  then  be  estimated  by  statisl ieally  eva'uating  1.hi«  set  of  invented 
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APPLICATION  TO  REAL  DATA 


For  a  first  attempt  at  applying  the  Q  inversion  technique  to  a 
real  data  set,  I  obtained  the  highest  quality  long-period  seismogram 
that  I  could  find,  a  500  hour  length  of  the  Isabella  strain  record  of 
the  May  22,  1960  great  Chilean  earthquake.  The  raw  time  series 
had  a  sampling  interval  of  one  minute  and  contained  32768  points 
I  pre-processed  the  raw  data  by  de-glitching  and  low-pass  filtering 
(to  inhibit  aliasing)  the  time  series.  This  allowed  me  to  decimate 
the  time  series  by  a  factor  of  eight  so  as  to  obtain  a  somewhat  more 
managable  length  of  4096  points  Fig  3.1  displays  the  pre- 
processed  seismogram  I  then  took  a  280  hour  long  section  of  the 
time  series  starting  at  28.4  hours  after  the  main  shock  and  applied 
a  Fast  Fourier  Transform  to  it  to  obtain  a  discretely  sampled,  finite 
length  complex  spectrum  Fig.  32  displays  the  modulus  of  the 
resulting  spectrum  Many  of  the  low  order  free  oscillations,  espe¬ 
cially  the  0S2,  0S3,  0S4,  and  0Ss  modes,  are  quite  evident  in  the  spec¬ 
trum  The  0T2.  oT3,  and  tT4  modes  are  not  quite  so  apparent 

In  order  to  apply  the  inversion  routine  to  this  data  set,  1 
required  estimates  of  the  splitting  parameters  that  define  the 
•inglet  frequencies  of  the  free  oscillations.  Buland,  Berger,  and  Gil¬ 
bert  (1979)  obtained  refined  estimates  of  the  singlet  frequencies 
for  0Sf  and  0S3  by  stacking  spectra  from  the  IDA  records  of  the  1976 
Sumbawa  Indonesia  earthquake.  From  these  frequencies,  they  com- 


Isabella  Strain  -  May  22,  1960  Chilean 


Fic  3  I  Pre-processed  (de-glitchcd  and  low-pass  filtered  to  inhibit  aliasing)  seismogram  of 
the  May  33.  I960  Chilean  earthquake.  Amplitude  is  relative  scaled  to  1  and  time 
is  in  hours  after  the  main  shock. 


Isabella  Strain  Spectrum  -  I960  Chilean  E^Q. 
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puted  the  three  splitting  parameters  u0,  / .  and  g,  that  define  the 
singlet  frequencies  via  the  relation  Um  =  oo+fm+gm*,  where  m  is  the 
azimuthal  order  of  the  singlet.  Table  1  lists  these  splitting  parame¬ 
ters  and  associated  singlet  frequencies  for  052  and  05s.  All  fre¬ 
quency  values  presented  in  this  thesis  are  in  units  of  cycles /minute 
(cpm)  unless  otherwise  noted.  Note  that  in  order  to  apply  the 
theory  developed  in  chapter  2,  these  values  must  be  converted  to 
angular  frequency  (e.g.,  radians/second). 


TABLE  1 

Assumed  Splitting  Parameters  and 
Associated  Singlet  Frequencies 


Mode 

Splitting 
Paras  (ccml 

Azimuthal 
_ Order _ 

Singlet 

Freo  (cDm) 

«0  =  .018588 

CM 

i 

II 

£ 

.018001 

-1 

.018280 

052 

/  =2.771xl(r' 

0 

.018568 

g  *  -3.253x1 0-8 

1 

.018842 

2 

.019109 

m.  =  -3 

.02*693 

u0  =  .028125 

-2 

.027845 

-1 

027989 

oS9 

/  =1.3l6xl0-4 

o 

.028125 

1 

.028252 

g  =  - 4.078x10-* 

2 

.028372 

3 

_ mm _ 

Using  these  splitting  parameters  as  the  assumed  splitting 
parameters  (the  splitting  parameters  used  to  compute  the  singlet 
frequencies  required  by  the  inversion  routine),  band-limited  sec¬ 
tions  of  the  complex  spectrum  in  the  0S2  and  0-?s  frequency  ranges 
were  processed  through  the  inversion  routine.  The  width  of  the 
spectral  bands  were  chosen  so  as  to  minimize  the  noise  power  in 
the  band  while  retaining  all  spectral  points  that  have  a  significant 
signal  component.  For  0S2,  the  band  selected  extended  from  .0170 
to  .0193  cycles  per  minute  For  0S3,  the  band  extended  from 
.027575  to  028575  cycles/min.  For  these  modes,  the  inversion  rou¬ 
tine  determined  a  single  minimum  in  the  difference  function  A  (see 
chapter  2)  for  a  range  of  Q‘ s  extending  from  50  to  2000.*  The 
resulting  inverted  Q's  were  399  for  0S2  and  358  for  „53.  The  associ¬ 
ated  complex  excitation  amplitudes  are  presented  in  Table  2.  In 
the  table,  the  amplitude  values  corresponding  to  a  given  multiple! 
are  all  scaled  by  the  same  factor  so  that  the  root-mean-square  of 
the  values  is  unity. 

In  order  to  evaluate  qualitatively  the  results  of  this  application 
of  the  inversion  routine  to  the  Chilean  data,  1  generated  some  plots 
that  display  a  comparison  of  the  actual  data  to  the  best-fitting 
model  Figs.  3.3  and  3.4  show,  for  each  mode,  a  time  domain  com¬ 
parison  of  the  envelopes  of  the  inverse  Fourier  transform  of  the 

•  I  used  the  criterion  of  a  single  minimum  in  A  for  $’s  between  50 
and  700  to  determine  which  modes  might  be  reliably  invertable 
eSg  and  3S3  were  the  only  modes  that  met  the  criterion 


Q=314 


comparison  as  displayed  in  Fig  3  3  for  o*^-  The  frequency  range 
,027:>7r)  -  .020fi7S  cpm. 


trum)  and  the  spectrum  of  the  best-fitting  model  in  the  frequency 
range  of  the  modes.  The  crosses  correspond  to  the  discrete  sam¬ 
ples  of  the  actual  spectrum,  and  the  continuous  curve,  to  the 
inverted  spectrum.  The  dashed  lines  indicate  the  location  of  the 
singlet  frequencies  used  in  the  inversion. 

The  plots  show  that  the  inversion  routine  is  yielding  model 
parameters  that  generate  envelope  functions  and  amplitude  spec¬ 
tra  that  fit  the  data  reasonably  well.  The  envelope  comparisons 
show  that  it  is  the  signal  in  the  earlier  portions  of  the  record  that  is 
primarily  being  fit  as  it  should  be,  while  the  noise  in  the  latter  por¬ 
tions  are  not.  The  0S3  envelope  fit  appears  to  be  quite  close  in  the 
first  150  hours.  The  closeness  of  the  fit,  however,  is  an  artifact  of 
the  large  number  of  inversion  parameters,  15  (7  complex  excitation 
amplitudes  and  <?),  relative  to  the  number  of  data  points,  17  real 
and  17  imaginary  components  of  the  spectrum.  The  number  of 
degrees  of  freedom  (=  #  of  data  points  -  #  of  parameters)  is  there¬ 
fore  only  21.  Note  that  for  0Sz.  with  37  degrees  of  freedom  (11 
parameters  and  48  data  points),  the  envelope  fit  is  not  as  good 
Thus,  the  apparently  good  0S3  envelope  fit  should  not  be  used  as  an 
indicator  of  especially  reliable  inversion  results.  The  spectral  com¬ 
parisons  also  show  that  the  spectral  points  that  contain  the  major- 
ity  of  the  information  on  the  signal,  i.e.,  the  highest  peaks  in  the 
spectrum,  are  fit  the  best.  The  comparison  displays  this  better 
than  the  9Sa  one.  Note  that  in  the  05e  comparison,  the  curve  show¬ 
ing  the  spectrum  of  the  best-fitting  model  falls  consistently  below 
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while  the  inverted  spectrum  is  plotted  as  a  continuous  curve.  The  vertical 
dashed  lines  indicate  the  locations  of  the  singlet  frequencies  used  in  the  inver¬ 
sion  Amplitude  is  relative  scaled  to  t  and  frequency  is  given  in  cycles/minute. 


Fig  3  6  Same  type  of  comparison  as  displayed  in  Fig.  3.5  for 
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the  actual  spectral  amplitudes.  This  is  a  result  of  the  fact  that  the 
inversion  fits  the  complex  spectrum,  i.e.,  the  real  and  imaginary 
parts,  while  the  plot  displays  the  modulus  of  the  spectrum.  Since 
the  modulus  is  the  sum  of  the  squares  of  the  real  and  imaginary 
parts,  if  the  best-fitting  spectrum,  on  the  average,  deviates  from 
the  actual  spectrum  equally  in  both  the  plus  and  minus  directions, 
the  modulus  will  tend  to  deviate  more  in  the  minus  direction,  thus 
causing  apparent  bias  in  the  plot.  As  with  the  envleopes,  the  CS3 
spectra  of  the  best-fitting  model  and  the  actual  data  compare  quite 
well.  Once  again,  this  is  likely  an  artifact  of  the  large  number  of 
inversion  parameters  relative  to  the  number  of  independent  data 
points  being  fit. 

The  primary  function  of  these  plots  is  to  provide  a  means  of 
comparison  of  the  results  of  different  applications  of  the  inversion 
to  a  given  data  set,  in  addition  to  providing  reassurance  that  a 
"reasonable  looking"  fit  is  being  obtained.  I  refrain  from  using 
them  to  compare  the  results  obtained  here  with  those  of  other 
researchers  since  there  is  no  common  ground  on  which  to  base 
such  a  comparison  Stein  and  Geller  (l978b)  make  use  of  plots  of 
smoothed  time  series  envelopes  to  display  the  results  of  their 
efforts  at  estimating  Q  However,  since  their  models  are  more 
heavily  constrained  than  mine  (being  constrained  by  the  assumed 
source  properties),  it  is  difficult  to  use  such  qualitive  comparisons 
as  the  plots  to  evaluate  the  relative  merit  of  the  two  techniques 
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A  somewhat  more  quantitative  evaluation  of  the  inversion 
results  can  be  made  by  appealing  to  normal  mode  excitation  theory 
to  supply  a  relationship  among  the  excitation  amplitudes  within  a 
given  multiplet  Stein  and  Geller  (1977)  derive  such  a  relationship. 
If  the  displacement  associated  with  an  earth  normal  mode  of  angu¬ 
lar  order  l  and  azimuthal  order  m,  Ui.m(r),  is  expressed  in  terms  of 
a  displacement  amplitude,  and  the  eigenfrequency.  uin,  vis  a 
vis  the  relation 

Ui.m(r.i)  =  +  (B,.m(r))*e 

then  the  displacement  amplitudes  within  a  multiplet  obey  the  fol¬ 
lowing  symmetry  relation: 

®1  .m  =  (  ,-m  ) 

where  •  denotes  complex  conjugation.  While  Stein  and  Geller 
derived  this  relation  for  a  point  double-couple  source,  by  using  the 
formalism  developed  by  Gilbert  (1971)  and  the  expressions  for  the 
normal  mode  strain  tensor  elements  given  by  Gilbert  and  Dziewon- 
ski  (1975),  it  is  possible  to  show  that  the  relation  holds  true  for  any 
point  source  Now,  for  gravimeter  records,  an  excitation  amplitude, 
Aim.  determined  by  the  inversion  routine  is  the  component  in  the 
recording  direction  of  the  seismograph  of  BhmOV}),  where  is  the 
reciever  location,  so  that  the  excitation  amplitudes  must  also 
satisfy  the  above  symmetry  relation.  For  strain  records  such  as  the 
Chilean  data,  the  Aim’s  are  a  component  of  VBk.m  *  H.mV.  Since  v  is  a 
purely  real  linear  operator,  V(Bi  m)*  +  (H ,m)*V  =  (7Bi  m  +  B(  so 


that  H.m  =  (ft.-**)*  implies  i4lm  =  (At  _m)*.  This  symmetry  leads  to  the 
following  constraints  on  the  amplitude  and  phase  components  of 
the  excitation  amplitudes: 

I  4j ,m  !  =  7t.m  =  ~~/l  ,-m 

where  |  I  and  y i,m  are  the  amplitude  and  phase  components  of 
Ai,m  defined  by  Alm  =  \At  n  These  symmetry  relations  provide 

a  means  of  evaluating  the  results  of  the  inversion  of  the  actual 
data  The  components  of  the  amplitudes  obtained  from  the  inver¬ 
sion  should,  within  their  uncertainty  limits,  satisfy  the  above  sym¬ 
metry  constraints.  By  refering  to  Table  2  and  to  the  plots  in  Figs 
4.5  to  4  8  in  the  following  section,  one  can  see  that  this  is  indeed 
the  case  for  every  excitation  amplitude  obtained  from  the  inversion 
of  the  Chilean  data  Where  the  uncertainty  is  relatively  small,  as 
with  the  m  =  ±  1  singlets  of  and  the  m  =  ±2  singlets  of  0-?2.  the 
estimates  come  quite  close  to  actually  meeting  the  symmetry  condi¬ 
tions.*  This  result  contributes  significantly  to  the  confidence  that  I 
have  in  the  reliability  of  the  results. 
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IV.  Signal  Detection  and  Analysis  Using  Quasi-Harmonic 
Decomposition  (QHD)  Methods 
Introduction 

In  many  research  efforts  in  geophysics,  time  series  analysis  plays  an 
Important  role.  The  data  that  the  seismologist  often  deals  with  is  often  quite 
complex,  consisting  of  many  distinct  components  (e.g.,  seismic  phases  of  vari¬ 
ous  types,  noise,  etc.)  superimposed  on  one  another.  Such  data  typically  must 
be  filtered  (in  a  variety  of  ways)  in  order  to  facilitate  the  extraction  of  the 
discrete  components  that  are  required.  The  data  components  of  interest  ordi¬ 
narily  consist  of  isolated  arrivals,  or  phases,  that  carry  the  basic  units  of  infor¬ 
mation.  One  such  filtering  scheme,  designed  specifically  to  isolate  individual 
arriving  phases  in  a  complex  time  series,  is  Quasi- Harmonic  Decomposition 
(QHD)*.  QHD  filtering  of  a  time  series  yields  estimates  of  the  arrival  times  of 
discrete  seismic  phases,  together  with  a  spectral  estimate  of  each  detected 
phase.  The  arrival  time  and  spectral  estimate  associated  with  a  given  phase  can 
then  be  used  to  reconstruct  the  arrival  in  the  time  domain-  The  net  result  is  a 
decomposition  of  e  complex  time  series  into  a  series  of  discrete  arrivals  that 
can  then  be  used  in  e  variety  of  ways  to  infer  properties  of  the  source  of  the 
phases  and  the  path  they  followed. 

Some  aspects  of  the  QHD  approach,  in  the  context  of  seismic  signal 
analysis,  are  discussed  by  Alexander  (1963),  Arcbambeau  at  at.  (1965),  Archam- 
beau  and  Flinn  (1965),  Dziewonski,  Block  and  Landisman  (1969),  Dziewonski. 
Mills  and  Block  (1972),  Instron,  Marshall  and  Blarney  (1971),  Cara  (1973)  and 
Denny  and. Chin  (1976),  among  others.  This  early  work  focused  on  narrow  band 

filtering  of  seismic  time  series  to  determine  the  dispersive  cheracteristics  of 

*Som«  of  the  terminology  u»ed  in  the  QHD  approach  ia  somewhat  specialised  and  therefore 
a  gloaaa rj  o!  the  moat  common  term*  and  aymhola  used  is  given  m  the  Append)*  l 
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surface  waves  and  on  polarization  filtering,  employing  the  multi-component 
information  contained  in  the  seismic  wave  field  to  extract  seismic  phases  hav¬ 
ing  particular  particle  motion  characteristics.  The  QHD  analysis  procedure  to 
be  described  here  involves  a  combination  of  dispersion  analysis,  threshold 
detection  and  polarization  analysis  (when  multi-component  data  is  processed), 
used  jointly  in  a  pattern  recognition  procedure  designed  to  automatically 
detect  and  isolate  body  wave  arrivals  of  a  particular  type.  (Other  versions  of 
the  computer  program  can  be  employed  to  isolate  and  analyze  surface  waves  of 
various  types.)  The  methods  and  analysis  procedures  are  described  in  Archam- 
beau  et  al.  (1982). 

Thus  the  QKD  analysis  procedure  used  in  the  computer  processing  is  capa¬ 
ble  of  identifying  body  waves,  of  prescribed  type,  in  a  background  of  noise  and 
other  signal  types,  and  provides  the  user  with  both  the  spectra  and  time 
domain  wave  forms  for  the  individual  signals  occurring  within  the  time  series. 
In  view  of  this  capability,  the  program  should  be  a  useful  tool  for  seismologists. 

Basic  Features  and  Methods  in  QHD  Analysis 

The  essential  feature  of  a  quasi-harmonic  decomposition  of  a  time  series, 
as  opposed  to  a  Fourier  decomposition,  is  that  6ofA  time  and  frequency  infor¬ 
mation  are  retained.  This  is  accomplished  by  parallel  filtering  of  the  time  series 
with  a  set  of  Gaussian  narrow  band  filters,  the  output  of  each  filter  correspond¬ 
ing  to  a  single  "quasi-harmonic  component"  of  the  original  time  series,  with 
each  of  these  having  a  modulated  sinusoidal  form.  It  is  the  modulation  of  the 
sinusoidal  output  that  contains  information  concerning  the  time  of  energy 
arrival  (group  times)  within  the  original  time  series,  while  the  sinusoid  (or  "car¬ 
rier”)  contains  the  spectral  information  related  to  these  arrivals.  (By  contrast, 
a  "Fourier  decomposition"  is  equivalent  to  filtering  with  a  set  of  delta  function 
filters,  centered  at  a  discrete  set  of  frequencies,  the  output  of  each  being  a 
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pure  sinusoid.  Since  the  sinusoids  are  not  modulated  they  contain  no  arrival 
time  information  in  themselves.) 

Thus  QHD  analysis  is  designed  to  provide  time  varying  spectral  information. 
While  moving  window  Fourier  analysis  also  provides  time  varying  spectral  esti¬ 
mates.  such  an  approach  suffers  from  poor  time  resolution  and  contamination 
due  to  truncation  effects,  which  can  be  severe.  On  the  other  hand,  QHD  pro¬ 
vides  time  varying  spectral  estimates  for  which  time  resolution  can  be  con¬ 
trolled  (automatically  by  a  computer  or  interactively  by  a  user)  and  truncation 
effects  do  not  occur.  However,  uncertainty  in  the  frequency  of  the  spectral 
estimate  obtained  does  occur,  and  is  the  price  to  be  paid  for  the  time  resolution 
achieved.  The  trade-off  between  time  and  frequency  resolution  in  QHD  is 
expressed  by  an  "uncertainty  principle"  having  the  form  AuAt  >  1/  2.  where 
and  &t  are  the  uncertainties  in  the  frequency  and  time  to  be  associated  with 
each  group  arrival  time  and  the  spectral  estimate  of  the  signals.  (For  Gaussian 
filters  Aodf  a  1/2.  while  for  non-Gaussian  filters  bu&t  >  1/2).  The  uncertainty 
in  frequency  resolution  results,  in  practice,  in  spectral  smoothing,  which  can 
usually  be  tolerated  and  is  often  even  desirable.  Since  Lu  is  nearly  equal  to  the 
bandwidth  of  the  Gaussian  filter  used,  then  this  uncertainty  principle  provides 
us  with  a  useful,  quantitative,  way  of  estimating  the  time  and  spectral  value 
uncertainties  associated  with  a  QHD  analysis. 

In  order  to  obtain  accurate  time  and  spectral  information  in  practice  (i.e. 
by  automated,  computer  controlled  analysis)  a  pattern  recognition  method  is 
used  to  detect  and  isolate  signals  of  specified  character.  The  pattern  recogni¬ 
tion  procedure  is  applied  in  the  signal  parameter  space  generated  from  the 
narrow  band  filtering  operations.  The  idea  is  to  have  the  computer  search  for 
regular  patterns  within  the  parameter  "space"  and  to  isolate,  or  pick  out.  those 
patterns  that  correspond  to  signals  of  interest  to  the  user.  For  example,  the 
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program  can  be  instructed  to  search  throughout  the  time  series  for  signals  with 
particular  dispersion  characteristics. 


Figure  1  -  Envelope  functions,  for  the  (reduced)  travel  time  interval  B  to  40  sec, 
plotted  on  the  time-frequency  plane.  The  original  seismic  time  series  is 
shown  in  the  panel.  Body  wave  signals  from  an  earthquake  in  a  high  noise 
background  arrive  at  (approximately)  15  sec,  and  can  be  detected  in  both 
the  original  seismogram  and  in  the  time-frequency  plane.  The  envelope 
amplitude  "ridges"  extending  to  high  frequency  are  discrete  body  wave 
arrivals  from  the  seismic  event  while  the  large  envelope  maxima  at  the 
lower  frequencies  are  (primarily)  due  to  the  noise  background. 

The  signal  parameters  are  generated  automatically  by  QKD  and  are 
derived  from  a  set  of  functions  associated  with  the  time  series  being  analyzed. 
In  particular,  at  each  narrow  band  filter  frequency  the  envelope,  or  modula¬ 
tion  function  of  each  filter  output  is  generated  along  with  the  instantaneous 
phase  and  frequency  associated  with  the  time  series.  Maxima  in  the  envelope 
functions,  corresponding  to  energy  group  arrival  times  at  the  frequencies 
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sampled  by  the  filters,  are  determined  by  the  program  and  this  data  set  is 
plotted  in  the  time-frequency  plane;  where  coherent  patterns  formed  by  the 
measured  group  arrival  time  data  define  the  dispersion  relations  for  "signals'' 
arriving  within  the  time  series  being  analyzed.  Figure  1  illustrates  the  time- 
frequency  plane  for  a  recorded  seismic  event  in  a  high  noise  background.  One 
function  of  the  pattern  recognition  component  of  the  QHD  program  is  to  locate 
and  "follow'’  coherent  "ridges"  of  envelope  maxima  in  the  t*f  plane,  and  to 
extract  (or  "isolate")  those  that  conform  to  "signals"  having  a  particular, 
desired,  dispersion  characteristic. 

In  addition  to  measuring  the  times  of  envelope  maxima,  the  program 
determines  both  the  instantaneous  phase  and  amplitude  of  the  envelope  at  the 
maxim.  (If  the  time  series  corresponds  to  a  spatial  component  of  the  wave 
field,  and  other  components  are  available  and  are  similarly  processed  by  QHD. 
then  the  polarization  is  also  computed  at  each  group  arrival  time.  The  polari¬ 
zation.  is  defined  to  be  the  phase  difference  between  the  vertical  and  horizon¬ 
tal  component  time  series,  at  a  point  of  observation.)  The  envelope  amplitude 
and  instantaneous  phase  at  the  group  arrival  times  are  directly  related  to  the 
Fourier  spectral  amplitude  and  phase.  Thus,  at  each  observed  group  arrival 
time-frequency  point  in  the  time-frequency  plane  a  signal  parameter  "vector" 
composed  of  the  estimated  Fourier  amplitude  and  phase,  and  polarization 
when  available,  is  defined  and  "attached"  to  each  measured  group  arrival  time 
value.  Thus,  in  addition  to  searching  for  coherent  patterns  of  arrival  time 
versus  frequency,  corresponding  to  those  of  interest  to  the  user,  the  program 
simultaneously  compares  the  spectral  amplitude  relative  to  noise  (end  the 
polarization  of  the  wave  field  when  available)  to  values  required  for  the  signal 
specified  by  the  user.  Thus,  the  program  will  perform  joint  dispersion  "fi'ter- 
ing".  threshold  detection  and  polarization  filtering.  The  actual  procedure  fol- 
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iowed  by  the  program  in  performing  this  "filtering-detection"  task  involves 
iterative  pattern  recognition  in  the  time-frequency  plane,  wherein  signal 
detections  are  "removed"  (spectrally  subtracted)  from  the  time  series  and  the 
whole  filtering-pattern  recognition  procedure  repeated  using  the  residual  time 
aeries.  This  iterative  procedure  is  designed  to  combat  interference  effects 
among  closely  time  spaced  signals  and/or  noise  "pulses". 

Once  a  set  of  candidate  detections  of  signals  have  been  made  on  the  basis 
of  the  criteria  indicated,  they  are  considered  as  a  group  by  the  program,  using 
criteria  for  association  of  the  isolated  "pulses"  detected,  into  seismic  phases. 
In  particular,  if  QKD  obtains  two  or  more  signal  detections  within  the  time 
resolution  range  of  the  narrow  band  filters  employed  in  the  analysis,  then 
these  detections  will  be  merged  to  form  a  single  seismic  phase.  (Such 
occurrences  are  fairly  common  and  are  due.  in  part,  to  the  multiple  pattern 
searches  using  residual  time  series.  In  this  case  a  signal  can  be  detected  more 
than  once  if  it  was  not  completely  removed  in  the  first  pass.)  However,  the 
instantaneous  frequency  variation  with  time  for  the  original  time  series  is  used 
to  check  for  indications  of  interference  within  the  filter  "resolution  time  win¬ 
dow",  which  may  occur  when  two  or  more  seismic  phases  are  present.  When 
significant  variation  of  the  instantaneous  frequency  from  the  filter  center  fre¬ 
quency  occurs  at  a  number  of  filter  frequencies  within  the  resolution  time  win¬ 
dow.  then  the  detections  within  the  window  are  not  merged,  but  are  considered 
to  be  separate  seismic  phases. 

Once  the  seismic  phases  have  been  defined  in  terms  of  a  set  of  signal  vec¬ 
tors  (containing  the  spectral  amplitude  and  phase)  associated  with  a  set  of 
group  arrival  times,  the  program  inverts  the  signal  spectral  information  to  the 
time  domain  and  thereby  generates  waveforms  for  each  signal  identified  The 
program  also  computes  the  envelope  functions  for  these  wave  forms  and 
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the  first  inflection  point  in  the  envelope  to  estimate  the  first  arrival  time  (sig¬ 
nal  arrival  time)  for  the  seismic  phase.  The  program  then  provides  a  final  list 
of  seismic  phase  detections.  The  associated  spectral  variables,  as  well  as  sig¬ 
nal  timing  data  for  each  detection,  are  written  to  an  output  file.  A  separate 
utility  program  will  then  plot  results  for  the  user. 

Program  Operations  and  Structure 

A  summary  of  the  program  operations,  in  the  order  in  which  they  are  per¬ 
formed  is  as  follows. 

(1)  Read  an  input  data  file  containing  the  time  series  to  be  processed. 

(2)  Preprocess  the  time  series:  Bandpass  filter,  detrend,  taper,  remove 
spikes,  etc. 

(3)  Obtain  a  noise  sample  from  the  data  record  and  compute,  or  select 
from  an  input  file  supplied  by  the  user,  appropriate  narrow  band 
filter  parameters. 

(4)  Decompose  the  time  series  into  quasi-harmonic  components  by 
Gaussian  narrow  band  filtering. 

(5)  Compute  associated  functions,  including  the  envelope,  instantane¬ 
ous  phase  and  instantaneous  frequency  functions,  for  each  filter. 

(6)  Determine  the  times  of  envelope  maxima  and  the  amplitudes  of  the 
envelope  functions  at  these  maxima  for  each  filter.  Also  determine 
the  instantaneous  phase  values  at  these  times.  (Also  compute  the 
polarization  and  emergence  angle  of  the  wave  field  at  the  group 
arrival  times,  when  multi-component  data  is  processed.) 

(7)  Search  for  and  identify  specified  types  of  signals,  based  on  pattern 
recognition  in  the  time-frequency  plane,  having  appropriate  disper- 
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iion  and  signal  power  levels  (relative  to  noise).  (Also  select  only 
those  with  proper  polarization  and  emergence  angles,  when 
multiple-component  data  is  being  processed.) 

(B)  Subtract  tbe  spectral  values  (amplitude  and  phase)  associated  with 
detected  signals  from  the  time  series  spectrum  and  invert  the  "resi¬ 
dual"  spectrum  to  obtain  b  "residual"  time  series".  Repeat  the  pro¬ 
cessing  steps  (4)  through  (B)  on  residual  time  series.  (Iterative  pat¬ 
tern  recognition.)  Terminate  when  step  (7)  produces  no  new  detec¬ 
tions  or  when  a  preselected  number  of  such  iterations  have  been 
performed. 

(9)  Associate  the  collection  of  signal  detections  produced  after  exit 
from  step  (8)  into  a  set  of  seismic  phases,  using  criteria  based  on 
the  resolution  (or  uncertainty)  time  of  the  Gaussian  filters  employed 
and  the  instantaneous  frequency.  (Termed  the  "culling  process'  ) 

(10)  Invert  the  isolated  seismic  phase  spectra  to  the  time  domain  to  give 
wave  forms  and  determinations  of  the  "first  motion"  times  of  each 

(11)  List  results  and  write  output  files. 

(12)  Return  to  read  a  new  data  file  or  terminate  run. 

The  actual  structural  elements  and  flow  of  the  program  operations  are 
indicated  in  Figure  2.  Detailed  operational  steps  within  the  input/output  sec¬ 
tions  and  the  various  execution  sections  indicated  in  the  diagram  are 
described  separately  in  the  Users  Guide  to  QHD  and  in  the  QHD  Program 
Operations  Manual. 

Examples  of  QHD  Analysis  Results 

Applications  of  the  QHD  analysis  approach  to  both  observed  and  syntheti¬ 
cally  produced  seismic  time  series  provide  the  best  means  of  evaluating  the 
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entire  "filtering- pattern  recognition"  approach  used  in  this  program.  With 
applications  to  synthetic  "data"  we  can  systematically  (and  quantitatively) 
determine  the  accuracy  and  sensitivity  of  the  method  as  a  function  of  noise 
and  signal  properties  using  entirely  known  signal  input  which  can  closely 
approximate  real  seismic  data.  On  the  other  band,  since  observed  seismic  data 
eon  differ  from  synthetic  data  in  many  (unknown)  ways,  it  is  also  important  to 
test  and  evaluate  the  method  using  a  variety  of  observed  data.  Usually  this 
involves  interpreting  the  results  of  the  processing  in  the  context  of  what  is 
known  about  the  source  and  path  of  propagation.  For  example,  large  explo¬ 
sions  observed  at  regional  and  teleseismic  distances  can  produce  relatively 
simple  seismograms,  with  signal  arrival  times  that  are  predictable  with  fair 
accuracy.  In  these  cases  the  QKD  results  can  be  compared  with  expected 
arrival  times,  Wi.ve  forms  and  spectra  of  the  most  prominent  seismic  phases 

a 

from  these  more  "controlled  events  and  "known”  paths,  to  evaluate  perfor¬ 
mance.  Similar  statements  apply  to  reflection  seismograms  from  explosions  or 
vibrator  sources  recoded  in  "test  areas"  where  the  seismic  velocity  structure 
is  reasonably  well  known. 

The  following  figures  provide  some  examples  of  QKD  processing  for  both 
synthetic  and  observed  data.  These  examples  constitute  only  a  very  small  and 
incomplete  part  of  the  systematic  testing  which  is  required  for  a  full  evalua¬ 
tion.  These  examples  are  actually  only  meant  to  demonstrate  the  kind  of  tests 
performed  and  the  general  nature  of  the  results  obtained  for  a  range  of 
"data";  from  simple  pulse  synthetics  to  very  complex  observed  earthquake 
data.  They  also  graphically  illustrate  the  steps  in  the  processing  and  the 
nature  of  the  output  obtained. 

Figure  3  illustrates  results  of  QHD  analysis  of  a  sequence  or  fairly  closely 
spaced  seismic  pulse  synthetics,  in  a  background  of  random  noise.  As  Figurr 
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Time  domain  synthesis  of  the  CKD  estimated 
undlaperaed  signal  spectra  (traces  labe.ed 
1  through  5).  Top  trace  is  the  origins, 
sequence  of  five  synthetic  pulses,  the 
second  trace  la  the  sum  of  the  five  Q?C 
isolated  pulses. 


Figure  3a.  A  staple  test  example  of  QHD  signal  Isolation  using  single  component  synthetic 
data.  The  signal  pulses  have  maximum  power  near  40  Hr.  and  approximate  pulse  forms 
generated  in  exploration  seismic  surveys.  The  QHD  pattern  recognition  in  the  time- 
frequency  plane  uses  threshold  detection  combined  with  "dispersion  filtering"  to 
automatically  isolate  the  five  undispersed  signal  pulses. 


i 

I 


-122- 


Figure  3b.  Spectral  estimates  of  each  of  the  five  Isolated  test  pulses.  The 

QMD  spectral  amplitude  estimate  Is  very  close  to  the  knovn  Fourier  amplitude, 
while  the  phase  is  sometimes  less  accurately  determined  but,  nevertheless, 
with  sufficient  accuracy  to  yield  time  domain  reconstructions  that  are  very 
close  to  the  original  time  domain  pulses. 
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3a  shows.  in  a  low  noise  environment  such  as  this,  the  isolation  of  multiple 
non-overlapping  undispersed  pulses  is  accomplished  with  high  accuracy  by 
QHD  pattern  recognition  in  the  time- frequency  plane,  and  the  automatic  tim¬ 
ing  of  the  pulse  arrivals  is  quite  accurate.  Figure  3b  shows  the  spectral  esti¬ 
mation  of  each  of  the  five  pulses,  obtained  from  the  amplitudes  of  the  envelope 
functions  at  their  maxima  in  tbe  t-f  plane  and  the  instantaneous  phase  values 
at  the  times  of  these  maxima.  As  noted  in  the  figure  caption,  the  spectral 
phase  is  not  as  accurately  determined  as  the  spectral  amplitude.  In  the  fre¬ 
quency  range  around  the  spectral  maximum,  however,  the  phase  values  are 
quite  accurate  and  hence  the  time  domain  wave  forms  generated  from  these 
spectral  estimates  closely  match  the  original  pulses. 

Figure  4  illustrates  QHD  isolation  of  two  overlapping  pulses  in  a  high  noise 
background.  In  this  example  the  noise  is  non-random,  having  the  dispersive 
structure  of  a  seismic  Rayleigh  wave.  Figure  4a  shows  the  QHD  determined 
dispersion  of  the  synthetic  Rayleigh  wave.  Figure  4b  shows  tbe  isolation  of  the 
two  pulses,  by  the  QHD  "dispersion  filtering".  In  this  moderately  demanding 
test  case  the  arrival  times  of  the  pulses  are  very  well  determined,  while  the 
spectral  estimates,  end  the  wave  forms  generated  from  them,  are  somewhat 
distorted  due  to  interference  between  the  two  pulses  themselves,  as  well  as 
from  the  Rayleigh  wave.  (No  analytical  corrections  were  made  for  interference 
in  this  example.) 

Figure  S  shows  a  similar  example,  with  five  variable  amplitude  pulses  dis¬ 
tributed  throughout  the  Rayleigh  wave  signal  train.  As  can  be  seen  from  the 
patterns  in  tbe  time-frequency  plane,  the  Rayleigh  wave  strongly  distorts  the 
pulse  dispersion.  However,  the  pattern  recognition  procedure  employed 
nevertheless  accurately  identifies  and  times  all  the  pulses.  It  can  also  be  seen 
that  the  spectral  estimates  and  associated  wave  forms,  while  resembling  the 
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Figure  4a.  Example  of  the  Isolation  of  two  closely  spaced  body  wave  signal  pulses 
in  a  Rayleigh  type  surface  wave  background.  The  pure  (unmixed)  signals  are 
shown,  along  with  the  QHD  generated  tiae -frequency  plane  results  for  the  pure 

Rayleigh  wave. 
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original*  reasonably  well,  are  somewhat  distorted  by  the  Rayleigh  wave  "noise". 

Figure  6  provides  a  reasonably  complete  example  of  the  application  of  the 
QHD  automated  analysis  to  a  complex  teleseismic  earthquake.  In  particular, 
this  example  illustrates  (1)  the  automated  selection  of  Gaussian  narrow  band 
filters  (which  is  based  on  the  event  signal  to  noise  ratio  as  a  function  of  fre¬ 
quency)  in  Figure  6a;  and  (2)  the  iterative  pattern  recognition  in  the  time- 
frequency  plane,  in  Figure  6b.  All  of  these  results,  in  Figures  6a  through  6d 
are  part  of  the  standard  output  from  the  program  and  serve  to  provide  a  user 
with  a  good  mean  of  evaluating  the  final  results,  which  are  the  spectral  and 
timing  information  in  Figure  6c  and/or  the  isolated  wave  form  results  in  Figure 
6d. 

The  final  example,  in  Figure  7.  is  an  application  or  this  analysis  method  to 
a  relatively  "high"  frequency  (0-20  Hz)  regionally  recorded  earthquake.  In  this 
example,  as  can  be  seen  simply  from  the  original  time  series  in  Figure  7a.  the 
noise  level  is  higher  than  in  the  previous  example.  The  "event"  is  also  very 
complex,  consisting  of  many  closely  spaced  seismic  phases  resulting  from  vari¬ 
able  energy  radiation  from  the  spatially  finite  source  itself,  plus  multiple 
arrivals  due  to  many  strong  reflectors  and  refractors  at  depth  between  the 
source  and  receiver.  Here  again,  many  of  the  important  QHD  program  opera¬ 
tions  enumerated  in  the  previous  section  are  illustrated,  in  particular  the  pat¬ 
tern  recognition  operations  which  are  shown  in  Figures  7b  and  7c. 

The  most  important  observation  to  be  made  regarding  this  method  is  that 
the  QHD  signal  isolation  provides  individual  spectral  information  for  each 
coherent  signal  within  the  complicated  event  wave  train  and  thereby  allows 
rapid,  quantitative  analysis  for  path  and  source  properties  to  be  undertaken. 
The  confidence  which  one  may  place  in  the  signal  identifications,  timing  and 
spectral  estimates  is  directly  proportional  to  the  accuracy  of  the  method  when 


Figure  6a 

(•) 


Shore  Period  ciUhuiIc  alfn*l»  fro* 
an  aarthquaka.  The  compression*)  wave 
crala  la  iIimi. 


r 


llfitl  wa.  Itolaa  l^adrt 
6a>Uti.f«  «#*  •  > 


Amplitude  apactra  froa  a  long  algnal 
window,  consisting  of  the  entire  P  wave 
train,  and  a  long  nolaa  aamrl*  The 
frequency  hand  within  which  OK?  spectral 
analvaia  la  to  taka  place  la  indicated, 
and  la  automatically  aat  by  the  program. 


t. 


*> 


Distribution  of  Gauaslan  filters  that 
have  been  automatically  dcsipnej  and 
distributed  within  the  ftequenev  bard 
aalectad  by  the  program. 


.  Example  of  the  sequence  of  QKD  operetlons  epplied  to  teleseismic  data. 
Design  of  the  Ceusslen  filters. 


•.mm** 

-a.m"** 


•.ti- 

«.«j  -i 

•.ti.l 

:•;:  *•  ‘  • 
A  *1  _  • 


QND  gsnarat ad  ciaa-f raquancy  plan*  data, 
with  laolatad  "arrival*',  corrasponomg  to 
undlaparaad  algnal  pulsts.  Tha  original 
tiaa  MfiN  (ahotm  at  the  top;  i*  procaaaad 
la  thla  "first  Ittrstion",  with  tha  apoctral 
•itlaatta  ohtalaad  for  body  wave  arrivals 
to  ho  subtract  ad  fro*  thl*  original  uat 
aarloa  to  "ran  dual"  tin  aariat. 

which  will  again  ha  aaarchtd  by  tha  progra* 
for  uadlaparaad  signal*. 


S.XT.f- *♦*  J 

s.a*Mt*aa 


i.»-  # 

!.«• 
It- 
1«- 
l.U- 
!-•*- 

•*-  a 

•-•i-  # 

•.71 - 

•-•>-  • 

Mi  -  % 

•  •  : 

Mi-  # 
•-*1 


mt  *•»■**% 

a 

1 \!\J> 

Ijl 

i 

t 

t 

'if"  l 

r,  ~ 

jl/i-' 

p 

l*  a  hi 

r 

i 

i 

-»  • 

QKD  ganaratai  titaa-fracutncv  plana  data  for 
tha  aacond  Itaratirr  rs»iiwal~tica  aana* 
fro*  which  pravlously  octacted  "aignala". 
laolatad  in  tha  first  dataction  natation, 
have  bean  (opact rally;  raorvad 


hA>  1 


»  •  »  1  » 


QKT>  immtd  r-cu-nev  »l«n.  d«t»  f'T 

t h.  third  it.r.tl.n  tJn.  im«  <*•'*  vP>  1  c ^ 
rr.viou.ly  dnt.ct.d  ".if-.l.  ,  i».-.«i»d 
In  ih»  llr.r  .nd  Mcone  «.«cti.'n  is.r.tlon. 

h*v«  h«tn  rnoiii. 


Flgura  6b.  QHD  operations  applied  to  teleseismic  dote,  (b)  Iterative  pattern 
recognition  in  the  tine-frequency  plane. 
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QKD  generated  tl«e-freue#ncy  plane  data, 
with  isolated  P  wave  "arrivals”,  corresponding 
to  uBdliperaad  signal  pulse*.  ?ha  on|iml 
t im  series  (shown  at  ch«  top)  la  processed 
la  tha  "first  iteration",  with  the  spectral 
eat lasts*  obtained  for  tha  P  wav*  arrivals 
to  ha  subtracted  fro*  this  original  tine 
•arias  to  generate  "residual"  tioa  aariaa. 
which  will  stain  ba  searched  by  tha  progras 
for  uedlapereed  signals. 
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Figure  7b.  QHD  operations  applied  to  regional  earthquake  date,  (b)  Iterative 
pattern  recognition  in  the  tiae-frequency  plane  (P  wave  aection). 
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Figure  7c.  QHD  operation*  applied  to  regional  earthquake  data,  (c)  Iterative 
pattern  raco*nltlon  in  the  tlaa-fraquancy  plane  (S  wave  section). 
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Figure  7d.  QHD  operations  applied  to  regional  earthquake  data,  (d)  Examples 
of  QHD  spectral  estimates  for  isolated  body  wave  signals  (F  and  S  waves). 
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Figure  7e.  QKD  operations  applied  to  regional  earthquake  data,  (e)  Tlae  domain 
signal  pulse  estimates. 
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applied  systematically  to  synthetic  "data",  such  as  that  illustrated  in  Figures  3, 
4  and  5,  and  similarly  when  applied  to  observed  data  from  tests  under  con¬ 
trolled  conditions.  In  addition  however,  the  closeness  of  the  reproduction  of 
the  original  time  series  by  the  summation  of  the  coherent  signals  obtained  by 
the  QHD  analysis  is  clearly  indicative  of  the  accuracy  and  validity  of  the 
results,  even  though  the  signal  content  of  the  original  time  series  may  be  unc¬ 
ertain  or  unknown.  As  illustrated  in  Figure  7e  (and  in  Figure  6d)  the  sum  of 
the  isolated  signals  does  closely  approximate  the  original  event  time  series,  so 
that  a  reasonable  degree  of  confidence  can  be  attached  to  the  detailed  results. 
With  additional  systematic  testing  and  extensive  applications  in  the  future  it 
should  be  possible  to  rigorously  quantify  this  degree  of  confidence. 
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Appendix  1 


GLOSSARY  OP  SYHBOLS,  ABBREVIATIONS  AND  TERMINOLOGY 


Tranter  Functions:  Spectrum  (amplitude  and  phase)  of  the  impulse  response  of 
a  niter.  Used  to  define  response  of  both  sensor  system  and  wave  propaga¬ 
tion  in  the  earth.  (Each  wave  type  has  a  different  transfer  function.) 

Quasi- Harmonic  Signal  Decomposition  (QHD):  Literally,  the  separation  or 
decomposition  of  a  time  series  into  a  set  of  nearly  harmonic  (or  "quasi- 
harmonic")  components,  from  which  the  envelope,  instantaneous  phase, 
frequency  and  polarization  time  series  are  generated.  Functionally, 
narrow-band  filters  are  used  to  generate  nearly  harmonic  output  at 
discrete  frequencies  (Le..  in  narrow  frequency  bands),  from  which  associ¬ 
ated  filter  envelope  functions;  instantaneous  phase,  frequency  and  polari¬ 
zation,  can  be  computed  as  functions  of  time  at  the  discrete  set  of  filter 
frequencies. 

Hatched  Filtering:  Cross  correlation  of  an  expected  or  sought  for  signal  with  a 
recorded  time  series,  with  a  "detection"  corresponding  to  a  pulse  (peak) 
In  the  correlation  time  series.  In  general:  the  operation  of  correlation  of 
expected  or  sought  for  signal  properties  with  "observed"  properties  of  a 
time  series. 

Pattern  Recognition:  Search  procedure,  usually  in  two  or  more  dimensions 
within  a  parameter  space,  aimed  at  the  selection  or  detection  of  a  series, 
or  pattern,  of  parameter  values  that  are  characteristic  of  a  signal.  The 
matched  filter  concept  and  methods  are  often  used  for  selection  or  detec¬ 
tion  of  sought  for  signal  within  the  parameter  space. 

DiaperMon:  General  term  referring  to  the  variation  of  group  velocity  (or  group 
arrival  time)  as  a  function  of  frequency  or  the  variation  of  phase  velocity, 
or  phase  arrival  time,  as  a  function  of  frequency. 

*NBF  -  Narrow- Band  Filter  Narrow-band  filter,  usually  of  Gaussian  type  with 
transfer  function: 

=  £7  V*  «P  ^  ^ 

where  fi  -  tn2/2.  Here  Qfl~Auj/u},  is  the  filter  quality  factor,  with 
A Uj  *  2nA/ j  equal  to  one  half  the  frequency  bandwidth  at  the  half  power 
points,  and  vj  representing  the  center  frequency  of  the  jth  filter  (j  =  1,2 
...)  from  a  set  of  such  filters. 


Gaussian  filter  Resolution  Time,  :  The  time  interval  between  two  signal  pulses 
that  is  generally  required  for  the  complete  time  resolution  of  the  signals 
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in  the  narrow  band  filtered  output.  For  a  Gaussian  filter,  of  half 
bandwidth  A»  at  the  filter  half  power  points,  then  btbu  c*  1/2. 

MHF- Signal.  s(t Jj),  and  NBF- Quadrature  Signal  *(*;/*)••  The  NBF  signal 
corresponds  to  the  direct  output  of  the  filter  (inverse  Fourier  transform 
of  the  product  of  the  signal  spectrum  and  the  filter  transfer  function  or 
the  time  convolution  of  the  signal  time  series  and  the  impulse  response  of 
the  filter)  while  the  quadrature  NBF-signal  corresponds  to  the  Hilbert 
transform  of  the  filter  output,  that  is  the  output  signal  shifted  in  phase  by 
sr/2.  The  quadrature  signal  may  be  obtained  by  multiplying  the  signal 
spectrum  by  exp  [— t  sgn  (o)tt/  2],  followed  by  multiplication  of  this  by  the 
filter  transfer  function  and  inverse  Fourier  transforming  the  result. 

Analytic  (or  complex)  NBF-signal.  *«(*;/>):  Complex  signal  representation 
formed  Jrom  the  NBF  signal  and  quadrature  signal; 

*«(*;//)  +  **(*'•/>)•  Viewed  in  the  complex  signal  plane  (looking 

down  the  time  axis)  a  vector  from  the  origin  to  a  point  on  yc  describes  the 
"instant  aneous  phase"  of  the  signal  by  the  angle  it  makes  with  the  real  yc 
aids.  The  rate  of  rotation  of  the  vector  about  the  origin,  with  time, 
describes  the  time  rate  of  phase  change  or  the  "instantaneous  fre¬ 
quency." 

binlope  or  Modulation  Function,  £(t  ;fj):  Filter  envelope  (or  modulation)  func¬ 
tion.  equal  to  the  values  of  the  extreme  variations  of  the  narrow-band  fil¬ 
tered  time  series.  The  envelope,  or  modulation  function,  is  formed  as  a 
time  series  from  the  filtered  signal  time  series,  s(t ;//).  and  its  quadrature 
time  series.  s(t;fj)  using:  E  -  [se  +  The  envelope  function  is  equal 

to  the  magnitude  of  the  analytic  signal. 

Group  Time.  Group  (or  energy)  arrival  time  corresponding  to  the  time  of 

a  peak  in  the  filter  envelope  function  E(t  :fj)  tor  a  "NBF’  filter  with  ctiiter 
frequency  fj.  The  group  velocity  £/(/)  is  the  velocity  at  which  energy  is 
propagated,  and  is  proportional  to  the  source  to  receiver  distance  divided 
by  t,. 


Group  Spectrum.  Ag(fj):  Envelope  amplitude  at  a  maximum  point,  occurring  at 
the  group  time  tt(fj).  This  amplitude  is  very  nearly  equal  to  the  Fourier 
spectral  amplitude  of  a  signal  at  the  frequency  / i  for  a  properly  designed 
(Gaussian)  filter. 

bstantanaoua  Phase.  p(t;/j):  Computed  at  each  time  point  as  the  inverse 
tangent  of  the  ratio  of  the  quadrature  NBF  signal  to  the  actual  NBF  signal. 
At  a  group  time,  tt.  it  is  very  nearly  equal  to  the  Fourier  phase  of  the  sig¬ 
nal. 

kiatantaneous  Frequency.  0(t;/<):  Corresponds  to  the  time  rate  of  change  of 
the  instantaneous  phase,  n  *  dp/dt .  Its  value,  for  a  NBF  output,  will  be 
close  to  the  center  of  frequency  of  the  particular  filter  output,  except  at 
times  when  two  signal  pulses  interfere,  at  which  time  it  may  have,  essen¬ 
tially.  any  value. 


ln»UnUnw>M  Polarisation,  P[t;fj):  Defined  as  the  instantaneous  phase  differ¬ 
ence  between  two  spatial  components  of  the  vector  displacement  (or  velo¬ 
city  or  acceleration)  field;  for  example:  Pgn{t \Jj)  =  t:fj)  -  pit (f ;/y)  for 
the  vertical  (Z)  and  radial  (R)  vector  components.  Only  the  value  at  the 
group  arrival  time  for  a  signal  is  a  good  approximation  to  the  signal  polar¬ 
ization  or  particle  motion. 

Time  Varying  Apparent  Emergence  Angle,  e  (t  ;/j):  Apparent  wave  number  vector 
orientation  in  the  vertical  plane  to  the  earth's  surface,  at  a  receiver  point. 
The  emergence  angle  is  defined  as  the  angle  between  the  upward  vertical 
and  the  wave  number  vector,  where  |k.  |  =  2n/\.  with  A  »  wavelength.  The 
apparent  emergence  angle  is  computed  from  the  observed  horizontal  dis- 
placement  amplitude.  Ah,  and  the  vertical  displacement,  Ag.  as: 
e  =  tan"1  [Ag/ Ay\  The  true  emergence  angle,  e.  is  related  to  e  and  can 
be  computed  from  it  when  the  elastic  properties  of  the  surface  material  at 
the  receiver  point  are  known.  The  amplitudes  Ah  and  Az  ere  the  appropri¬ 
ate  group  spectral  values,  A,  (/^). 

Deterministic  "Noise"  Correction,  6Ah‘  Correction  applied  to  a  particular  signal 
group  spectral  amplitude  Ag(fj),  computed  from  other  nearby  group 
arrivals,  which  can  be  noise  or  other  (interfering)  signal  arrivals.  The 
correction  is  computed  analytically  using  the  first  order  theoretical 
response  of  a  narrow  band  Gaussian  filter  to  pulse  like  signal  input.  The 
corrected  group  spectral  value  is  A}  *  Ag  +  6As.  where  6A#  is  algebraic 
and  complex. 


Statistical  Noise  correction.  AAjt-  Noise  related  uncertainty  estimate  applied  to 
the  signal  amplitude  spectrum,  i.e..  A?  =  Ag  ±  LAs-  This  uncertainty  esti¬ 
mate,  at  a  particular  frequency,  is  taken  to  be  equal  to  the  mean  noise 
level.  As(f),  times  the  probability.  Pit,  of  occurrence  of  an  unresolved 
noise  pulse.  For  an  NBF  time  series  P*  is  taken  to  be  equal  to  the  number 
of  noise  pulses  occurring  per  unit  time,  times  twice  the  resolution  time, 

-  f  v 

At.  of  the  narrow  band  filter.  Thus,  AAjyr  =  As(f  )  Pn  ~  2Af  1 


T0 


As 


where  N  are  the  total  number  of  noise  arrivals  in  a  sample  time  segement 

Ta 


Signal  Arrival  Time  ("Time  of  First  Motion"),  t*:  The  time  of  onset  of  a  signal 
pulse.  Corresponds  to  the  travel  time  associated  with  a  "Signal  velocity" 
or,  equivalently,  with  the  ray  theoretic  arrival  time. 
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Appendix 

P  Wave  Amplitude  Variations  at  WWSSN 
Stations  in  the  U.  S.  and  Canada 
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Figure  Al.  Mean  relative  amplitudes  from  earthquakes  in  the  Kurile  Islands 
are  compared  with  all  data  to  a  northwest  azimuth  from  the 
United  States. 


RELATIVE  AMPLITUDE 
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Figure  A2.  Mean  relative  amplitudes  from  earthquakes  in  Japan  are  compared 
with  all  data  to  a  northwest  azimuth  from  the  United  States. 
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RELATIVE  AMPLITUDE 
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Figure  A3.  Mean  relative  amplitudes  from  earthquakes  in  the  Bonin  Islands 
are  compared  with  all  data  to  a  northwest  azimuth  from  the 
United  States. 
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Figure  Ai .  Mean  relative  amplitudes  from  other  northwest  earthquakes  no: 

located  in  the  Kuriles,  Japan,  or  Bonin  Islands  are  compared 
with  all  data  to  a  northwest  azimuth  from  the  United  States. 
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Table  A22 


Estimated  Sediment  Amplification  Corrections 
for  Wi'SSN  and  SDCS  Stations 


Station 

Amplitude 

Correction* 

Station 

Amplitude 

Correction* 

AAM 

1.4 

JCT 

mm 

ALQ 

•  * 

LON 

mm 

ATL 

•  • 

LUB 

1.8 

BEC 

1.8 

MDS 

•  * 

BKS 

1.3 

MSO 

•  • 

BLA 

0B2NV 

•  * 

B02 

•  * 

OGD 

•  * 

COR 

*• 

OXF 

1.8 

DAL 

1.6 

RCD 

1.6 

DUG 

•  • 

RKON 

•  » 

FLO 

1.5 

SCP 

»* 

GEO 

•  • 

SHA 

1.8 

GOL 

•  • 

TUC 

•  » 

CSC 

•  • 

IVES 

•  • 

KNME 

•  • 

’Station  amplitude  to  be  divided  by  correction. 

"No  sediment  correction  necessary,  in  the  case  of  a  limestone 
of  unknown  hardness,  no  correction  assumed. 


Estimated  Sediment  Amplification  Corrections 
for  Cleary  (196?)  LRSM  Stations 


Station 

Log  Amplitude 
Correction* 

Station 

Log  Amplitude 
Correction* 

ARWS 

*• 

KCAZ 

•  ft 

ATNV 

•  • 

KNUT 

(-0.20) 

AYSD 

(-0.26) 

LCNM 

•  ft 

BDPA 

(-0.11) 

MMTN 

(-0.15) 

BFCL 

•  * 

MNNV 

•  ft 

BLWV 

(-0.15) 

MPAR 

(-0.11) 

BUQB 

•  • 

MV  CL 

•  ft 

CKCB 

(-0.15) 

X1ZAR 

(-0.11) 

CP  CL 

•  • 

NDCL 

•ft 

CTOK 

(-0.26) 

NGWS 

•ft 

CVTN 

(-0.15) 

TCP  NT 

(-0.26) 

DENY 

(-0.11) 

PFM1 

•  ft 

DRCO 

•  • 

pmwy 

•  ft 

EYXV 

(-0  15) 

PT0R 

•  * 

FML’T 

•  • 

RK0N 

•  • 

FRMA 

(-0.26) 

RTNM 

(-0.20) 

FSAZ 

•  • 

RYND 

(-0.26) 

CDVA 

(-0.15) 

SEMN 

•  ft 

G1MA 

(-0.26) 

SJTX 

(-0.26) 

CON'B 

(-0.26) 

SSTX 

•  ft 

CVTX 

(-0.20) 

STNV 

(-0.04) 

HBOK 

(-0.15) 

TFCL 

(-0.20) 

HKND 

(-0.26) 

TKWA 

•  ft 

HKWY 

(-0.26) 

TUPA 

(-0.11) 

HL1D 

•  • 

WINY 

(-0.1B) 

HMHC 

•  • 

WNSD 

(-0.26) 

HNME 

•  • 

WTTX 

(-0.15) 

HTMN 

WY.'UT 

•  ft 

*Log  amplitude  correction  to  be  added  to  station  log  amplitude. 

**No  sediment  correction  necessary;  in  the  case  of  a  limestone  site  of  u 
hardness,  no  correction  assumed. 


Table  A24 

Inferred  Qa  for  LRSM  Stations  from  Evernden  and  Clark  (1970) 
Which  Were  Not  Covered  in  Cleary  (1967) 


Station 

Log 

Amplitude’ 

Estimated 

Sediment 

Correction 

Bedrock 

Log  Amplitude 

Qa 

(Inferred) 

(-0.15) 

(-0.21) 

89 

(-0.11) 

(-0.02) 

133 

ess 

(-0.15) 

(-0.20) 

(-0.35) 

72 

CUNY 

(-0.21) 

•  * 

(-0.22) 

B7 

DUOK 

(-0.32) 

(-0.20) 

(-0.52) 

59 

EDMT 

(0.15) 

•  * 

(0.15) 

243 

EKNV 

(-0.13) 

(-0.15) 

(-0.2S) 

79 

LSNH 

(0.032) 

•  • 

(0.032) 

154 

SKTX 

(0.C53) 

(-0.15) 

(-0.057) 

113 

’includes  (-0.02)  correction  from  least-squares  fit  to  common  data  in  Cle 
(1967)  plus  (-^-0.01)  correction  to  V.Ti'SSN-SDCS  baseline 
No  sediment  correction  necessary;  in  the  case  cf  a  limestone  site  of  unit 
hardness,  no  correction  assumed. 


Table  A25 


Inferred  <?0  for  LRSM  Stations  from  Booth  et  al.  (1974) 
Which  Were  Not  Covered  in  Cleary  (1967) 


Station 

Log 

Amplitude* 

Estimated 

Sediment 

Correction 

Bedrock 

Log  Amplitude 

(Inferred) 

AXAL 

(0.10) 

(-0.26) 

(-0.16) 

97 

BEFL 

(0.33) 

•  • 

(0.33) 

2500 

BRPA 

(0.02) 

(-0.11) 

(-0.09) 

112 

FKCO 

(0.19) 

(-0.26) 

(-0.07) 

117 

KYMA 

(-0.32) 

(-0.26) 

(-0.58) 

55 

JELA 

(0.29) 

(-0.26) 

(0.02) 

149 

JRAZ 

(-0.16) 

•  • 

(-0.16) 

97 

KCMO 

(0.12) 

•  * 

(0.12) 

212 

LYLA 

(0.19) 

(-0.26) 

(-0.07) 

117 

M01D 

(0.41) 

(-0.20) 

(0.21) 

337 

SWMA 

(0.04) 

(-0.15) 

(-on) 

107 

includes  (-0.01)  correction  from  least-squares  fit  to  common  data  in  Cleary 
(1967)  plus  (+0.01)  correction  to  WWS5N-SDCS  baseline,  therefore  no  correction. 
**No  sediment  correction  necessary;  in  the  case  of  a  limestone  site  of  unknown 
hardness,  no  correction  assumed. 


Table  A1 


Northern  Novaya  Zemlya  Explosions 
P  Wave  Amplitude 
Nanometers 


STATION 

10/27/66 

10/21/67 

11/07/65 

10/14/69 

10/14/70 

(••vest) 

BKS 

506 

110 

176 

275 

2113 

COR 

374 

132 

220 

264 

LON 

165 

247 

CSC 

57 

159 

MSO 

DL'C 

355 

66 

115 

166 

1100 

TIC 

134 

195 

235 

BOZ 

ALQ 

187 

57 

82 

121 

721 

COL 

225 

57 

79 

125 

RCD 

1145 

LIB 

209 

242 

418 

JCT 

166 

233 

DAL 

635 

FLO 

2267 

OXF 

990 

203 

305 

451 

2373 

SHA 

1321 

305 

396 

616 

3699 

AAM 

1453 

330 

418 

506 

2654 

ATL 

652 

137 

2055 

BLA 

269 

412 

SCP 

192 

305 

440 

CCO 

132 

256 

415 

2267 

OCD 

159 

214 

291 

WES 

132 

192 

275 

1541 

BEC 

305 

1 10 

(cast) 

( 


Table  A1 


Northern  Novaya  Zemlya  Explosions  (continued) 
P  Wave  Amplitude 
Nanometers 


STATION 

9/27/71 

8/28/72 

9/12/73 

8/29/74 

8/23/75 

10/21/’ 

(west) 

BKS 

B03 

2752 

616 

605 

990 

COR 

418 

572 

528 

B80 

LON 

412 

638 

682 

781 

CSC 

605 

MSO 

561 

737 

DUG 

379 

339 

TUC 

726 

407 

BOZ 

ALQ 

203 

242 

291 

COL 

RCD 

990 

LUB 

2763 

JCT 

DAL 

2466 

1365 

FLO 

OXF 

1425 

3313 

1001 

1012 

1420 

SHA 

836 

5064 

1 1B9 

1497 

1717 

AAM 

2069 

1100 

ATL 

946 

550 

BLA 

SCP 

CEO 

770 

1343 

OCD 

933 

693 

605 

749 

V.'ES 

990 

572 

13  nc 

517 

319 

385 

319 

(vast) 

Table  A2 


Southern  Novaya  Zemlya  Explosions 
P  Wave  Amplitude 
Nanometers 


STATION 

10/27/73 

11/02/74 

10/18/75 

(•west) 

BKS 

220 

2730 

968 

COR 

176 

1078 

LON 

267 

GSC 

203 

MSO 

DUG 

119 

TUC 

121 

BOZ 

ALQ 

81 

GOL 

44 

481 

ROD 

LUE 

242 

JCT 

754 

DAL 

FLO 

OXF 

255 

3071 

1315 

SKA 

3434 

1233 

AAM 

242 

4051 

ATL 

1387 

BLA 

123 

2036 

SCP 

214 

CEO 

OGD 

88 

561 

WES 

750 

BEG 

(east) 

T 


Table  A3 


Semipalitinsk  East  Explosions 
P  Wave  Amplitude 
Nanometers 


STATION 

1/15/65 

11/30/69 

1 1/02/72 

7/23/73 

12/14  '73 

(west) 

BKS 

66 

55 

126 

142 

55 

COR 

142 

276 

315 

LON 

71 

146 

63 

CSC 

52 

60 

MSO 

33 

DUG 

TUC 

11 

7 

15 

17 

4}£ 

BOZ 

66 

ALQ 

19 

15 

36 

10 

COL 

55 

42 

BO 

31 

RCD 

86 

LUB 

31 

63 

63 

JCT 

4 

6 

10 

DAL 

78 

116 

FLO 

35 

31 

OXF 

39 

102 

150 

SHA 

AAM 

118 

63 

165 

244 

110 

ATL 

31 

23 

59 

BLA 

23 

19 

55 

74 

SCP 

55 

102 

CEO 

71 

67 

OCD 

96 

130 

39 

WES 

43 

55 

BEC 

63 

(east) 

l 


Table  3 


Semipalitinsk  East  Explosions  (continued) 
P  Wave  Amplitude 
Nanometers 


STATION 

5/31/74 

7/04/76 

(west) 

BKS 

39 

COR 

7B 

66 

LON 

52 

55 

CSC 

26 

23 

MSO 

22 

DUG 

100 

TUC 

5 

BOZ 

ALQ 

11 

COL 

34 

RCD 

LUB 

JCT 

7 

DAL 

FLO 

OXF 

SKA 

AAM 

47 

ATL 

23 

BLA 

17 

19 

SCP 

GEO 

47 

67 

OGD 

39 

27 

WES 

23 

15 

BEC 

(cast) 

Table  A4 


Additional  Semipalitinsk  East  Explosions 
P  Wave  Amplitude 
Nanometers 


STATION 

11/23/76 

12/07/76 

5/29/77 

(west) 

BKS 

61 

55 

COR 

178 

111 

132 

LON 

CSC 

33 

33 

33 

0B2NV 

63 

MS  0 

SO 

44 

35 

DUG 

147 

123 

TUC 

B 

7 

BOZ 

ALQ 

13 

16 

15 

GOL 

44 

45 

LUB 

JCT 

DAL 

RKON 

231 

326 

273 

OXF 

SHA 

AAM 

ATL 

33 

33 

BLA 

44 

41 

22 

SCP 

77 

50 

50 

CEO 

no 

OCD 

44 

50 

49 

WES 

HNME 

132 

112 

(east) 

Table  A5 


Semipalitinsk  West  Explosions 
P  Wave  Amplitude 
Nanometers 


STATION 

10/19/66 

4/20/67 

10/17/67 

9/29/69 

6/23/70 

(west) 

BKS 

31 

27 

15 

39 

47 

COR 

7B 

110 

86 

LON 

53 

31 

37 

65 

51 

CSC 

14 

14 

5 

MSO 

DUG 

63 

TUC 

3  X 

1 

4 

6 

BOZ 

39 

42 

ALQ 

B 

6 

7 

13 

13 

COL 

22 

16 

15 

37 

34 

RCD 

31 

LUB 

15 

11 

JCT 

4H 

3 

3H 

DAL 

31 

FLO 

OXF 

29 

17 

35 

SHA 

A  AM 

31 

23 

39 

47 

ATL 

13 

4 

21 

21 

BLA 

20 

SCP 

11 

11 

6 

19 

27 

CEO 

OCD 

11 

WES 

BEC 

(uasl) 

Table  A5 


Semipatitinsk  West  Explosions  (continued) 
P  Wave  Amplitude 


Nanometers 


STATION 

3/22/71 

4/25/71 

12/30/71 

2/20/7: 

(west) 

BKS 

31 

63 

43 

COR 

110 

LON 

43 

7B 

65 

CSC 

17 

MSO 

23 

DUG 

TUC 

B 

3H 

B02 

ALQ 

B 

B 

COL 

24 

53 

27 

25 

RCD 

LUB 

15 

39 

IB 

jcr 

5 

13 

4 

3 

DAL 

55 

31 

FLO 

35 

OXF 

31 

35 

21 

SKA 

AAM 

39 

47 

31 

51 

ATL 

27 

15 

BLA 

19 

13 

B 

SCP 

19 

43 

CLO 

ocn 

WES 

BHC 

11 

(east) 

Table  A6 


Western  Kazakh  Explosions 
P  Yt'ave  Amplitude 
Nanometers 


STATION  12/06/69  12/12/70  12/23/70 


(west) 

BKS 

17 

COR 

63 

LON 

35 

CSC 

3H 

MSO 

DUG 

19 

TUC 

i  H 

BOZ 

ALQ 

10 

COL 

21 

RCD 

LUB 

JCT 

DAL 

FLO 

33 

OXF 

83 

SKA 

AAM 

77 

ATL 

24 

BLA 

30 

SCP 

110 

CEO 

70 

OCD 

22 

WES 

I3EC 

(e.isl) 

2S 

21 

70 

71 

60 

55 

30 

25 

5 

6 

17 

15 

32 

39 

33 

25 

7 

11 

40 

150 

123 

78 

78 

35 

39 

<36 

79 

175 

213 

150 

157 

60 

89 

49 

59 

Table  A7 


Kurile  Earthquakes 
P  Wave  Amplitudes 

Corrected  for  Geometric  Spreading  to  A  =  60° 
Nanometers 


STATION 

11/22/66 

3/20/67 

8/10/67 

2/10/63 

4/23/69 

(v.-est) 

BKS 

112 

101 

100 

COR 

453 

253 

LON 

161 

24 

GSC 

0B2NV 

417 

66 

67 

MSO 

DUG 

513 

79 

125 

14B 

TUC 

340 

36 

59 

71 

BOZ 

509 

64 

143 

238 

ALQ 

299 

25 

65 

52 

GOL 

202 

18 

46 

53 

15 

Ll’B 

450 

120 

45 

JCT 

440 

60 

DAL 

744 

95 

RKON 

OXF 

973 

161 

184 

319 

SHA 

496 

AAM 

475 

197 

305 

35 

ATL 

462 

76 

154 

182 

55 

BLA 

407 

83 

174 

262 

65 

SCP 

351 

44 

111 

139 

57 

CEO 

27] 

OGD 

222 

WES 

HNMC 

(east) 

435 

104 

191 

191 

Table  A  7 


Kurile  Earthquakes  (continued) 

P  Wave  Amplitude 

Corrected  for  Geometric  Spreading  to  £  =  60° 
Nanometers 


STATION* 

7/25/63 

10/25/76 

3/19/77 

(west) 

BKS 

133 

912 

COR 

242 

13S8 

LON 

195 

CSC 

710 

924 

0B2NV 

MSO 

B42 

127B 

DUG 

165 

1179 

1197 

TUC 

122 

505 

724 

BOZ 

ALQ 

149 

• 

GOL 

233 

LUB 

251 

1545 

JCT 

DAL 

1049 

RKON 

OXF 

SKA 

AAM 

4634 

ATL 

765 

2174 

BLA 

146 

2715 

SCP 

CCO 

1404 

OGD 

273 

WES 

1 62 

HNME 

(east) 

Table  AB 


Japanese  Earthquakes 
P  Wave  Amplitude 

Corrected  for  Geometric  Spreading  to  A  =  60“ 
Nanometers 


STATION  12/31/76  1/01/77  2/18/77  6/12/77 


(west) 


BKS 

196 

734 

COR 

1B6 

323 

LON 

116 

CSC 

7 

56 

325 

19 

0B2XV 

16 

39 

MSO 

18 

53 

310 

30 

DUG 

19 

93 

4B3 

51 

TUC 

8 

16 

15 

B02 

ALQ 

12 

37 

30 

COL 

9 

LUB 

JCT 

DAL 

RKON 

44 

50 

OXF 
SKA 
A  AM 
ATL 
BLA 
SCP 
CEO 
OGD 
WES 
KNME 
(cast) 


Table  A9 


Bonin  Islands  Earthquakes 
P  Wave  Amplitude 

Corrected  for  Geometric  Spreading  to  A  »  60° 
Nanometers 

STATION  9/22/76  12/05/76  12/22/76  1/05  '77 


(west) 


BKS 

212 

99 

339 

COR 

106 

LON 

43 

37 

CSC 

149 

79 

157 

0E2.NV 

210 

MSO 

33 

B7 

DUG 

277 

115 

124 

TUC 

51 

22 

132 

60 

BOZ 

ALQ 

26 

50 

149 

COL 

52 

12 

62 

LUB 

170 

JCT 

DAL 

RKON 

179 

315 

OXF 

SKA 

AAM 

ATL 

BU 

SCP 

GEO 

OGD 

IVES 

H.N.V.E 

(east) 


Table  AID 

Other  Earthquakes  to  a  Northwest  Azimuth  from  the  United  States 

P  Wave  Amplitude 

Corrected  for  Geometric  Spreading  to  A  =  50° 
Nanometers 


STATION 

10/22/75 

4/22/77 

4/23/77 

7/20/77 

S/C7/77 

(west) 

EKS 

44 

53 

85 

CCR 

112 

Lon 

21 

CSC 

44 

0B2NV 

52 

26 

63 

54 

MSO 

65 

14 

26 

105 

DUG 

45 

70 

21 

TUC 

55 

22 

7 

44 

47 

BOZ 

ALQ 

31 

GOL 

20 

13 

LUB 

JCT 

56 

45 

DAL 

97 

RKOX 

30 

146 

OXF 

SKA 

AAM 

ATL 

35 

BLA 

130 

SCP 

14 

1C5 

GEO 

50 

OGD 

49 

37 

WES 

HNME 

74 

(east) 

..'*1 


Table  All 


Soulh  American  Earthquakes  (continued) 

P  Wave  Amplitude 

Corrected  for  Geometric  Spreading  to  A  =  60° 
Nanometers 


STATION 

9/30/76 

12/03/76 

12/04/76 

3/09/77 

3/13/77 

(west) 

BKS 

120 

B6 

107 

COR 

214 

186 

143 

LON 

59 

GSC 

165 

95 

126 

29 

59 

0B2NV 

MSO 

100 

49 

75 

62 

DUG 

150 

45 

121 

55 

53 

TUC 

127 

61 

70 

20 

BOZ 

ALQ 

54 

42 

GOL 

41 

48 

LUB 

103 

JCT 

27 

DAL 

257 

RKON 

96 

OXF 

SHA 

AAM 

ATL 

59 

37 

65 

54 

BLA 

60 

76 

SCP 

CEO 

96 

OGD 

34 

13 

WES 

HNME 

(east) 

37 

Table  A1 1 


South  American  Earthquakes  (continued) 

P  Y/ave  Amplitude 

Corrected  for  Geometric  Spreading  to  A  *  80* 
Nanometers 


STATION 

4/15/77 

6/02/77 

6/05/77 

6/09/77 

6/19/77 

(v.-est) 

BX5 

177 

62 

151 

203 

COR 

119 

177 

LON 

32 

11 

26 

363 

CSC 

29 

66 

161 

215 

OB2NV 

131 

56 

74 

156 

MSO 

93 

13 

99 

116 

DUG 

33 

69 

176 

225 

TL'C 

116 

17 

63 

149 

130 

BOZ 

ALQ 

GOL 

Ll’B 

JCT 

DAL  67 

RKOX  221  60  171 

OaF 
SKA 
AAM 


ATL 

72 

75 

96 

BLA 

83 

69 

264 

SCP 

CEO 

72 

69 

77 

OGD 

40 

34 

YiES 

KNME 

125 

65 

(cast) 


Table  A.12 


rm 


Northern  Novaya  2emiya 


STATION 

MEAN 

S.E.M. 

N 

(west) 

BKS 

0.B9 

0.C4 

10 

COR 

0.52 

0.04 

B 

LON 

0.52 

0.04 

6 

CSC 

0.65 

0.13 

3 

MSO 

0.77 

0.03 

2 

DUG 

0.53 

0.02 

7 

TUC 

0.51 

0.04 

5 

BOZ 

ALQ 

0.37 

0.01 

B 

GOL 

0.35 

0.01 

4 

RCD 

1.84 

0.00 

2 

LUB 

1.23 

0.10 

4 

JCT 

1.10 

0.06 

2 

DAL 

2.63 

0.11 

3 

FLO 

1.15 

1 

ONF 

1.43 

0.03 

10 

SHA 

1.B9 

0.07 

10 

AAM 

1.95 

0.13 

7 

ATL 

1.02 

0.03 

5 

BLA 

1.86 

0.02 

2 

SCP 

1.37 

0.01 

3 

GEO 

1.19 

0.07 

6 

OGD 

0.92 

0.04 

7 

VfES 

0.92 

0.04 

6 

BEC 

0.51 

0.02 

6 

(cast) 

Table  A13 


Southern  Novaya  Zemlya 


STATION 

MEAN 

S.E.M. 

N 

(west) 

BKS 

1.39 

0.09 

3 

COR 

1.31 

0.09 

2 

LON 

1.94 

1 

CSC 

1.41 

1 

MSO 

DUG 

0.93 

1 

TUC 

0.94 

1 

BOZ 

ALQ 

0.56 

1 

GOL 

0.27 

0.03 

2 

RCD 

LUB 

1.67 

1 

JCT 

0.99 

1 

DAL 

FLO 

OXF 

1.67 

0.07 

3 

SKA 

1.66 

0.06 

2 

AAM 

1.95 

0.19 

2 

ATL 

0.69 

1 

BLA 

0.94 

0.03 

2 

SCP 

1.46 

1 

CEO 

OGD 

0.67 

0.06 

2 

WES 

0.33 

1 

BEC 

(cast) 

*7 


Table  A15 


Semipalatinsk  East  Including  Additional  Events 


STATION 

MEAN 

S.E.M. 

N 

(west) 

BKS 

1.26 

0.05 

8 

COR 

3.26 

0.56 

6 

LON 

1.52 

0.06 

5 

CSC 

0.77 

0.03 

7 

OB2NV 

1.34 

1 

MSO 

0.62 

0.05 

5 

DUG 

2.63 

0.12 

3 

TUC 

0.15 

0.01 

B 

BOZ 

1.47 

1 

ALQ 

0.32 

0.02 

6 

GOL 

0.69 

0.04 

7 

LUB 

0.67 

0.05 

3 

JCT 

0.12 

0.03 

4 

DAL 

1.21 

0.13 

2 

RKON 

5.75 

0.37 

3 

OXF 

1.14 

0.13 

3 

SKA 

AAM 

1.90 

0.19 

6 

ATL 

0.62 

0.03 

6 

BLA 

0.59 

0.05 

9 

SCP 

1.16 

0.09 

5 

GEO 

1.69 

0.21 

5 

OGD 

0.99 

0.05 

6 

WES 

0.53 

0.03 

4 

KNME 

2.38 

0.02 

2 

(east) 

<**?%*i 


Table  A16 


Semipalitinsk  IVest 


STATION' 

MEAN 

S.E..M. 

N 

(west) 

BKS 

1.55 

0.06 

8 

COR 

3.74 

0.25 

4 

LON 

2.33 

0.22 

8 

CSC 

0.76 

0.10 

4 

MSO 

1.15 

1 

DUG 

8.25 

1 

TUC 

0.16 

0.01 

7 

BOZ 

3.39 

0.S2 

2 

ALQ 

0.47 

0.C4 

7 

COL 

1.20 

0.05 

9 

RCD 

1.51 

1 

LUB 

0.S1 

0.04 

5 

JCT 

0.22 

0.03 

7 

DAL 

1.19 

0.C5 

3 

FLO 

0.30 

1 

OXF 

1.33 

0.10 

6 

SKA 

A  AM 

1.53 

0.16 

6 

ATL 

0.63 

0.05 

6 

BLA 

0.52 

0.07 

4 

SCP 

0.75 

0.07 

7 

CEO 

OCD 

0.40 

1 

WES 

BEC 

(cast) 

0.57 

1 

Table  A17 


Western  Kazakh 


STATION 

MEAN 

S.E.M. 

N 

(west) 

BKS 

0.65 

0.07 

3 

COR 

2.05 

0.34 

3 

LON 

1.44 

0.06 

3 

CSC 

0.15 

1 

MSO 

DUG 

0.73 

0.C7 

3 

TUC 

0.11 

0.02 

3 

BOZ 

ALQ 

0.41 

0.02 

3 

COL 

O.BS 

0.05 

3 

RCD 

LUB 

0.71 

0.11 

2 

JCT 

022 

C.04 

2 

DAL 

0.99 

1 

FLO 

1.44 

1 

OXF 

3.42 

0.24 

3 

SKA 

AAM 

2.35 

0.47 

3 

ATL 

0.95 

0.05 

3 

BLA 

1.62 

0.10 

3 

SCP 

4.72 

0.23 

3 

GEO 

3.50 

0.23 

3 

OGD 

1.52 

0.33 

3 

WES 

1.10 

1 

BEC 

1.40 

1 

(east) 

Other  Northwest  Data 


STATION 

MEAN 

S.E.M. 

N 

(west) 

BKS 

0.95 

0.16 

3 

COR 

1.44 

1 

LON 

1.39 

1 

CSC 

0.7? 

1 

032NV 

1.16 

0.21 

4 

MSO 

1.04 

0.22 

4 

DUG 

1.17 

0.27 

3 

TUC 

0.62 

0.07 

5 

BOZ 

ALQ 

0.44 

1 

GOL 

0.26 

0.00 

2 

LUB 

JCT 

2.07 

0.57 

2 

DAL 

1.24 

1 

RKON 

2.23 

0.29 

2 

OXF 

SHA 

AAM 

ATL 

0.50 

1 

BLA 

2.25 

1 

SCP 

1.39 

0.44 

2 

GEO 

0.90 

1 

OGD 

0.69 

0.03 

2 

VES 

HNME 

1.29 

1 

(east) 

;>f  fyar0*,i 


